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^ Abstract 

r— I The quest to comprehend how nuclear processes influence astrophysical phenomena 
is driving experimental and theoretical research programs worldwide. One of the 

I— I main goals in nuclear astrophysics is to understand how energy is generated in 
^ stars, how elements are synthesized in stellar events and what the nature of neutron 
stars is. New experimental capabilities, the availability of radioactive beams and 
increased computational power paired with new astronomical observations have 
advanced the present knowledge. This review summarizes the progress in the field 

^ of nuclear astrophysics with a focus on the role of indirect methods and reactions 
involving beams of rare isotopes. 
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1 Introduction 



Nuclear reactions in stars and stellar explosions are responsible for the ongoing synthesis of 
the elements [1-9] . Nuclear physics plays an important role as it determines the signatures of 
isotopic and elemental abundances found in the spectra of stars, novae, supcrnovae, and X-ray 
burst as well as in characteristic 7-ray radiation from nuclear decays, or in the composition 
of meteorites and prcsolar grains (see [10] for a review on the nuclear structure input to 
nuclear astrophysics). 

The rapid neutron capture process (r process) is responsible for the existence of about half 
of the stable nuclei heavier than iron; yet a site that can produce the observed elements 
self-consistently has not been identified [11,12]. Capture cross sections for most of the nuclei 
involved are hard if not impossible to measure in the laboratory and indirect experimental 
approaches have to be employed to gather the relevant nuclear structure information. Nu- 
clear masses and /3-decay half-lives are among the few direct observables that are input for 
calculations that model nucleosynthesis in the r process. X-ray bursts provide a unique win- 
dow into the physics of neutron stars. They arc the most frequent thermonuclear explosions 
known. The brightness, frequency and opportunity to be observed with different telescopes 
makes them unique laboratories for explosive nuclear burning at extreme temperatures and 
densities [13-15]. The reaction sequence during an X-ray burst proceeds through nuclei at 
or close to the proton drip line mainly by (p, 7) and (p, a) reactions and /3 decays (ap and 
rp process) [16] . Most rp-process reaction rates are still based exclusively on theory. Energy 
in explosive hydrogen burning events such as X-ray bursts is initially generated in the CNO 
cycle and as the temperature increases, ct-capture on unstable oxygen and neon (^^O and 
^^Ne) leads to a break out and an ensuing chain of proton captures that can go as far as tin. 
Supernovae play a crucial role in the understanding of the universe as they are the major 
source of nucleosynthesis and possibly of cosmic rays. Core-collapse supernovae [17] are one 
of the proposed sites of the r process. Thermonuclear supernovae (type la) are powered by 
explosive carbon and oxygen burning of a white dwarf that has reached the Chandrasekhar 
mass limit. For both types of supernovae, the driving processes are not well understood and 
weak interaction rates play a key role [18]. Temperatures and densities are so high that 
electron captures on unstable nuclei become crucial. 

Most aspects in the study of nuclear physics demand beams of energetic particles to induce 
nuclear reactions on the nuclei of target atoms. It was from this need that accelerators were 
developed. Over the years, many ways of accelerating charged particles to ever increasing 
energies have been devised. Today we have ion beams of all elements from protons to uranium 
available at energies well beyond those needed for the study of atomic nuclei. The quantities 
used in nucleosynthesis calculations are reaction rates. A thermonuclear reaction rate is a 
function of the density of the interacting nuclei, their relative velocity and the reaction cross 
section. Extrapolation procedures are often needed to derive cross sections in the energy 
or temperature region of astrophysical relevance. While non-resonant cross sections can be 
rather well extrapolated to the low-energy region, the presence of continuum, or sub-threshold 
resonances, can complicate these extrapolations. We will mention some of the important 



1 



examples. 

In the Sun, the reaction ^Be(p,7)'^B plays a major role for the production of high energy 
neutrinos from the /3-decay of ^B. These neutrinos come directly from the center of the Sun 
and are ideal probes of the Sun's structure. John Bahcall frequently said that this was the 
most important reaction in nuclear astrophysics [19]. Our understanding of this reaction 
has improved considerably with the advent of rare-isotope beam facilities. The reaction 
^^C(a, 7)^^0 is extremely relevant for the fate of massive stars. It determines if the remnant 
of a supernova explosion becomes a black- hole or a neutron star [20]. These two reactions 
are only two examples of a large number of reactions, which are not yet known with the 
accuracy needed in astrophysics. 

In this review, we summarize recent developments and achievements in nuclear astrophysics 
with a focus on theoretical approaches and experimental techniques that are applicable 
to or utilize rarc-isotope beams, respectively. Section 2 will cover reactions within stars. 
Section 3 is devoted to nuclear reaction models, Section 4 reviews the effect of environment 
electrons. Section 5 outlines approaches with indirect methods and Section 6 summarizes 
recent nuclear astrophysics experiments with rare-isotope beams. Finally, in Section 7 we 
present our outlook for the present and future of this field. 



2 Reactions within stars 



2.1 Thermonuclear cross sections and reaction rates 



The nuclear cross section for a reaction between a nuclear target j and a nuclear projectile 
k is defined by 

number of reactions target~^sec~^ r/n^ 



flux of incoming projectiles rikV ' 

where the target number density is given by rij, the projectile number density is given by n^, 
and V is the relative velocity between target and projectile nuclei. The number of reactions 
per unit volume and time can be expressed as r = avrijUk, or, more generally, by 

J cr\vj-Vk\(fnj(fnk. (2) 



The evaluation of this integral depends on the type of particles and their distributions. For 
nuclei j and k in an astrophysical plasma, obeying a Maxwell-Boltzmann distribution (MB), 

Eq. (2) simplifies to rj^k = {o'v) rijUk, where the reaction rate (av) is the average of av over 
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the temperature distribution in (3). More specifically, 

= (c^v)^-,fc"'i"'fc (4) 

Here rrijk denotes the reduced mass of the target-projectile system. 
2.1.1 Photons 

When in Eq. (2) particle A; is a photon, the relative velocity is always c and there is no need 

to integrate quantities over d^rij. Thus, one obtains rj = Aj-.^n^ where A^^^ results from an 
integration of the photodisintegration cross section over a Planck distribution for photons of 
temperature T 

E'^ 1 

d^n^ — , ^ ^, dE^, (6) 

^ 7r^{ch)^exp{E^/kT) -1 ^' 

which leads to 



r 



1 r ^ ca{E^)E^ 

— Ao-vfTjn,- = d rij , ^ dE^. (7) 

^'^^ ' ^ TT^ichYJ ^ Jo exp E^/m-l ^ 



^io exp(E^/A;T) - 1 



There is, however, no direct need to evaluate photodisintegration cross sections, because, 
due to detailed balance, they can be expressed by the capture cross sections for the inverse 
reaction Z + m — > j + 7 [21] 

where m„ = mi2Q/12 is the mass unit, Q/,„ is the reaction Q-value, T is the temperature, 
(o"v)^. ^ is the inverse reaction rate, ^{T) = J2i{2Ji + 1) cxp{—Ei/kT) are partition functions, 
and A are the mass numbers of the participating nuclei in a thermal bath of temperature T. 



2.1.2 Electron, positron and neutrino capture 

The electron is about 2000 times less massive than a nucleon. Thus, the velocity of the 
nucleus j is negligible in the center of mass system in comparison to the electron velocity 
{\vj — Ve\ ~ |fe|), and there is no need to integrate quantities over d^rij. The electron capture 
cross section has to be integrated either over a Boltzmann or a Fermi distribution of electrons, 
depending on the astrophysical scenario. The electron capture rates are a function of T and 
the electron number density, rig = YepN^ [22]. In a completely ionized plasma, = J2i ZiYi, 
i.e., the electron abundance is equal to the total proton abundance in nuclei. Here denotes 
the abundance of nucleus i defined by Yi = ni/{pNA), where Ui is the number density of 



3 



nuclei per unit volume and Na is Avogadro's number. Therefore, 



(9) 



This treatment can be generalized for the capture of positrons, which are in a thermal 
equilibrium with photons, electrons, and nuclei. At high densities (p > 10^^ g/cm^) the 
neutrino scattering cross sections on nuclei and electrons are large enough to thermalize the 
neutrino distribution. Inverse electron (neutrino) capture can also occur and the neutrino 
capture rate can be expressed similarly to Eqs. (7) or (9), integrating over the neutrino 
distribution. 



2.1.3 Beta- decay 

For normal decays, like /? or a decays with half-life ri/2, we obtain an equation similar to 
Eqs. (7) or (9) with a decay constant Aj — ln2/ri/2 and 

Tj = XjUj. (10) 



2.1.4 Charged particles 

Nuclear cross sections for charged particles are suppressed at low energies due to the Coulomb 
barrier. The product of the penetration factor and the Maxwell-Boltzmann (MB) distribution 
at a given temperature yields an energy window in which most of the reactions occur, known 
as the Gamow window. 



Experimentally, it is more convenient to work with the astrophysical S factor 

S{E) = Ea{E) exp(27r77), (11) 

with 77 being the Sommerfeld parameter, describing the s-wave barrier penetration, rj — 
ZjZke^/hv, and energy E and v the relative velocity of the ions. In this case, the steep 
increase of the cross section is transformed into a rather flat, energy dependent function. 
One can easily see the two contributions of the velocity distribution and the penetrability in 
the integral 

1/2 

(12) 



kT £;V2 



where b = 27rr]E^^^ — {2mjky^'^7re'^ZjZk/h and mjk is the reduced mass in units of m„ 
(unit atomic mass). Experimentally it is very difficult to perform direct measurements of 
fusion reactions involving charged particles at very small energies. The experimental data at 
higher energies can be guided by a theoretical model for the cross section, which can then 
be extrapolated down to the Gamow energy. However, the extrapolation can be inadequate 
due to the presence of resonances and subthreshold resonances, for example. 
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Fig. 1. Schematic representation of the energy dependence of a fusion reaction involving charged 
particles (Courtesy of C. Spitaleri). 

A simple result can be obtained by assuming a constant 5'-factor, i.e., S{E) = SEq. In this 
case, the first derivative of the integrand in Eq. (12) yields the location Eq of the Gamow 
peak, and the effective width A of the energy window, i.e. 

A = (^l^^y^' = 0.749(Z|Z>,,T|)V6keV, (13) 

carrying the dependence on the charges Zj, Zk, the reduced mass rrijk of the involved nuclei 
in units of rriu, and the temperature Tg given in 10^ K. 



Figure 1 outlines one of the main challenges in astrophysical reactions with charged particles. 
The experimental data can be guided by a theoretical model for the cross section, which can 
then be extrapolated to the Gamow energy. The solid curve is a theoretical prediction, which 
supposedly describes the data at high energies. Its extrapolation to lower energies yields the 
desired value of the S-factor, or cross section, at the Gamow energy £"0- The extrapolation 
can be complicated by the presence of unknown resonances. 



2.1.5 Neutron-induced reactions 

For neutron-induced reactions, the effective energy window for s-wave neutrons (/ = 0) is 
given by the location and width of the peak of the MB distribution function. For / > 0, the 
penetrability of the centrifugal barrier shifts the effective energy Eq to higher values. For 
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neutrons with energies less than the height of the centrifugal barrier one gets [23] 

Eo ^ 0.172r9 (^l + MeV, A ^ 0.19479 + ^) ^ MeV, (14) 

Usually, £^0 is not much different (in magnitude) from the neutron separation energy. 

2.2 Reaction networks 



The time evolution of the number densities, rii, of each of the species i in an astrophysical 
plasma (at constant density) is obtained by solving equations of the type 

V / p=const j j,k j,k,l 

where the N^. can be positive or negative numbers that specify how many particles of species 
i are created or destroyed in a reaction x. The reactions x fall in three categories: 

(1) decays, photodisintegrations, electron and positron captures and neutrino-induced re- 
actions, rj — XjUj, 

(2) two-particle reactions, rj^k — ("^v)^. n^nfe, and 

(3) three-particle reactions, rjki — i k like the triple-o; process {a + a + a — > 
The A^*'s are given by: 

Ni Ni 
Nj = Ni, Nj^^ = - and Nj^f.^i = -— — — -, 

1 lm=l I ^ ^Jm I • 1 lm=l \-^^3m\- 

where the products in the denominators run over the different species destroyed in the 
reaction and avoid double counting when identical particles react with each other. 

In terms of the nuclear abundances, Yi — ni/{pNA) such that for a nucleus with atomic 
weight Ai, AjXi represents the mass fraction of this nucleus, Y.-^iXi — 1 ^-nd the reaction 
network equations can be rewritten as 

dY 

^ = E N]^jY3 + E ^IkPNA (av)^. , Y,Y, + ^ Nl,^,p'Nl (av)^. Y,Y,Y,. (16) 

3 3,k j,k,l 

The energy generation per unit volume in a time interval At is expressed in terms of the 
mass excess AMjC^ of the participating nuclei 

Ae = -^Ay,7V^AM,c^ ^ = - ^ ^TV^AM.c^ (17) 
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Fig. 2. Example of reaction networks (pp-chains and CNO-cycles). A particular nucleus on the 
Segre chart can take different paths along the reaction network, as shown in the inset at the lower 
right side. (Courtesy of S. Typel) 

The solution of the above group of equations allows to deduce the path for the r process 
until the heavier elements are reached. The relative abundances of elements are also obtained 
theoretically by means of these equations by using stellar models for the initial conditions, 
as the neutron density and the temperature. Nuclear physics has to contribute with /5-decay 
half-lives, electron and positron capture rates, photo-nuclear and neutrino cross sections. 

Simple examples of reaction networks are shown in Figure 2 for typical pp-chains and CNO- 
cycles. On the right we show a particular nucleus on the Segre chart from where different 
paths can start along the reaction network. 



3 Nuclear reaction models 



Explosive nuclear burning in astrophysical environments produces short-lived, exotic nuclei, 
which again can be targets for subsequent reactions. In addition, it involves a very large num- 
ber of stable nuclei, which are still not fully explored by experiments. Thus, it is necessary to 
be able to predict reaction cross sections and thermonuclear rates with the aid of theoretical 
models. Especially during the hydrostatic burning stages of stars, charged-particle induced 
reactions proceed at such low energies that a direct cross-section measurement is often not 
possible with existing experimental techniques. Hence extrapolations down to the stellar en- 
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ergies of the cross sections measured at the lowest possible energies in the laboratory are 
usually applied. To be trustworthy, such extrapolations should have as strong of a theoreti- 
cal foundation as possible. Theory is even more mandatory when excited nuclei are involved 
in the entrance channel, or when unstable, very neutron-rich or neutron-deficient nuclides 
(many of them being even impossible to produce with present-day experimental techniques) 
have to be considered. Such situations are often encountered in the modelling of explosive 
astrophysical scenarios. 

3.1 Potential and DWBA models 

Potential models assume that the physically important degrees of freedom arc the relative 
motion between structureless nuclei in the entrance and exit channels. The only microscopic 
information is introduced in terms of spectroscopic factors and parameters of the optical 
potential. The weakness of the models is that the nucleus-nucleus potentials adopted for 
calculating the initial and final wavefunctions from the Schrodinger equation cannot be 
unambiguously defined. Single-particle wavefunctions are calculated using nuclear potentials 
of the form 



where V^)(r) and Vs{r) are the central and spin-orbit interactions, respectively, and Vc{r) 
is the Coulomb potential of a uniform distribution of charges. The potentials Vo(r) and 
Vs{t)^ are usually given in terms of a Woods-Saxon (WS) parameterization. The parameters 
of the potentials - their depth, range and difi^useness, are chosen to reproduce the ground 
state energy Eb (or the energy of an excited state). For knockout reactions, they are also 
adjusted to reproduce the orbital radius of the nucleon. Most often, the same parameters 
do not reproduce the proper continuum wavefunctions, do not yield location and widths 
of resonances, etc. These can be obtained by readjusting the strengths of the potentials, 
effectively increasing the number of parameters at hand. 

The WS parameterization is well suited to describe any reaction of interest, except perhaps 
for those cases in which one of the partners is a neutron-rich halo nucleus. Then, the extended 
radial dependence leads to unusual forms for the potentials. Also, for capture reactions in 
which the light partner is either a deuteron, triton, a-particle or a heavier nucleus, folding 
models are more appropriate. The central part of the potential is obtained by a folding of 
an effective interaction with the ground state densities, pa and pb, of the nuclei A and B: 



with s = |r -|- r2 — ri|. Aq is a normalization factor which is close to unity. 

Folding models are based on an effective nucleon- nucleon interaction, t'e//, and nuclear den- 
sities, pi, which are either obtained experimentally (not really, because only charge densities 
can be accurately determined from electron scattering), or calculated from some microscopic 
model (typically Hartree-Fock or relativistic mean field models). The effective interactions 



V{v)^V^{r) + Vs{T) (l.s) + Fc(r) 



(18) 




(19) 
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as well as the nuclear densities are subject of intensive theoretical studies. 



Potential models have been applied to all kinds of calculations for nuclear astrophysics. For 
simplicity, let us consider radiative capture reactions involving a target nucleus and a nucleon. 
The wavefunctions for the nucleon (n) + nucleus (x) system are calculated by solving the 
radial Schrodinger equation 



2m,- 



_/(/ + !) 
dr^ r 



2 



Ua (r) + V{r)ua (r) = E^Ua (r) . (20) 



The nucleon n, the nucleus x, and the n + x = a-system have intrinsic spins labeled by 
s = 1/2, Ix and J, respectively. The orbital angular momentum for the relative motion of 
n+x is described by I. Angular momenta are usually coupled as 1 + s= j and j + 1^= J, where 
J is called the channel spin. In Eq. (18), for V one uses s.l = + 1) — l{l + 1) — 3/4] /2 
and a in Eq. (20) denotes the set of quantum numbers, — {Eh,lh,jb, Jb} for the bound 
state, and etc = {Ec, Ic, jc, Jc} for the continuum states. 

The bound-state wavefunctions are normalized to unity, / dr\ua^{r)f = 1, whereas the 
continuum wavefunctions have boundary conditions at large distances given by 



m 



u^Sr - oo) = [Ht^ (r) - S^M"-^ (r)J e^'^^'^ (21) 

where Sa^ = exp [2?(5q^ (E)], with (E) and ai {E) being the nuclear and the Coulomb 
phase shifts, respectively. In Eq. (21), Hi"^^ (r) = Gi{r) ± iFi (r), where Fi and Gi are the 
regular and irregular Coulomb wavefunctions. For neutrons, the Coulomb functions reduce to 
the usual spherical Bessel functions, ji (r) and rii (r). With these definitions, the continuum 
wavefunctions are normalized as (ue'Jue^) — S {E'^ — Ec) Saa'- 



For n + X a + J and vrL (tt = E, (M) = electric (magnetic) L-pole) transitions, the cross 
sections are obtained from 

radxap. ^ ( E^x + E, \^'^+^ 2(2/, + 1) L+l 

\ he ) (24 + l)(24 + l)L[(2L + l)!!]2^^^; ^+ > 

X<;-^"^'^'^ \{lcJc\\O.L\\lbJb)\\ (22) 

Jb jb L 

where Ei, is the binding energy and (/cjc II ^ttlH hjb) is the multipole matrix element. For 
electric multipole transitions, 

{Icjc WOelW hjb) = (-l)'^+'c-.c+L-i/2^^(2L + l){2j, + 1) ^ 

v47r \l/2 -1/2 

/■oo 

X / dr r^ui,{r)uc{r), (23) 
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Fig. 3. Potential model calculation [24] for the reaction ^^0(p,7)^''F. The dotted line and the 
dashed line are for the capture to the ground state and to the first excited state respectively. The 
experimental data are from Refs. [25-27]. The dotted-dashed lines are the result of shell model 
calculations published in Ref. [28] . 



where is the effective charge, which takes into account the displacement of the center-of- 
mass, cl = ZnC {—rrin/ma)^ + Z^e {m^/ma)^ ■ In comparison with electric dipole transitions, 
the cross sections for magnetic dipole transitions are reduced by a factor of f^/c^, where v 
is the relative velocity of the n + x system. At very low energies, w -C c. Ml transitions will 
be much smaller than the electric transitions. Only in the case of sharp resonances, the M 1 
transitions play a significant role. 

The total radiative capture cross section is obtained by adding all multipolarities and final 
spins of the bound state {E = Enx), 

a^^^-'^^^-{E) = Y.{SF)j,afX{E) , (24) 
where {SF)j^ are spectroscopic factors. 

As an example. Figure 3 shows a potential model calculation [24] for the S'-factor of the 
^^0(p,7)^^F reaction. The rate of this reaction influences sensitively the ^^O/^^O isotopic 
ratio predicted by models of massive (> 4:Mq) AGB stars, where proton captures occur at 
the base of the convective envelope (hot bottom burning). A fine-tuning of the ^^0(p,7)^''F 
reaction rate may account for the measured anomalous ^''O/^^O abundance ratio in small 
grains which are formed by the condensation of the material ejected from the surface of AGB 
stars via strong stellar winds [29]. The agreement of the potential model calculation with 
the experimental data seen in Figure 3 is very good and comparable to more elaborated 
calculations [28]. 
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Fig. 4. Microscopic calculations for the reaction ^Be(p, 7)^6. The dashed line is the no-core shell- 
model calculation of Ref. [33] and the dotted line is from the resonant group method calculation of 
Ref. [34]. Experimental data are from Refs. [35-37,32,38-41]. 



3.2 Microscopic models 



In microscopic models, nucleons are grouped into clusters and the completely antisym- 
metrized relative wavefunctions between the various clusters are determined by solving the 
Schrodinger equation for a many-body Hamiltonian with an effective nucleon-nucleon inter- 
action. Typical cluster models are based on the Resonating Group Method (RGM) or the 
Generator Coordinate Method (GCM). They are based on a set of coupled integro-differential 
equations of the form 

E / d'r' [7^i^(r, r') - r')] 9air') = 0, (25) 

a' 

where 

7^i^(r,r') = (vl/^(a,r)|/7|vl/5(a',r')) and MtS{v,r') = A{a,r)\^ sia' ,r')) . 

In these equations H is the Hamiltonian for the system of two nuclei (A and B) with the 
energy E, '^a,b is the wavefunction of nucleus A (and B), and gai^) is a function to be found 
by numerical solution of Eq. (25), which describes the relative motion of A and B in channel 
a. Full antisymmetrization between nucleons of A and B are implicit. 

Modern nuclear shell-model calculations, such as the Monte-Carlo shell model, or the no- 
core shell model, are able to provide the wavefunctions \E'a,b for light nuclei. But so far 
they cannot describe scattering wavefunctions with a full account of anti-symmetrization. 
Moreover, the road to an effective NN interaction which can simultaneously describe bound 
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and continuum states has not been an easy one. Thus, methods based on Eq. (25) seem to 
be the best way to obtain scattering wavefunctions needed for astrophysical purposes. Old 
interactions, such as Volkov interactions, are still used for practical purposes. It is also worth 
mentioning that this approach has provided the best description of bound, resonant, and 
scattering states of nuclear systems [30] . 

As an example of applications of this method, we again give a radiative capture reaction. The 
creation and destruction of ^Be in astrophysical environments is essential for the description 
of several stellar and cosmological processes and is not well understood. also plays an 
essential role in understanding our Sun. High energy Uf. neutrinos produced by ^B decay in 
the Sun oscillate into other active species on their way to earth [31]. Precise predictions of 
the production rate of ^B solar neutrinos are important for testing solar models, and for 
limiting the allowed neutrino mixing parameters. The most uncertain reaction leading to 
^B formation in the Sun is the '''Be(p, 7)^B radiative capture reaction [32]. Additionally, the 
Coulomb dissociation method, discussed later in this review, has given some new insights 
about the electromagnetic matrix elements for this reaction. Figure 4 shows a comparison of 
microscopic calculations for the reaction ^Be(p, 7)^B with experimental data. The dashed- 
dotted line is the no-core shell-model calculation of Ref. [33] and the dotted line is for 
the resonant group method calculation of Ref. [34]. Experimental data are from Refs. [35- 
37,32,38-41]. It is evident that both, theory and experiment, need improvements for this 
important reaction. 

3.2.1 Asymptotic normalization coefficients 

Although the potential model works well for many nuclear reactions of interest in astro- 
physics, it is often necessary to pursue a more microscopic approach to reproduce exper- 
imental data. Instead of the single-particle wavefunctions one often makes use of overlap 
integrals, /^(r), and a many-body wavefunction for the relative motion, ^c(i")- Both /^(r) 
and ^c(r) might be very complicated to calculate, depending on how elaborate the micro- 
scopic model is. The variable r is the relative coordinate between the nucleon and the nucleus 
X, with all the intrinsic coordinates of the nucleons in x being integrated out. The radiative 
capture cross sections are obtained from the calculation of cr^^j^'^^''' oc | ^/^(r) | Ir^yil l^cl?")) 

The imprints of many-body effects will eventually disappear at large distances between the 
nucleon and the nucleus. One thus expects that the overlap function asymptotically matches 
the solution of the Schrodinger equation (20), with V = Vc for protons and V = for 
neutrons. That is, when r — > oo. 



where the binding energy of the n + x system is related to k by means of -E5 = /i^K^/2m„^, 
Wp^q is the Whittaker function and is the modified Bessel function. In Eq. (26), Q is 




for protons 



for neutrons 



(26) 
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Fig. 5. Comparison of various radial overlap integrals r^/^(r) for ^''F*(l/2+) with the normalized 
Whittaker function (dashed curve). Most of the contribution to the rms radius comes from the 
region outside the core, with radius r^- 

the asymptotic normalization coefficient (ANC). In Figure 5 we show the comparison of the 
ANC for ^''F (1/2"'') as a function of the distance r, with the Whittaker function, Eq. (26). 
As can be seen, most of the contribution to the rms radius comes from the region outside 
the core. 



In the calculation of o^l j^^^' above, one often meets the situation in which only the asymptotic 
part of Ihir) ci^nd ^ df) contributes significantly to the integral over r. In these situations, 
\&c(?") is also well described by a simple two-body scattering wave (e.g. Coulomb waves). 
Therefore, the radial integration in (f^j'^^'^' can be done accurately and the only remaining 
information from the many-body physics at short distances is contained in the asymptotic 
normalization coefficient Cj, i.e. cr^^j^^^''' oc Cf. We thus run into an effective theory for 
radiative capture cross sections, in which the constants Ci carry all information about the 
short-distance physics, where the many-body aspects are relevant. It is worthwhile to mention 
that these arguments are reasonable for proton capture at very low energies because of the 
Coulomb barrier. 



The asymptotic normalization coefficients, Ca, can also be obtained from the analysis of 
peripheral transfer and breakup reactions. As the overlap integral, Eq. (26), asymptotically 
becomes a Whittaker function, so does the single-particle bound-state wavefunction Mq,, 
calculated with Eq. (20). If we label the single-particle ANC by 6j, then the relation between 
the ANC obtained from experiment, or a microscopic model, with the single-particle ANC 
given by {SF)ibj = Cf (this becomes clear from Eq. (24)). The values of {SF)i and bi 
obtained with the simple potential model are useful telltales of the complex short-range 
many-body physics of radiative capture reactions. One can also invert this argumentation 
and obtain spectroscopic factors if the Ci are deduced from a many-body model, or from 
experiment, and the hi are calculated from a single-particle potential model [42]. 
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Microscopic calcTilations of ANCs rely on obtaining the projection, or overlap, of the many- 
body wave functions of nuclei A and A — 1. The overlap integral {{A — = /f,(r) must 
have correct asymptotic behavior with respect to the variable r which is the distance between 
the nucleon N and the cm. of the nucleus A — \. The most common methods are: (a) the 
resonating group method (RGM), as described above, (b) the Fadeev method for three-body 
systems, (c) a combination of microscopic cluster method and 7^-matrix approaches, to be 
discussed later, (d) Green's function Monte- Carlo method, (e) no-core shell model, or (f) 
hyperspherical functions method. As an example, early applications of the ANC method 
have obtained S'i7(0) = 15.5 eV b for the ^Be(p,7)^B reaction using ANCs calculated with 
^hjj oscillator wave functions and M3Y(E) effective A^A^ potential as a model for |^) and 
((A — 1) A^l [43]. The M3Y A^A^ interaction is an effective interaction constructed as in 
Eq. (19), with Vf,jj given in terms of sums of (3) Yukawa functions. 



3.2.2 Threshold behavior and the r-process 

The threshold behavior of radiative capture cross sections is fundamental in nuclear astro- 
physics because of the small projectile energies in the thermonuclear region. For example, 
for neutron capture near the threshold, the cross section can be written as [44] 

where Cq is the logarithmic derivative for the s wave at a channel radius. Since Cq is only 
weakly dependent on the projectile energy, one obtains for low energies the well-known 1/v- 
behavior. 

With increasing neutron energy, higher partial waves with / > contribute more significantly 
to the radiative capture cross section. Thus the product av becomes a slowly varying function 
of the neutron velocity and one can expand this quantity in terms of v or y/E around zero 
energy, 

av = -5^") (0) + S^""^ iO)VE + (0) ^ + . . . . 

The quantity S^^'^E) = av is the astrophysical .S-factor for neutron-induced reactions and 
the dotted quantities represent derivatives with respect to i.e., 

^H = 2v^^ and ^(")=4E^ + 2^. 

dE dE^ dE 

The astrophysical S'-factor for neutron-induced reactions is different from that for charged- 
particle induced reactions. In the astrophysical ^'-factor for charged-particle induced reac- 
tions also the penetration factor through the Coulomb barrier has to be considered (Eq. (11)). 
Inserting this into Eq. (5), one obtains for the reaction rate of neutron-induced reactions 

(av) = S{0) +[^y 3(0) (kBT) ^ + ^S{0)kBT + .... (27) 
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In most astrophysical neutron-induced reactions, neutron s-wavcs will dominate, resulting 
in a cross section showing a l/-;;-bchavior (i.e., (t{E) oc 1/^/E). In this case, the reaction 
rate will become independent of temperature, R = const. Therefore it will suffice to measure 
the cross section at one temperature in order to calculate the rates for a wider range of 
temperatures. The rate can then be computed very easily by using 

R = {av) = {a)rpVT = const. , (28) 

with vt = {2kT/mf''^. 

The mean lifetime of a nucleus against neutron capture, i.e., the mean time between 
subsequent neutron captures, is inversely proportional to the available number of neutrons 
rin and the reaction rate i?n7, Tn = (nn-Rn7)~^- If this time is shorter than the /3-decay half- 
life of the nucleus, it will be hkely to capture a neutron before decaying (r process). In this 
manner, more and more neutrons can be captured to build up nuclei along an isotopic chain 
until the /?-dccay half-life of an isotope finally becomes shorter than Tn. With the very high 
neutron densities encountered in several astrophysical scenarios, isotopes very far off stability 
can be synthesized. 



3.2.3 Halo nuclei 

For low values of the binding energy e.g. for halo-nuclei., the simple 1/f-law does 

not apply anymore. A significant deviation can be observed if the neutron energy is of the 
order of the |£^B|-value. For radiative capture to weakly-bound final states, the bound-state 
wave function Uij{r) in Eq. (23) decreases very slowly in the nuclear exterior, so that the 
contributions come predominantly from far outside the nuclear region, i.e., from the nuclear 
halo. For this asymptotic region, the scattering and bound wave functions in Eq. (23) can 
be approximated by their asymptotic expressions neglecting the nuclear potential, ui{kr) oc 
ji{kr), and Mio(r) oc h\^\iC,r), where ji and are the spherical Bessel, and the Hankel 
function of the first kind, respectively. The separation energy \Eb\ in the exit channel is 
related to the parameter ^ by \Eb\ — /i^^^/(2m„-r). 

Performing the calculations of the radial integrals in Eq. (23), one readily obtains the energy 
dependence of the radiative capture cross section for halo nuclei [45]. For example, for a 
transition s ^ p it becomes 

(rc) / N I {E + 3\EB\y , . 

'^(El)i^ ^ P) OC ^ z:. , 1^ I , (29) 



while a transition p s has the energy dependence 



'E 
E -\- \Eb\ 



a;g)(p^s)^ ^ ; ,^ , . (30) 



If ii^ -C \Eb\-i the conventional energy dependence is recovered. From the above equations 
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one obtains that the reaction rate is not constant (for s-wave capture) or proportional to T 
(for p-wave capture) in the case of small lE^sl-values. 

In the case of charged particles, S{E) is expected to be a slowly varying function in energy 
for non-resonant nuclear reactions. In this case, S (E) can be expanded in a McLaurin series, 

as was done to obtain Eq. (27). Using the expansion in Eq. (12) and approximating the 
product of the exponentials exp {—E/k^T) and exp [Inr] [E)] by a Gaussian centered at the 
energy Eq, Eq. (12) can be evaluated as [46] 



{av) = ( 



2 y/2 



rriah 



A 



^eff {Eq) exp 



3^ 
' kT 



(31) 



with 



S,s{Eq)^S{Q) 




(32) 



The quantity Eq defines the effective mean energy for thermonuclear fusion and is given by 
Eq. (13). The quantity r is given by r = SEo/kT, and A is given by Eq. (13). 



3.2.4 Resonances 

For the case of resonances, where Er is the resonance energy, we can approximate a (E) by 
a Breit-Wigner resonance formula, 

TT^l^ (2J« + 1) T,T, 

""'^^^ 2„E (2J„ + 1)(2J, + 1) (£;, - Ef + (rtot/2)^ ' ^ ^ 

where Jr, Ja, and are the spins of the resonance and the nuclei a and 6, respectively, and the 
total width Ftot is the sum of the particle decay partial width Fp and the 7-ray partial width 
F^. The particle partial width, or entrance channel width, Fp, can be expressed in terms of 
the single-particle spectroscopic factor SFi and the single-particle width Fg.p. of the resonance 
state, Fp = SFi x Fg.p.. The single-particle width Fg.p. can be calculated from the scattering 
phase shifts of a scattering potential with the potential parameters being determined by 
matching the resonance energy. The 7 partial widths F^ are calculated from the reduced 
electromagnetic transition probabilities -B( Jj — > Jf]L) which carry the nuclear structure 
information of the resonance states and the final bound states. The reduced transition rates 
are usually computed within the framework of the nuclear shell model. 

Most of the typical transitions are Ml or E2 transitions. For these, the relations are 

FE2[eV] = 8.13 X 10-^ E^^ [MeV] B{E2) [e^fm^] (34) 

and 

FMi[eV] = 1.16 X 10-2 E^^ [MeV] B{Ml) [/x^] . (35) 
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For the case of narrow resonances, with width F <^ Ej., the MaxweUian exponent exp (—E/kT) 
can be taken out of the integral, and one finds 



(av) 



where the resonance strength is defined by 



For broad resonances, Eq. (12) is usually calculated numerically. An interference term has 
to be added. The total capture cross section is then given by [47] 

aiE) = a^^iE) + a,{E) + 2 [a^,{E)a,{E)]^^'' cos[5^{E)] . (38) 

In this equation 5-^{E) is the resonance phase shift. Only the contributions with the same 
angular momentum of the incoming wave interfere in Eq. (38). 



3.3 TZ-matrix theory 



Reaction rates dominated by the contributions from a few resonant or bound states are 
often extrapolated to energies of astrophysical interest in terms of TZ- or K -matrix fits. The 
appeal of these methods rests on the fact that analytical expressions can be derived from 
underlying formal reaction theories that allow for a rather simple parametrization of the data. 
However, the relation between the parameters of the 7?.-matrix model and the experimental 
data (resonance energies and widths) is only quite indirect. The X-matrix formalism solves 
this problem, but suffers from other drawbacks [48]. 



3.4 Elastic and inelastic scattering reactions 



In the 7^-matrix formalism, the eigenstates of the nuclear Hamiltonian in the interior region 
of a nucleus are denoted by Xa, with energy £^a, and are required to satisfy the boundary 
condition 

ar 

at the channel radius r = R, where the constant 6 is a real number. The true nuclear 
wavefunction ^ for the compound system is not stationary, but since the Xx form a complete 
set, it is possible to expand ^ in terms of Xx, i.e. 

R 



^^Y^AxXx, where Ax^f Xx'lfdr. 
X Jo 
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Fig. 6. Total S factor data (filled-in circles) [52] for ^^C(a, 7)"'^^0 compared with El (open trian- 
gles) and E2 (open squares) contributions [53,54]. The solid line represents the sum of the single 
amplitudes of an 7^-matrix fit [26] (the dotted and dashed lines are the El and E2 amplitudes, 
respectively). In addition, the TZ matrix fit of [55] to their data (dotted-dashed line) is shown. 
(Adapted from Ref. [56]). 



The differential equations for \I/ and Xx are (for s-wave neutrons) 



2m dr"^ 
2m dr"^ 



+ VXx = ExXx 



r < R . 



(39) 
(40) 



Multiplying Eq. (39) by \E' and Eq. (40) by Xx, subtracting and integrating, we have 

1 



A, 



E- Ex 2mR 



Xx (R) [R^'{R) + b^{R)] 



where the prime indicates the differentiation with respect to r. This result, together with 
the definition of gives 

^{R) = n [R-^'{R) + b^{R)] (41) 
where the function TZ relates the value of "^{R) at the surface to its derivative at the surface: 



h'^^XxiR)XxiR) 



2mR 



E,-E 



(42) 



Rearranging Eq. (41) we have R'$'{R)/'${R) = (1 — bTZ)/TZ, which is just the logarithmic 
derivative which can be used to determine the S'-matrix element 5*0 in terms of the TZ 
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function. This gives 



^0 



1 + 



2ikRn 
1 - (6 + ikR)n 



-2ikR 



Finally, we assume that E is near to a particular Ex, say Ea, neglect all terms A 7^ a in Eq. 
(42), and define 



„- and A„ 

m 



2kR 



r 



SO that the ^'-matrix element becomes 



^0 



1 + 



{E^ + Aa-E)-irj2_ 



-2ikR 



(43) 



and the scattering cross section is 

^2ikR _ ]^ _|_ 



TT 



A;2 



{E^ + Aa-E)- iTj2 



(44) 



We see that the procedure of imposing the boundary conditions at the channel radius leads 
to isolated s-wave resonances of Breit-Wigner form. If the constant b is non-zero, the position 
of the maximum in the cross section is shifted. The level shift does not appear in the simple 
form of the Breit-Wigner formula because E^ + is defined as the resonance energy. In 
general, a nucleus can decay through many channels and when the formalism is extended 
to take this into account, the 7?.-function becomes a matrix. In this 7?.- matrix theory the 
constant b is real and Xx (R) and Ex can be chosen to be real so that the eigenvalue problem 
is Hermitian [49]. 



The 7^-matrix theory can be easily generalized to account for higher partial waves and spin- 
channels. If we define the reduced width by 7I = fi^Xl {R)/2mR, which is a property of a 
particular state and not dependent of the scattering energy E of the scattering system, we 
can write 



7Aa'7Aa 

Ex — E 



where a is the channel label. 7^0, Ex, and b are treated as parameters in fitting the experi- 
mental data. If we write the wavefunction for any channel as ^ ~ / -|- SaO, where / and O 
are incoming and outgoing waves, Eq. (41) means 



R 



r (R) + SaO' (R) _i- bn 

I (R) + SaO (R) ~ n ' 



Thus, as in Eq. (43), the ^'-matrix is related to the 7^- matrix and from the above relation 
we obtain that. 



Sa 



0{R) 



i-an 



(45) 
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The total cross sections for states with angular momenta and spins given hy I, s and J is 

O-aa' ^ X! gj\SaJls,a'J'l's'f , a ^ a' , (46) 

'^a JJ'll'ss' 

where gj are spin geometric factors. 

In the statistical model, it can be argued that because the jS'-matrix elements vary rapidly 
with energy, the statistical assumption implies that there is a random phase relation between 
the different components of the ^'-matrix. The process of energy averaging then eliminates 
the cross terms and gives 

(7abs = ^«a' ^ ^Y9j[^ - \SaJlsf] = ^^QjT/s («) , (47) 



a'j^a,J'l's' '"a Jls '"a Jls 



where the symmetry properties of the S'-matrix in the form Yla'yta,j' Saj,a' j' S^j a'j' — ^ 
with J = Jls, J' = J'l's' have been used, and we have introduced the general definition of 
the transmission coefficient 

Tii{a) = l-\S^ju\\ (48) 



3.5 Radiative capture reactions 



Consider an 7?.-matrix calculation of the radiative capture reaction n + a;—>a + 7 to a state 
of nucleus a with a given spin Jf. The cross section can be written as [50] ctj^ = J2ji ^JiJ^i 
with 2 J ■ + 1 

^ B (2J„ + 1)(2J, + 1) U ^'^"''f'' ' ■ ^^^^ 

Here, Jj is the total angular momentum of the colliding nuclei n and x in the initial state, J„ 
and J-c are their spins, and /, k, and li are their channel spin, wave number and orbital angular 
momentum in the initial state. Tji-j^.j. is the transition amplitude from the initial continuum 
state {Ji,I,li) to the final bound state {Jf,I). In the one-level, one-channel approximation, 
the resonant amplitude for the capture into the resonance with energy Er^ and spin Jj, and 
subsequent decay into the bound state with the spin Jf can be expressed as 

^lUJfJi - ^e » » — - —Tj—- l^UJ 

Here we assume that the boundary parameter is equal to the shift function at resonance 
energy and 0;. is the hard-sphere phase shift in the Zjth partial wave. 



0/. = arctan 



(51) 



where Fi^ and Gi. are the regular and irregular Coulomb functions, is the channel radius, 
and CT;. is the Coulomb phase factor, = Xlfe=i arctan(?7i//c), where rji is the Sommerfeld 
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parameter. r^'ji.(E) is the observable partial width of the resonance in the channel n + x, 
^'^j^ (E) is the observable radiative width for the decay of the given resonance into the bound 
state with the spin Jf, and Tj. J2 T^j;. is the observable total width of the resonance level. 



I 

The energy dependence of the partial widths is determined by 



and 




^nik(E) = -j^f^SXiERj (52) 

where r^\i.{ER^) and r'!^j^{Eii^) are the experimental partial and radiative widths, £/ is 
the binding energy of the bound state in nucleus a, and L is the multipolarity of the 7-ray 
transition. The penetrability Pii{E) is expressed as 



The non-resonant amplitude can be calculated by 



<NR 



_(2)3/2^ii+i-i/+lpj 



L+l/2 



-lUJ.j,- K^} ^ - - ' ' L[(2L + 1)!!] 

X G,, ( A;, re) W^i, (2/^rc) V^(ii OLO I i/0) C/(Li/ /; J/) Ji (ii i/) , 
where, is the effective charge, ?7 is a geometric coefficient, and 



(L + 1)(2L + 1) 



CjfiifFi.{k,rc) 



,L+1 



ar r- 



F,(A;,r) Gi^{k,r) 



FiXk,r,) GiXk,r,) 



(55) 



(56) 



Wi{2Kr) is the Whittaker hypergeometric function, k — ^/2iinx£f and If are the wave number 
and relative orbital angular momentum of the bound state, and k^ — {E + ef)/hc is the wave 
number of the emitted photon. 



The non- resonant amplitude contains the radial integral ranging only from the channel radius 
Tc to infinity since the internal contribution is contained within the resonant part. Further- 
more, the 7^-matrix boundary condition at the channel radius implies that the scattering 
of particles in the initial state is given by the hard sphere phase. Hence, the problems related 
to the interior contribution and the choice of incident channel optical parameters do not 
occur. Therefore, the direct capture cross section only depends on the ANC and the channel 
radius Tc. 

The 7^-matrix method described above can be extended to the analysis of other types of 
reactions, e.g. transfer reactions [57]. The goal of the 7^-matrix method is to parameterize 
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some experimentally known quantities, such as cross sections or phase shifts, with a small 
number of parameters, which are then used to extrapolate the cross section down to astro- 
physical energies. One example is given in Figure 6 which shows the experimental data and 
7?.-matrix fits for the cross section of the reaction ^^C(q;, 7)^^0 cross section, of relevance to 
helium burning [56]. 



3.6 Statistical models 



A large fraction of the reactions of interest proceed through compound systems that exhibit 
high enough level densities for statistical methods to provide a reliable description of the 
reaction mechanism. The theoretical treatment of nuclear reactions leading to formation and 
decay of compound nuclei was developed by Ewing and Weisskopf [58], based on two ideas: 
(a) the compound nucleus formation independence hypothesis as proposed by Niels Bohr [59], 
and (b) the reciprocity theorem, or time-reversal properties of the underlying Hamiltonian. 
This allows one to relate capture and decay cross sections, usually expressed in terms of 
transmission probabilities, defined in Eq. (48). 

Later, the Ewing- Weisskopf theory was extended to include angular momentum dependence 
by Hauser and Feshbach [60] . The Hauser-Feshbach (HF) model has been widely used with 
considerable success in nuclear astrophysics. Explosive burning in supernovae involves in 
general intermediate mass and heavy nuclei. Due to a large nucleon number, they have 
intrinsically a high density of excited states. A high level density in the compound nucleus 
at the appropriate excitation energy allows for the use of the statistical-model approach for 
compound nuclear reactions [60] which averages over resonances. 

A high level density in the compound nucleus also allows for the use of averaged transmission 
coefficients T, which do not reflect resonance behavior, but rather describe absorption via 
an imaginary part of the (optical) nucleon-nucleus potential as described in Ref. [51]. This 
leads to the expression derived 

(^'«'^^.-) = (2jr + i)(2j, + i) 

X ^(2 J + 1) Zn^iJ^^ <)To'{E: J, vr, E;;, J^, tt^;) 



J,7r 



Ttot{E, J, tt) 



for the reaction 0)171" from the target state to the excited state of the flnal 
nucleus, with a center of mass energy E^j and reduced mass /lij. J denotes the spin, E the 
corresponding excitation energy in the compound nucleus, and tt the parity of excited states. 
When these properties are used without subscripts they describe the compound nucleus, 
subscripts refer to states of the participating nuclei in the reaction i'^(j, 0)171" and superscripts 
indicate the specific excited states. Experiments measure o"°^(j, o; i?jj), summed over all 
excited states of the final nucleus, with the target in the ground state. Target states fi in an 
astrophysical plasma are thermally populated and the astrophysical cross section a* (j, o) is 
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Fig. 7. Cross-section data for two sets of '^^As{p,n) measurements (squares: [62]; circles: [63]) in 
comparison with Hauser-Feshbach predictions (solid line: [64]). Also shown is the experimental cross 
section of ^^Rb{p, n) (crosses with error bars: [65]) in comparison with HF predictions (dashed line: 
[64]). (Adapted from Ref. [61]). 



given by 



EMJt + 1) expi-Ei'/kT) Eu crrU, o; E, 



(5^ 



E,(2Jr + l)exp(-Er/fcT) 
The summation over v replaces T^{E, J, vr) in Eq. (57) by the total transmission coefficient 



T„(E,j,7r) = ^T;(E,j,7r,E;;,j;;, 



E — Sm,o 



+ / X! To{E, J,TT, Em, Jm,T^m)p{Em, Jm,TTm)dEm- 



(59) 



Here Sm,o is the channel separation energy, and the summation over excited states above the 
highest experimentally known state z/^ is changed to an integration over the level density p. 
The summation over target states /i in Eq. (58) has to be generalized accordingly. 



The important ingredients of statistical-model calculations, as indicated in the above equa- 
tions, are the particle and 7-transmission coefficients T and the level density of excited states 
p. Therefore, the reliability of such calculations is determined by the accuracy with which 
these components can be evaluated (often for unstable nuclei). 

Figure 7 - adapted from Ref. [61] - shows the cross-section data for two sets of ''^As(p, n) 
measurements (squares: [62]; circles: [63]) in comparison with Hauser-Feshbach predictions 
(solid line: [64]). Also shown is the experimental cross section of *^Rb(p, n) (crosses with 
error bars: [65]) in comparison with HF predictions (dashed line: [64]). The experimental 
results are on average lower than the HF predictions. 
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3.7 Spin-isospin response 



Beta-decay, electron capture and neutrino scattering involve similar operators and nuclear 
matrix elelments. Wc thus consider only the case of neutrino scattering. In the following, = 
{p£, El} and Qi, = {q, Ei,} are the lepton and the neutrino momenta, k = Pi — Pj = {k, ^0}, 
is the momentum transfer. Pi and Pf are momenta of the initial and final nucleus, M is the 
nucleon mass, m^ is the mass of the charged lepton, and Qv, Qa, Qm and Qp are, respectively, 
the vector, axial-vector, weak-magnetism and pseudoscalar effective dimensionless coupling 
constants. Their numerical values are typically given hy Qv — Qa — ^■'^^i Qm — >^p — >^n — 
3.70, and = ^^(2Mm^)/(A;2 + m^). 

For the neutrino-nucleus reaction the momentum transfer is k = pi — q^,, and the correspond- 
ing cross section reads 

a{Ei, Jf) = ^^^F{Z ± 1, E,) £ d(cos^)r,(q, J^), 

where 

is the Fermi ftmction {Z± = Z + 1, for neutrino, and Z ~ 1, for antineutrino) , ^ = q • p is 
the angle between the incident neutrino and ejected lepton, and the transition amplitude for 
initial (final) angular momentum Jj (J/) is 

Z -h i g^^g^ Mi,Mf 

One can cast the transition amplitude in the compact form [66] ^ 



(J/||O0j|| Ji) |^£0 + ^ I (^/||OMj||^j) |^>Cm 

M=0±1 



.23fJ(|(J/||O0j||J,)(J/||Ooj||J.))/:0 



(60) 



where G = (3.04545 ± 0.00006) x 10 is the Fermi coupling constant (in natural units), and 



|p|cosg f q.p. 



£„.A = l + |f -l!^, l. = i-|^±(|-|)s., (61) 



^ The indices and z denote the time-component and the third-component of four-vectors, respec- 
tively. 
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with 

g.^fe-q^ ^-^'P'^"^^-^-\ p, = Lp=M^:^^^, (62) 

K, K 

being the z-components of the neutrino and lepton momenta, and = ±1 for neutrino 
scattering and antineutrino scattering, respectively. 

Exphcitly, the operators in (60) arc 

O0J = QvM] + 2ig^M^ + i{g^ + gp,)Mtj, 

Omj = i{5Mz9p2 - 9a + Mg^)MM^ + 2^v^mj - ^MzVyM] , (63) 
where k = k/«;, k,= |k|, and the following short notation 

^^"^^2M' ^^"^^2M' ^^"^^^ + ^^W'^""^^2mS' ^^'^^"TMme' ^^^^ 
has also been introduced. The elementary operators are given by 

M]=J^{p)Ury, = K''Up)Ur){a ■ V); 

■Mmj = E^'"'"' F^uJiip) [YL{i) ® ; Al^j = ac"^ E^'"'""'^mljJl(p) [>!(?) ® V]j, 

L L 

(65) 

where p — nr. Notice that the initial and final states of the matrix elements in Eq. (60) 
involve an isospin unit change, and implicitly contain isospin operators t±. 



3.7.1 Fermi and Gamow-Teller matrix elements 

Most reactions in typical stellar scenarios involve small momentum transfer such that p <^ 
1. In this case, the angular dependence of the above operators becomes irrelevant. Using 
jmip) ~ ^mo for the spherical Bessel functions in Eqs. (65) and after some algebra one can 
show that (for charged-current) 



Jf) ~ C [I {Ji\\T,Ur±{k)U) r + 9a\ {Ji\\T.U<k)T^m\Ji) V\ , (66) 

where C is a function depending on the lepton and neutrino energies. The r+ operator 
corresponds to /3~ decay and the t_ to decay, so that t+ | n) = | p) and t_| p) = | n), 
changing a neutron into a proton and vice- versa. 

The spin-independent and spin-dependent operators appearing on the right-hand-side of the 
above equation are known as the Fermi and Gamow-Teller operators. The Fermi operator is 
the isospin raising/lowering operator: in the limit of good isospin, which typically is good to 
5% or better in the description of low-lying nuclear states, it can only connect states within 
the same isospin multiplet. That is, it is capable of exciting only one state, the state identical 
to the initial state in terms of space and spin, but with (T, Mr) = {Ti, Mn ± 1) for /3~ and 
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Fig. 8. Number of events obtained from the convolution of the neutrino fluxes with the cross 
section obtained with different nuclear structure models: Hybrid (dashed-dot line) [70], quantum 
phase approximation (QRPA) (dashed line), QRPA^ (dashed-dot dot line) [71], projected QRPA 
(PQRPA) (sohd line) [67], and renormalized QRPA (RQRPA) (dashed line) [72]. (Adapted from 
Ref. [67]). 



/3+ decay, respectively. 



Eq. (66) is only appropriate at the lowest-order expansion in p, when the nucleus responds like 
an elementary particle. Then we can characterize its response by its macroscopic quantum 
numbers, the spin and charge. The next-to-leading-order term of the expansion in powers 
of p probes the nucleus at shorter length scales. The operators in Eq. (65) are obtained by 
an expansion of the plane-wave lepton wavefunction, which, for not too large k, becomes 
exp(zk-r) ~ l-|-zk-r. Thus, the next term in the expansion includes a "first forbidden" term 
X^ili riT3(i) and similarly for the spin operator Yl!i=^i ® ''■(0]j=o,i, 2^3(0- These operators 
generate collective radial excitations, leading to the so-called "giant resonance" excitations, 
with a typical excitation energy of 10-25 MeV. They tend to exhaust the Thomas-Reiche- 
Kuhn sum rule, 

El(/lEK0^3(^)K)P~-^~y (67) 

where the sum extends over a complete set of final nuclear states. The first-forbidden opera- 
tors tend to dominate the cross sections for scattering the high energy supernova neutrinos 
{Ufj^s and z/,-s), with Ei, ~ 25 MeV, off light nuclei. It also follows from Eq. (67) that the cross 
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sections per target nucleon are roughly constant. This conclusion changes when high energy 
neutrinos, with high energy transfers, are considered. 

The number of events detected for supernova explosions can be calculated as. 



where the index a = Ug, z/g, and [u^ = i^r, ^^i, ^r) indicates the neutrino or antineutrino 
type, Nf is the number of target nuclei, Ta^E^,) is the neutrino flux, a{Ei,) is the neutrino- 
nucleus cross section, e{Ei,) is the detection efficiency, and E^, is the neutrino energy. 

In Figure 8 (from reference [67]), we show calculations for the number of events of detected 
supernova neutrinos due to (z/g + z/g) interactions on ^^Fe, such as those performed in the 
KARMEN collaboration [68]. The A^e and Ne are calculated as a function of the neutrino 
temperatures T^^ and T,^^, folding (Te{E^) from different nuclear structure models with the 
neutrino fluxes Ti,^{E^,T,^^) and T^^{E^,T^^), respectively [69]. The fluxes depend on the 
distance to the supernova, the neutrino energy, and the neutrino effective temperature. The 
number of events are obtained from the convolution of the neutrino fluxes with the cross sec- 
tion obtained with different nuclear-structure models: Hybrid (dashed-dot line) [70] , quantum 
phase approximation (QRPA) (dashed line), QRPA5 (dashed-dot dot hne) [71], projected 
QRPA (PQRPA) (sohd line) [67], and renormahzed QRPA (RQRPA) (dashed hne) [72]. 
One clearly sees that the differences between the calculated cross sections with different nu- 
clear models increase as a function of the neutrino temperatures. This is an example of the 
limitations of nuclear models in describing weak-interaction processes in stars. 



3.8 Field theories 

Field theories adopt a completely independent approach for nuclear physics calculations in 
which the concept of nuclear potentials is not used. The basic method of field theories is to 
start with a Lagrangian for the fields. Prom this Lagrangian one can "read" the Feynman 
diagrams and make practical calculations, not without bypassing well-known complications 
such as regularization and renormalization. Quantum chromodynamics (QCD) is the proper 
quantum field theory for nuclear physics. But it is a very hard task to bridge the physics 
from QCD to the one in low-energy nuclear processes. Effective field theory (EFT) tries to 
help in this construction by making use of the concept of the separation of scales. One can 
form small expansion parameters from the ratios of short and long distance scales, defined 



long distance scales 
and try to explain physical observables in terms of powers of e. 

In low-energy nuclear processes, the characteristic momenta are much smaller than the mass 
of the pion, which is the lightest hadron that mediates the strong interaction. In this regime. 




by 



short distance scales 



(68) 
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one often uses the pionless effective field theory, in which pions are treated as heavy particles 
and are integrated out of the theory [76] . In this theory, the dynamical degrees of freedom are 
nucleons and the pion and the delta resonance degrees of freedom are hidden in the contact 
interactions between nucleons. The scales of the problem are the nucleon-nucleon scattering 
length, a, the binding energy, B, and the typical nucleon momentum k in the center-of-mass 
frame. Then, the nucleon-nucleon interactions are calculated perturbatively with the small 
expansion parameter 

{l/a,B,k) 

P = (69) 

which is the ratio of the light to the heavy scale. The heavy scale A is set by the pion mass 
(m^ ~ 140 MeV). 

The pionless effective Lagrangian will only involve the nucleon field = {p, n) and its 
derivatives. It must obey the symmetries observed in strong interactions at low energies, 
such as parity, time-reversal, and Galilean invariance. The Lagrangian can then be written 
as a series of local operators with increasing dimensions. In the limit where the energy goes 
to zero, the interactions of lowest dimension dominate. To leading order (LO), the relevant 
Lagrangian [h = c = 1) is given by 

£ = I^Q^ + Z! ) ^ _ Co(^^P^)(^^P^)^ (70) 
Y 2m J 

where m is the nucleon mass [73]. The projection operators V enforce the correct spin and 
isospin quantum numbers in the channels under investigation. For spin-singlet interactions 
'P'i = o'2'r2Ti/ VS, while for spin-triplet interactions Vi = 0-2(^1^2/ V8- 

The Feynman-diagram rules can be directly "read" from the Lagrangian at hand. In the case 
that the scattering length a is large, i.e., a 3> 1/A, as it is in the nucleon-nucleon system, 
the full scattering amplitude T is obtained from an infinite sum of such Feynman diagrams 
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(see Figure 9) , leading to a geometric series that can be written analytically as 

^^^^ = I-cIapV -^^"^ ^ I WTE-cVm + ie^ ^^'^ 

where E = p^/m is the total center-of-mass energy (sec, e.g., Ref. [75]). The integral is 
linearly divergent but is finite using dimensional regularization. One gets 

Tfi 

where fi is the regularization parameter. The scattering amplitude T[p) has then the same 
structure as the s-wave partial- wave amplitude, 



p cot 5 — ip^ 

and one obtains the effective range expansion for the phase shift S, 

pcotd = h To— H 

a 2 

in the zero-momentum limit when the coupling constant takes the renormalized value 

Co(/i) = (72) 
m [1/a — /I) 



In leading order, we see that the effective range vanishes, or tq = 0. The small inverse 1/a 
scattering length is given by the difference between two large quantities. For example, in 
proton-neutron scattering we have ap„ = —23.7 fm in the pn spin-singlct channel. Choosing 
the value ^ = m^r for the regularization parameter, one obtains Co = 3.54 fm^. Physical 
results should be independent of the exact value of the renormalization mass /i as long as 
1/a < n <^ rriT^. 



The theory is now ready for practical applications. For example, this procedure has been 
applied to obtain the electromagnetic form factor of the deuteron, the electromagnetic po- 
larizability and the Compton scattering cross section for the deuteron, the radiative neutron 
capture on protons, and the continuum structure of halo nuclei. Based on the same effective 
field theory, the three-nucleon system and ncutron-deuteron scattering have been investi- 
gated [75]. Better agreement with data can be obtained in higher orders (next-to-leading 
order [NLO], next-to-next-to leading order [N'^LO], etc.). For nuclear processes involving 
momenta p comparable to mTr, the starting effective, pionfuU, Lagrangian is more compli- 
cated. But the basic field theoretic method remains the same. The EFT unifies single-particle 
approaches in a model-independent framework, with the added power counting that allows 
for an a priori estimate of errors. Concepts of quantum field theory, such as regularization 
and renormalization, are key ingredients of the theory. 
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Fig. 10. The cross section for d np. The curves correspond to EFT cacluations for cold np ^ d 
and the dashed lines denote the a 3% theoretical uncertainty. Ref. [74] has further reduced this 
uncertainty to below 1%. (Adapted from Ref. [74]). 

In nuclear astrophysics, this theory has been applied to np d'y for big-bang nucleosynthesis 
[77,74]; ud reactions for supernovae physics [79] and the solar pp fusion process [80]. EFT has 
also been used to deduce observables in reactions with halo nuclei and loosely bound states, 
with promising applications to astrophysics [81-83]. So far, perhaps the most enlightening 
application of EFT for nuclear physics is the np d'y cross section, specially because there 
is no data at the energies of relevance for the big bang nucleosynthesis (BBN). EFT has 
provided a calculation with 1% error [74] in the energy range relevant to BBN. The EFT 
predictions also agree with a very recent measurement of the inverse process in the same 
energy region (see Figure 10). 



4 Effects of environment electrons 

The form of the astrophysical S factor given in Eq. (11) assumes that the electric charges 
of nuclei are "bare" charges. However, neither at very low laboratory energies, nor in stellar 
environments this is the case. In stars, the bare Coulomb interaction between the nuclei is 
screened by the electrons in the plasma surrounding them. If one measures reaction rates in 
the laboratory, using atomic targets (always), then atomic electrons screen as well. But the 
screening is different from the screening in the stellar plasma. Therefore we discuss these two 
problems separately in the following subsections. 

4-1 Stellar electron screening problem 

In astrophysical plasmas with high densities and/or low temperatures, effects of electron 
screening are very important, as will be discussed later. This means that the reacting nuclei, 
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due to the background of electrons and nuclei, feel a different Coulomb repulsion than in the 
case of bare nuclei. Under most conditions (with non- vanishing temperatures), the general- 
ized reaction-rate integral can be separated into the traditional expression without screening 
(4) and a screening factor 

= fscriZj, Zk, p, T, Yi) (crv)^. , (73) 
in terms of the nuclear abundances, defined in section 2.2. 



This screening factor is dependent on the charge of the involved particles, the density, temper- 
ature, and the composition of the plasma. At high densities and low temperatures, screening 
factors can enhance reactions by many orders of magnitude and lead to pycnonuclear ignition. 

Consider a concentration of negative and positive charges with neutral total charge, that is, 
J2i ZiCCio = 0, where Qo is the spatially uniform concentration of positive {i = +) or negative 
{i = — ) charges. Because of the interaction between the charges, these concentrations are no 
more spatially uniform, with smaller charges tending to concentrate around larger charges. 

The concentrations around the charges are populated according to the statistical distribution 
of the individual charge energies in the presence of a Coulomb field V{r), yet to be found. 
Assuming Boltzmann statistics, this argumentation implies that 



c+{r) = c+oexp 



Z+eV{r) 
kT 



and c_ (r) = c_o exp 



Z_eV{r) 
kT 



(74) 



If the ion close to which we are considering the screening is positive, then V{r) > and 
c+(r) < c+o, or c_(r) > c_o, and the reverse is true if V{r) < 0. 



The charge density at position r is given by 

P(^) = Yl ^i^^i ^ ZiCCio exp 



ZieVjr) 
kT 



(75) 



If ZieV{r)/kT -C 1 {weak screening), then p(r) — —{e^V{r)/kT)Y,iZ'^CiQ. 



To obtain the potential V{r) one has to solve the Poisson equation for the potential V{r) 
which, for the above charge distribution, becomes 



1 d 

^2 



dr 



where the Debye radius Rd is defined by 



Rd — 



47rp(r) 



kT 



R 



D 



V, 



(76) 



ATre-'EiZfcio' 

Since V{r) — > as r — >• cxd, the solution of this equation is V{r) = {A/r)exp{—r/RD)- The 
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Fig. 11. Scematic view of the Debye-Hiickel sphere. An ion at the center of the sphere is surrounded 
by a cloud of ions, with the ions of opposite charge (electrons) agglomerating closer to it. 

normalization constant is fixed by the condition V{r) Zic/r as r — » 0. Thus, 

V(r)^fexp(-^) (77) 

Screening modifies the Coulomb potential between the nuclear radius R and the classical 
turning point Rq^ and consequently modifies the barrier penetration. For weak screening 
Rd ^ R, Ro- other words, we can expand V{r) around r = 0. To first order, the barrier 
energy for an incoming projectile with charge is V{r) = ZiZ2e^ /r + f/(r), where the 
Debye-Hueckel screening potential, U{r) = U{0) = const., is given by Uq = —ZiZ2e'^/RD. 

The impact of the screening potential on the barrier penetrability and therefore on the as- 
trophysical reaction rates can be approximated through a screening factor / = exp {Uo/kT), 
which, in the weak screening limit, becomes / ~ 1 — Uo/kT. 

In summary, for the weak screening limit, the reaction rates are modified according to 

{a V) screened = f{crv)bare (78) 

where 

;/ y 1/2 f 1/2 

/ = 1 + 0.188 ' ^^3^,^ , where ^ = Y.{Zf + Z,yY,. (79) 

The Debye-Hiickel model (Figure 11) is an ideal plasma since the average interaction energy 
between particles is smaller than the average kinetic energy of a particle [87,88]. In the case 
of strong screening, in low density plasmas, the potential energy cannot not be described by 
the Debye-Hiickel model since the probability of finding other charged particles in a Debye 
sphere almost vanishes. For a strongly coupled plasma, the ion-sphere model [89,90] is more 
suitable. The ion-sphere model is equivalent to the Wigner-Seitz sphere used in condensed- 
matter theory. It assumes an ion having Z^ bound electrons, positioned at r = 0, and Zf 
free electrons (Z^ + Zf = Z) occupying the rest of the ion sphere volume. The plasma effects 
are taken into account by confining the ion and the Zf electrons inside the ion sphere. To 
obtain the potential V{r) one adds to the bare ion potential, V = Ze'^/r, the potential due 
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to bound electrons, Vb, and that due to free electrons, Vf. A Slater type, or Kohn-Sham type, 
exchange-potential is also added. To obtain the bound and free electron densities one solves 
the Schrodinger equation, or Dirac equation, with V{r). From this, one builds the bound 
and free electron densities which are then used to calculate the new potentials 14 and Vf. 
This process is done iteratively until convergence is reached. 



The plasma density enters the ion-sphere model through the boundary conditions imposed 
on the potential, that is, through the neutrality conditions of the ion-sphere. Approximate 
schemes to obtain the ion-sphere potential have been developed. A widely used approximation 
for the potential energy of a single free electron electron inside the ion-spehre is given by 



V{v,r') 



+ 



r — r 



1 - 




0{Ri-r), 



(80) 



where r and r' are the positions of the bound electron and the projectile ion, respectively, 
9{Ri — r) is the step function, the ion-sphere radius Ri{= [3{Z — l)/47rne]^'^^) is given by 
the plasma electron density Ug since the total charge within the ion-sphere is neutral. This 
hydrogenic ion-sphere potential (80) can be generalized to Zf free electrons inside the ion- 
sphere. The ion-sphere model also has its limitations. For instance, charge transfer processes 
in collisions between positive ions in strongly coupled plasmas has not been fully explored 
and could modify the range of vahdity of the model [91]. 



For screening in plasmas with intermediate densities, i.e., when UeR ~ 1, where R = Ri or 
R = Rd, more complicated models are necessary and are still under theoretical scrutiny. 
This is based on the simple observation that in the stars along the main sequence, there 
is only about 1-3 ions within the Debye sphere. Thus, in principle, the Debye screening 
model should not be applicable to screening in these environments. Also, static models as 
the Debye-Hiickel and ion-sphere models do not contain dynamical effects due to the fast 
motion of free electrons. Dynamical fluctuations, due to the fast motion of the electrons, 
and non-spherical effects could modify the screening in non-static models. The possibility of 
dynamic effects was first mentioned in Ref. [92] and studied in Ref. [93]. But, the existence of 
dynamic effects was criticized since the reacting particles are in thermodynamic equilibrium 
and hence such an effect is not expected [94]. According to Ref. [95], higher order effects, 
beyond the Debye-Hiickel approximation, modify the screening enhancement in solar fusion 
reactions by only a very small amount of about 1%. This conclusion is not in accord with 
results obtained by other authors [96,91]. 



More recently, there have been additional claims that mean field models cease to be valid 
under the conditions prevailing in stellar cores in general and in the Sun, because particle 
fluctuations within the Debye-Hiickel sphere are percent-wise large. These claims have been 
substantiated with molecular dynamics calculations [96]. However, it has not been pursed 
further and has not been verified independently. This certainly deserves further theoretical 
studies. 
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Fig. 12. S'-factor data for the ^He((i,p)^He reaction from Ref. [105]. The dashed curve represents 
the /S-factor for bare nuclei and the solid curve that for screened nuclei with Up = 219 eV. 



4-2 Laboratory atomic screening problem 



Laboratory screening has been studied in more detail experimentally, as one can control dif- 
ferent charge states of the projectile+target system in the laboratory [97-99]. Experimental 
techniques improve steadily and one can measure fusion cross sections at increasingly lower 
energies where the screened Coulomb potential can be significantly smaller than the bare 
Coulomb potential. The deviation from the bare Coulomb potential is seen as an increase in 
the astrophysical S'-factor extracted at the lowest energies (see Figure 12). This enhancement 
has been experimentally observed for a large number of systems [100-104]. The screening 
effects of the atomic electrons can be calculated [97] in the adiabatic approximation at the 
lowest energies and in the sudden approximation at higher energies with a smooth transition 
in between [106]. 

In the adiabatic approximation one assumes that the velocities of the electrons in the target 
are much larger than the relative motion between the projectile and the target nucleus. 
In this case, the electronic cloud adjusts to the ground-state of a "molecule" consisting 
of two nuclei separated by a time-dependent distance -R(t), at each time instant t. Since 
the closest-approach distance between the nuclei is much smaller than typical atomic cloud 
sizes, the binding energy of the electrons will be given by the ground-state energy of the 
Zp + Zt atom, B{Zp + Zt)- Energy conservation implies that the relative energy between 
the nuclei increases by Ue = B{Zp + Zt)—B{Zt). This energy increment enhances the fusion 
probability because the tunneling probability through the Coulomb barrier between the 
nuclei increases accordingly. In other words, the fusion cross section measured at laboratory 
energy E represents in fact a fusion cross section at energy E + Ue, with Ue being known as 
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the screening potential. Using eq. (11), one gets for non- resonant reactions 



a {E + Ue) ^ exp 



nri{E) ^ 



a{E) , (81) 



where we assumed that the factor S{E)/E varies much slower with E, as compared to the 
energy dependence of exp [—2Trrj{E)] . 

The exponential factor on the right-hand-side of cq. (81) is the enhancement factor due 
to screening by the atomic electrons in the target. For light systems, the velocity of the 
atomic electrons is comparable to the relative motion between the nuclei. Thus, a dynamical 
calculation is more appropriate to study the effect of atomic screening [106]. However, the 
screening potential obtained from a dynamical calculation cannot exceed that obtained in 
the adiabatic approximation because the dynamical calculation includes atomic excitations, 
which reduce the energy transferred from the electronic binding to the relative motion. 
The adiabatic approximation is thus the upper limit of the enhancement due to laboratory 
screening. 

The experimental value of Ue needed to reproduce the experimental data by using Eq. (81) 
are systematically larger than the adiabatic model by a factor of 2 [97-99]. For example, the 
cross section of the ^He((i,p)^He reaction was studied over a wide range of energies [105]: 
the results led to Ug — 219 ± 15 eV, significantly larger than the adiabatic hmit from atomic 
physics, Uad = 119 cV. Many theoretical attempts to explain this puzzle have been carried 
out (see, e.g., refs. [106-110]). The fusion cross sections change exponentially with a small 
variation of the relative energy between the nuclei. Many small effects have been considered 
theoretically and, as shown in Ref. [107], they are not able to explain the differences between 
the experimental and theoretical values of Ue- The calculated fractional change in the cross 
sections involving light nuclei at astrophysical energies are: (a) vacuum polarization (10~^), 
(b) relativity (10~^), (c) Bremsstrahlung (10^^), (d) atomic polarization (10~^) and nuclear 
polarization (10^"'^'^) [107]. In Ref. [99] effects due to thermal motion, vibrations inside atomic, 
molecular or crystal system, and due to finite beam energy width were considered. All these 
effects are marginal at the energies, which are presently measurable (at the level of 10~^, or 
below) . 

A possible solution of the laboratory screening problem was proposed in Refs. [111,112]. Ex- 
perimentalists often use the extrapolation of the Andersen- Ziegler tables [113] to obtain the 
average value of the projectile energy due to stopping in the target material. The stopping 
is due to ionization, electron-exchange, and other atomic mechanisms. However, the extrap- 
olation is challenged by theoretical calculations, which predict a lower stopping. Smaller 
stopping was indeed verified experimentally [99]. At very low energies, it is thought that 
the stopping mechanism is mainly due to electron exchange between projectile and target. 
This has been studied in Ref. [114] in the simplest situation; proton-|-hydrogen coUisions. 
The calculated stopping power was added to the nuclear stopping power mechanism, i.e. to 
the energy loss by the Coulomb repulsion between the nuclei (Rutherford scattering). The 
obtained stopping power is shown to be smaller than the extrapolations from the Andersen- 
Ziegler, as verified experimentally by [115]. 
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The stopping power in atomic He"'"4-He collisions using the two-center molecular orbital basis 
was calculated in Ref. [116], and a good agreement with the data of Ref. [115] at low energies 
was obtained. In particular, it was found that a threshold effect exists, sharply decreasing 
the stopping power at lower energies due to the disappearance of resonant tunneling in the 
electron-exchange mechanism. The agreement with the data disappears when the nuclear 
recoil is included. In fact, an unexpected "quenching" of the nuclear recoil was observed 
experimentally in Ref. [117], for stopping of deuteron projectiles on deuteron targets. But 
this cannot be explained in terms of a threshold cutoff effect and it seems to violate a basic 
principle of nature, as the nuclear recoil is due to Coulomb repulsion (Rutherford scattering) 
between projectile and target atoms. Energy loss due to Rutherford straggling is though 
to be well described theoretically [113]. The fusion reaction d{d,p)t was recently studied 
in deuterated metals and insulators, i.e. for 58 samples across the periodic table, where a 
dramatic increase was observed for all the metals [118]. The experimentally determined values 
of the screening energy are about one order of magnitude larger than the value achieved in 
a gas-target experiment and significantly larger than the theoretical predictions [119]. This 
result has not been proven independently. 

The present status of the laboratory screening of nuclear fusion reactions is rather confusing, 
and more research is obviously needed. However, for the solar fusion reaction chains, the 
effects of laboratory and stellar screening by electrons seem to be under control. 



5 Solutions with indirect methods 



5.1 Elastic scattering 

Elastic scattering of nuclei is sensitive to their matter distribution. This is due to the depen- 
dence of the optical potential on the matter distribution of nuclei. Folding optical potentials, 
often used in the analysis of experiments, depend on the nuclear projectile, P, and target, 
T, densities as 



where Vgff is the effective nucleon-nucleon interaction and a is a parameter to fit the imagi- 
nary normalization of the optical potential. Nuclear densities are a basic input in theoretical 
calculations of astrophysical reactions at low energies. 

Elastic scattering in high-energy collisions essentially measures the Fourier transform of the 
matter distribution. Considering for simplicity the one-dimensional case, for light nuclei one 
has 



J J + a 

where q = 2ksin6/2, for a cm. momentum k, and a scattering angle 6. For heavy nuclei. 



Uopt = (l + ia) / pp{ri)pT{r2)veff{ri,r2)d^ridh2, 



(82) 




(83) 
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Fig. 13. Elastic scattering data for ("'^^C,^^C) is shown [85,86], together with data for the in- 
elastic scattering for the excitation of Gamow- Teller states in the ^^C("'^^C,"'^^N)-'^^B reaction at 
E = 135 MeV/nucleon. 

the density p is better described by a Fermi function, and 



/ 



iqx 

^^r^^(;^-(4-)s-(gi?)e--, (84) 



for i? >> a, and qa » 1. A similar result emerges from the elastic scattering amplitude, 
F{q), for the momentum transfer q, calculated in the eikonal approximation (see below), i.e., 

/ei(g) ~ / dbbJo{qb)[l - e^^(^)] ~ / dbb L^,. ~ ^JiiqR) exp(-vrga), 

J J l + exp[('^)] q 

where J„ are Bessel functions of order n. Thus, the distance between minima in elastic 
scattering cross sections measures the nuclear size, while its exponential decay dependence 
reflects the surface diffuseness. 

During recent years, elastic proton scattering has been one of the major sources of information 
on the matter distribution of unstable nuclei at radioactive beam facilities. Information on the 
matter distribution of many nuclei important for the nucleosynthesis in inhomogeneous Big 
Bang and in r-processes scenarios could also be obtained from elastic scattering experiments. 
Due to the loosely-bound character and small excitation energies of many of these nuclei, 
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high energy resolution is often necessary. But the data deviate from this simple behavior 
as soon as the energy transfer, or Q-value, differs from zero, as manifest in the inelastic 
scattering data. 



For low-energy nucleus-nucleus scattering, the distorted wave Born approximation (DWBA) 
amplitude is given by 

/dwba W = / 4^{k',r)U^,{r)xi^\k,r)d\ (85) 

where fi is the reduced mass, x^/s ^ (xi^"*) ^'^'s incoming (outgoing) distorted waves in channel 
a {(3), and is the optical potential. 



Reactions with secondary beams have been studied at relatively high energies, Ei^^^ > 
30 MeV/nucleon. The distorted waves can be approximated by eikonal waves. This is vahd 
for small-angle scattering and the scattering amplitude in the eikonal approximation is 



U{e) = ik / hdh Mqb) h - e"^^''> 



(86) 



where k — \/2/lE /his the relative momentum, and 



<j>{b)^<f>c{b) + Mb) , Mb)^ 



hv 



dzU 



(87) 



is the nuclear eikonal phase and Xc{b) is the Coulomb eikonal phase 



^c{b) ^^-^In (kb). 



Optical potentials are usually parameterized in the form 

UN{r) = -VRfv{r)-iWRfw{r)+^iaiVi^fw{r) + 2 f^") l^[Vs fs{r)] (l-s) + Koui. 



where ^ 

^'^''^ ^ l + exp[(r-i?,)M] ^^^^ 
for i = V,W and S ; with Ri = rj A^^^ . The first (second) term is the usual real (imaginary) 
part of the optical potential. The third term is peaked at the surface of the nucleus and is 
used to simulate a stronger absorption of the incoming nucleon at the surface of the nucleus. 
It is a correction due to the Pauli blocking effect. The last term in Eq. (88) is a spin-orbit 
correction. It causes interference between the scattering from opposite sides of the nucleus. 

Sometimes one has to go beyond the optical model description of inelastic scattering and 
introduce the effects of polarization potentials. Under the action of a small interaction, a 
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wavefunction is modified in lowest order to 




(90) 



If we assume tfiat {ipn) is tlie ground state, tliis equation says tliat during tfie action of the 
potential U, the wavefunction acquires small components from excited states. At the end of 
this process, the wavefunction can return to its initial state again. The modification of the 
wavefunction during the action of the potential is called "polarization" . 

While elastic scattering data is considered a simpler way to access the sizes, density profile, 
and other geometric features of nuclei, inelastic scattering requires many more pieces of infor- 
mation about the intrinsic properties of nuclei. This is shown in Figure 13, where the elastic 
scattering data for (^^C,^^C) is shown [85,86] together with data for the inelastic scattering for 
the excitation of spin-dipole states in the ^^C(^^C,^^N)^^B reaction at = 135 MeV/nucleon. 
The beautiful exponential decrease of the cross section with the nuclear diffuseness is clearly 
seen for the elastic scattering data. The inelastic data is much more sensitive to the models 
used to describe the nuclear excitation. Sometimes, further complications, such as polariza- 
tion effects, must be taken into account. The coupling to other inelastic channels also has to 
be considered. This is often necessary and complicates the nice feature of elastic scattering 
as a probe of the nuclear geometry and density profiles. 



5.2 Coulomb excitation and dissociation 

In low-energy collisions, the theory of Coulomb excitation is very well understood [120] and 
has been used to analyze experiments on multiple excitations and reorientation effects [121]. 

At the other end of the beam energy scale - Coulomb excitation of intermediate-energy or 
relativistic heavy ions - the kinematics is characterized by straight-line trajectories [122]. In 
the experiment, the selection of impact parameters b exceeding the sum of the radii of the 
two colliding nuclei by several fm (via restrictions on scattering angles) keeps the colliding 
systems at "safe" distances, minimizing the nuclear contribution to the excitation also in 
reactions above the Coulomb barrier. Most Coulomb-excitation experiments at rare-isotope 
beam facilities to date have been performed at intermediate bombarding energies of 50- 
100 MeV/nucleon. It has been a very successful tool to extract precious information on 
electromagnetic properties of nuclear transitions with relevance to nuclear structure as well 
as nuclear astrophysics [123]. 



5.2.1 Coulomb excitation 

Following multipole expansion of the electromagnetic field of a nucleus with charge Z2 as it 
evolves along a classical Rutherford trajectory, and with first order time-dependent pertur- 
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Ey [MeV] 

Fig. 14. Total number of virtual photons for the El multipolarity, "as seen" by a projectile passing 
by a lead target at impact parameters bmin = 12.3 fm and larger, for three typical bombarding 
energies. (Adapted from T. Glasmacher [123]). 



bation theory, the Coulomb excitation cross section is given by [120] 



dVt \dVt , V (2A + 1 

\ / cl TTA/t ^ 

where B{TTX,Ii — > //) is the reduced transition probabiUty of the projectile nucleus, ttA = 
El, E2, Ml, ... is the multipolarity of the excitation, and fi = —A, — A + 1, . . . , A. The orbital 
integrals S{'n'X,fi) contain the information on the dynamics of the reaction [124]. Inclusion 
of absorption effects in S{TTX,fi) due to the imaginary part of an optical nucleus-nucleus 
potential were worked out in Ref. [125]. These orbital integrals depend on the Lorentz factor 
7 = (1 — v^/c^)~^/^, with c being the speed of light, on the multipolarity vrA/i, and on the 
adiabacity parameter ^(6) = Ufib/jv < 1, where Ufi = {Ef — Ei) /h = E^/h is the excitation 
energy (in units of K) and b is the impact parameter. 

Because, the Coulomb excitation process is an external process, i.e., it occurs when the 
nucleons from one nucleus are outside the nuclear matter distribution from the other nucleus, 
the matrix elements for Coulomb excitation are the same as those for excitation by real 
photons (except for E{) Coulomb excitations, which are extremely small). Therefore, Coulomb 
excitation cross sections are directly related to the photonuclear cross sections by means of 
the equation [124] 

= Y. '^^^h (E.) + T. '^^<>. (E.) , (92) 

"-C/a; EX MX 

where cr^;^ {Ex) are the photonuclear cross sections for the multipolarity vrA and E^ is the 
excitation energy. 

The photonuclear cross sections are related to the reduced matrix elements, for the excitation 
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Fig. 15. Jupiter and comet Shoemaker-Levy 9, as imaged by the Hubble Space Telescope (HST), on 
May 18, 1994, when the giant planet was at a distance of 420 million miles (670 million km) from 
Earth. The gravitational interaction of Jupiter with the comet has broken it up into many pieces 
(picture from NASA). Much stronger tidal forces occur when nuclei come close to each other due 
to their mutual electromagnetic fields. The insert shows the dissociation of a fast nuclear projectile 
passing by the Coulomb field of a high-Z projectile. 



energy E^., through the relation [124] 

, (27r)3(A + l) (E,\''-' dB{nX,E,) 

^""^'^ = A[(2A + l)!!f l^J dE^ 

where dB/dE^ are the electromagnetic response functions, such that 



dB{nX,E,) 



dE,. 



For differential cross sections one obtains 



(94) 



where Q denotes to the solid scattering angle. 



Due to the use of high-energy projectiles in radioactive beam facilities, it is important to 
account for the strong absorption properly, as it occurs at small impact parameters. The 
eikonal formalism developed in Ref. [125] is appropriate for this purpose. The virtual photon 
numbers in Eq. (92) become 



A 



(2A + 1)!! 



2 



n^A(^x) = ^l" X. .X , A l^nXml'^ 9m{Ex) , (95) 

[Ztt)'^ (A + 1) ^ 
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and 

9m{E.) = 2n / db b X^(^) exp{-2 Mb)} , (96) 

where 0/(6) is the imaginary part of the eikonal phase 0(6), Eq. (87). The functions GT^\m{c/v) 
are given in Ref. [122]. 

In Figure 14 we show a calculation (with = E^) of the virtual photons for the El 
multipolarity, "as seen" by a projectile passing by a lead target at impact parameters equal 
to and exceeding 6 = 12.3 fm, for three typical bombarding energies. As the projectile energy 
increases, more virtual photons of large energy are available for the reaction. This increases 
the number of states accessed in the excitation process. 



5.2.2 Coulomb dissociation 

Coulomb dissociation is a process analogous to what happened to the comet Shoemaker-Levy 
as it disintegrated during its approximation to Jupiter in 1994 (see Figure 15). Approximately 
1.5 to 2.2 hours after closest approach, the comet (which was presumably a single body at 
the time) was broken apart by tidal forces into at least 21 pieces. The pieces continued to 
orbit Jupiter with a period of approximately 2 years. Due to gravitational forces from the 
Sun, which changed the orbits slightly on the next approach to Jupiter, the pieces impacted 
the planet. Much stronger tidal forces occur when nuclei come close to each other due to their 
mutual electromagnetic field. This leads to their dissociation, especially for weakly-bound 
nuclei (see insert to Figure 15). 

The idea behind the Coulomb dissociation method is relatively simple. The (differential, or 
angle-integrated) Coulomb breakup cross section for a + A^b-\-c-\-A follows from Eq. 
(91). It can be rewritten as 

= F-\E,- e- 0) . ^ ,^AE,). (97) 

where E^ is the energy transferred from the relative motion to the breakup, and cr^+a ^ ^j^J^E^) 
is the photo-dissociation cross section for the multipolarity ttA and photon energy E^. Time 
reversal allows one to deduce the radiative capture cross section b -\- c a -|- 7 from 

^ 6+c(-^7)' 



2(2i„ + 1) k 



2 



<+c-.a^,{E,) - (2^.^ + i)(2i, + l)^^^^+« - ^^^^ 
where k"^ — 2m},c{E^ — S) with S equal to the separation energy, and = E^/hc. 

This method was proposed in Ref. [126] and has been tested successfully in a number of 
reactions of interest to astrophysics. The most celebrated case is the reaction ^Be(p, 7)^6, first 
studied in Ref. [127], followed by numerous experiments in the last decade (see Section 6.2.3 
for the details of Coulomb dissociation experiments). 
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Ean (MeV) E,, (MeV) 

Fig. 16. Upper figure: Electric dipole response function for "He. The shaded area represents the 

experimental results from a Coulomb dissociation experiment [128]. The dashed and dotted lines 
correspond to results from three-body decay models from Refs. [129,130]. Lower figure: Measure- 
ment of two-body correlations in the three-body decay of ^He. The lower panels show the ratio 
between the measured a — n and n-n relative-energy spectra (upper panels) and the spectra simu- 
lated (histograms) according to standard phasespace distributions [128]. (Courtesy of T. Aumann). 



The two-neutron capture on ^He could perhaps play a role in the post-collapse phase in type- 
II supernovae. The bottleneck in this nucleosynthesis scenario is the formation of nuclei with 
A>9 from nucleons and a-particles. In principle, the reaction ^He(2n,7)^He could be rele- 
vant in bridging the instability gap at A = 5, although it is believed that this reaction cannot 
compete with the {an, 7) process in a type-II supernova scenario. Experiments with Coulomb 
dissociation have been used to study this question, as shown in the example presented in 
Figure 16. The upper part of the figure displays the electric dipole response function for ^He. 
The shaded areas represent the experimental results from a Coulomb dissociation experiment 
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[128]. The dashed and dotted hnes correspond to results from three-body decay models from 
Refs. [129,130]. In the lower figure we show the measurement of two-body correlations in 
the three-body decay of ''He. The lower panels display the ratio between the measured a-n 
and n-n relative-energy spectra (upper panels) and the spectra simulated (histograms) ac- 
cording to standard phasespace distributions [128] . Prom the analysis of this experiment it 
was found that 10% of the dissociation cross section proceeds via the formation of ^He. A 
rough estimate yields 1.6 mb MeV for the photoabsorption cross section for ^He(7, n)'''He, 
which agrees with theoretical calculations [131]. From this experiment one concludes that 
the cross sections for formation of ^He and ^He via one (two) neutron capture by ^He are 
not large enough to compete with the (an, 7) capture process (for more details, see ref. 
[132]). Nonetheless, this and the previously mentioned examples, show the relevance of the 
Coulomb dissociation method to assess some of the basic questions of relevance for nuclear 
astrophysics. 

5.2.3 Higher-order effects and nuclear dissociation 

Eq. (97) is based on first-order perturbation theory. It further assumes that the nuclear 
contribution to the breakup is small, or that it can be separated under certain experimental 
conditions. The contribution of the nuclear breakup has been examined by several authors 
(see, e.g. [125]). For example, has a small proton separation energy 140 keV). For 
such loosely-bound systems it was shown that multiple-step, or higher-order effects, are 
important [133,134]. These effects occur by means of continuum-continuum transitions. The 
role of higher multipolarities (e.g., E2 contributions [135] in the breakup reactions and the 
coupling to high-lying states) also has to be investigated carefully. Investigations related to 
the effects of relativity have attracted much theoretical interest as well [136-139]. 

5.3 Transfer reactions 

Transfer reactions, ^4(0, 6)S, are effective when a momentum matching exists between the 
transferred particle and the internal particles in the nucleus. Therefore, beam energies should 
be in the range of a few 10 MeV/nucleon [140]. One of the many advantages of using transfer 
reaction techniques as surrogates for direct measurements in nuclear astrophysics is to avoid 
the treatment of the screening problem [141]. 

5.3.1 Trojan horse method 

Low-energy reaction cross sections of astrophysical interest can be extracted directly from 
breakup reactions A-\- a ^ h -\- c-\- C hj means of the Trojan horse method, as proposed 
in Ref. [142] (see Figure 17). If the Fermi momentum of the particle x inside a — [h + x) 
compensates for the initial projectile velocity the low-energy reaction A + x = C + cis 
induced at very low (even vanishing) relative energy between A and x (for more details, see 
Refs. [143,144]). 
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Without loss of generality, we neglect the effects of spin, and write the cross section for 
A + X ^ c + C at an energy close to the threshold as 

aA+.^B+c= ^T.i'^l+IM:,]^ (99) 

'^x I 

where is the wave number in the incident channel, and Six is the matrix element for the 
1-th partial wave. 



In DWBA, the cross section for the breakup reaction A + a ^ b + c + C is 



2 



mambrric k^kc 



with 



Im 



2 



(100) 



Tin, = {xi ^YUi iHxl xtA.) , (101) 
where the x's denote the scattering waves and the interaction between b and x is denoted 
by Vbx- The internal projectile wavefunction corresponding to this potential is (p^x- 

The radial motion of particles x + A is governed by the function fi which asymptotically is 
given by 

fi"^ YT^ [Gi{ri, kxr) + iFi{r], kxv)] , 

KxT 

where Fi and Gi are the regular and irregular Coulomb wavef unctions, respectively. 

The threshold behavior of Ex for the reaction A-\- x = C + c is well known: since IjS^j;! ~ 
exp(— 27r?7). Thus, 

(TA+x^c+c - ^ E(2^ + 1) exp(-27rr;), (102) 

l^x I 

independent of the partial wave I. Thus, it follows that ua+x^b+c ~ (l/^D exp(— 27777). 
addition to the threshold behavior of Six, the three-body cross section given by Eq. (100) 
is also governed by the threshold behavior of fi{r). This is determined from the behavior 
of the irregular Coulomb function Gi{r],kxr) for /c^; — > (77 — > 00). The combined threshold 
behavior of fi{r), which for r — > 00 is given by 

fl. ~ ihr)'^' exp(7r7y) K,i+,{0: 

where Ki denotes the modified Bessel function of the second kind. The quantity ^ is inde- 
pendent of kx and is given by ^ = (Sr/os)^/^, where qb — jmZA'Zx^ is the Bohr length. 

Based on the arguments given above, the cross section for the reaction A^ a ^ b ^ C 



^ In a simple two-body reaction, T is related to the scattering amplitude, Eq. (85), by / = 



45 




Fig. 17. (a) Charged-particle reactions at low energies, A + x ^ c + C, are strongly hindered by 
Coulomb repulsion, (b) But, if particle x is brought into the reaction zone of nucleus A inside a 
projectile a = {b + x) with velocity Va, it can induce a reactions at low energies, corresponding to 

Vx ^ Va- V Fermi < Va- 



close to the threshold of A + x, with x "hidden" in a, is given by [142] 



dQhdQrdEi 



(E^. ^ 0) ^ const. (103) 



That is, the coincidence cross section tends to be a constant which, in general, will be different 
from zero. This is in striking contrast to the threshold behavior of the two particle reaction 
A + X = c + C . The strong barrier penetration effect on the charged particle reaction cross 
section is canceled completely by the behavior of the factor Tim for rj ^ oo. Thus, from 
a measurement of reaction A + a^b + c + C and a theoretical calculation of the factors 
Tim in Eq- (100), one extracts the matrix elements Sic needed for A + x = c + C . Basically, 
this technique extends the method of transfer reactions to continuum states. Very successful 
results using this technique have been reported by many authors, e.g. Ref. [149,150,141]. The 
problem with the method is that the x + A scattering is off-t he-energy shell. The initial- and 
final-state interactions should be taken into account to get a correct absolute value of the 
extracted astrophysical factor. 

The main characteristic of the Trojan horse method is the suppression of both Coulomb 
barrier and screening effects in the off-shell two-body cross section [143,144]. However, the 
quasifree A + x process can occur at very low, sub-Coulomb energies, even lower than the 
simple composition of projectile velocity -|- Fermi velocity, thanks to the role of the {A, x)- 
binding energy in compensating for the A + a relative motion. This is a different approach to 
the Trojan horse method compared to the original idea of Ref. [142]. In particular, in plane- 
wave impulse approximation (PWIA) the cross section of the A + a— >c + C + 6 three-body 
reaction can be factorized as given by 



d^ 




(104) 



dflcdflcdEc 

where {da/dQ)off is the off-energy-shell differential cross section for the two-body reaction 
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Fig. 18. S'-factor obtained with the Trojan horse method (fuh red dots) for the reaction -^^N(p, a)^^C 
(from Ref. [141]) using the ^H(-'^^N,a^^C)n reaction at Eiat = 60 MeV. The direct data from Refs. 
[145-147] are also shown as open symbols (circles, squares, and triangles, respectively). The red 
line represents a fit to the data. For comparison, a Breit-Wigner parametrization is also displayed 
as the black line [145]. 

A + X c + C and KF is a kinematical factor given by 



where B = A + x. In Eq. (104), ^{k^b) is the Fourier transform of the bound-state wave- 
function of a = 6 + a;. Of course, the plane-wave impulse approximation is an idealistic 
approximation and corrections to this method are certainly necessary. For example: (a) the 
plane wave might be replaced by a Coulomb, or distorted, wave in the final channel; (b) the 
initial and final channels might be treated differently; (c) higher-order processes might be 
relevant and so on. These corrections have been discussed, e.g. in Ref. [143,144]. 

Figure 18 shows the ^'-factor obtained with the Trojan horse method (full red dots) for the 
reaction ^^N(p,Q;)^^C (from Ref. [141]) using the '^}i{^^N,a'^'^C)n reaction at Eiab = 60 MeV. 
This reaction is important for nucleosynthesis in AGB stars. The direct data from Refs. 
[145-147] are also shown as open symbols (circles, squares, and triangles, respectively). The 
red line represents a fit to the data. For comparison, a Breit-Wigner parametrization is also 
shown as the black line [145]. 

Figure 19 shows in the left panel the bare Sb{E) factor data (filled circles) for the reaction 
^Li{d, a)a obtained with the trojan horse method [148]. It also shows the screened S{E) factor 
data (open diamonds) from direct experiments [101]. The dashed curve is a polynomial fit 
to S{E) and the solid curve includes the effects of electron screening with Ue = 340 eV. 




(105) 
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Fig. 19. Left: Bare Si,{E) factor data (filled circles) for the reaction ^Li(d, a)a from the trojan horse 
method [148] and screened S{E) factor data (open diamonds) from direct experiments [101]. The 
dashed curve is a polynomial fit to S{E) and the solid curve includes the effects of electron screening 
with Ue = 340 eV. Right: S{E) for ^Li{p, a)^He obtained with the trojan reaction ^Li(d, n)^He (solid 
circles) [152] compared to direct data (open triangles [153] and open circles [101]). The line shows 
the result of a second-order polynomial fit to the trojan horse reaction data [152]. 

On the right panel one sees S{E) for ^Li(p, a)^He obtained with the trojan horse reaction 
^Li{d, nYB.e (solid circles) [152] compared to direct data (open triangles [153] and open circles 
[101]). The line shows the result of a second-order polynomial fit to the trojan horse reaction 
data [152]. It is apparent that an independent measurement of the screening potential can 
be obtained in experiments with the trojan horse method [154]. The results agree with the 
direct measurements for f/g [154]. The method also shows that bare ^-factors can be obtained 
directly with this method. 



5.3.2 Asymptotic normalization coefficients 

The asymptotic normalization coefficient (ANC) method relies on the fact that the amplitude 
for the radiative capture cross section b + x—>-a + '~f is given by 

M^x = (/^.(rbx) |aA(rbx)| ^J^Hrbx)) , (106) 
as was described in previous sections. The overlap integral 

/^.(rbx) =< M^b, r^M.iUM^b) > (107) 

corresponds to the integration over the internal coordinates and of b and x, respec- 
tively. For low energies, the overlap integral I^^ is dominated by contributions from large 
rbx- Thus, what matters for the calculation of the matrix element A^ttA is the asymptotic 
value which, for charged particles is according to Eq. (26), /^^ ~ C^^ W_rit^^i/2{2Kbxrbx)/rbx, 
where C^^ is the ANC and W is the Whittaker function. This coefficient is the product of 
the spectroscopic factor and a normalization constant, which depends on the details of the 
wavefunction in the interior part of the potential. Thus, C^^ is the only unknown factor 
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Fig. 20. Left: Schematic diagram for a radiative capture reaction. Right: The upper hah of this 
diagram contains part of the information of the diagram on the left. If the cohision is peripheral, 
only the tail of the bound-state wavefunction of p+B is involved, allowing the extraction of its 
ANC by a DWBA analysis of the experimental data. 

needed to calculate the direct capture cross section. 

The normalization coefficients can be found from: 1) analysis of classical nuclear reactions 
such as elastic scattering by extrapolation of the experimental scattering phase shifts to the 
bound-state pole in the energy plane, or 2) peripheral transfer reactions whose amplitudes 
contain the same overlap function as the amplitude of the corresponding astrophysical ra- 
diative capture cross section. This is shown schematically in Figure 20, where the left panel 
shows a schematic diagram for a radiative capture reaction. The upper half of the diagram 
on the right contains part of the information of the diagram on the left. If the collision is 
peripheral, only the tail of the bound-state wavefunction of p+B is involved, allowing the 
extraction of its ANC by an DWBA analysis of the experimental data. This method was 
proposed in Ref. [43] and has been used with success for many reactions of astrophysical 
interest, as discussed, e.g. in Ref. [150] and mentioned in Section 6.2.6. 



To illustrate this technique, let us consider the proton transfer reaction A{d,a)B, where 
d = a + p, B = A + p. Using the asymptotic form of the overlap integral, the DWBA cross 
section is given by 



da 

dn 
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a 



(108) 



where a is the reduced cross section, not depending on the nuclear structure. Pap (Pap) 
are the asymptotic normalization of the shell-model bound-state proton wave functions in 
nucleus d{B), which are related to the corresponding ANC's of the overlap function as 
i^ap)^ = ^apPap- Here Sap is the spectroscopic factor. Suppose the reaction A{d,a)B is 
peripheral. Then each of the bound-state wavefunctions entering a can be approximated by 
its asymptotic form and a oc ^ApPap- Hence 
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f^Apf^ap 



(109) 



is independent of P\p and P^p- Thus for surface-dominated reactions, the DWBA cross section 
is actually parameterized in terms of the product of the square of the ANC's of the initial 
and the final nuclei (C'Ap)^(^ap)^ rather than spectroscopic factors. This effectively removes 
the sensitivity in the extracted parameters to the internal structure of the nucleus. 
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Now consider the elastic a + p scattering amplitude. It has a pole in the momentum plane 
[151] 

^'^^^y"^- 2ik 2ikok-ko' ^ ^ 

where, for the bound state of dior kg — in and for a resonance for kg — ku, with ku — kr—i ki, 
and where 5';^^^ is the elastic S'-matrix. The residue in the pole Wi^j^ is 

= -(-1)'^ ^e^"""' (C'^pidJ^ ^0 = ^ (111) 

Wujd = -(-1)''^ i (C^pidUmf^ ^0 = kn. (112) 
For narrow resonances, kj « kr, 

(C,(«))' = (-1)^"* ^e-'' e'^'^ui^^^ r,,,,. (113) 

Here, 7]r is the Sommerfeld parameter for the resonance at momentum kr, Si^j^{kr) is the 
potential (non- resonant) scattering phase shift taken at the momentum kr [151]. 

For elastic scattering close to the threshold, the bound-state pole behaves (see Section 3.3) 
as 

/..W^=-^e-'^---' ^^J:,r^^, . (114) 

where = 2P/^(£')7^, and Pi^{E) is the Coulomb barrier penetrability. In this equation, 
0/^ is the hard-sphere phase shift in the partial wave Id and cr^^ = tan~ (r^rf/n) ,7^ is the 

n=l 

reduced width: 

2_ 1 ^-r,,,;rf+l/2(2/^re) 2 

- W; V^aplajdivc)) ' i^^^J 

^H'ap I c 

where Tc is the channel radius. 

Thus, close to threshold and in the presence of a bound state, the scattering amplitude 
behaves as the high-energy tail of the resonance located at energy E — —Sc, also called 
a "subthreshold" resonance. It is not a resonance because the real resonance is located at 
complex energies on the second energy sheet, while the subthreshold resonance is a bound 
state located on the first energy sheet at negative energy. At negative energies (positive 
imaginary momenta), Eq. (114) reduces to Eq. (110). Thus, the elastic scattering amplitude 
(phase shift) offers another possibihty to determine the ANC by extrapolating Eq. (110) to 
the bound-state pole [155]. 



5.4 Nucleoli knockout reactions 



Modern shell-model calculations are now able to include the effects of residual interactions 
between pairs of micleons, using forces that reproduce the measured masses, charge radii 
and low-lying excited states of a large number of nuclei. For very exotic nuclei, the small 
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Fig. 21. Coordinates used to describe knockout reactions in the text. (Adapted from J. A. Tostevin). 



additional stability that comes with the filling of a particular orbital can have profound effects 
upon their existence as bound systems, their lifetime and structure. Thus, the verification 
of the ordering, spacing and the occupancy of orbitals is essential in assessing how exotic 
nuclei evolve in the presence of large neutron or proton excess and to what extent theories 
have predictive power. Such spectroscopy of the single-particle structure in short-lived nuclei 
typically uses direct nuclear reactions. 



Heavy-ion induced single-nucleon knockout reactions - performed at intermediate energies 
and in inverse kinematics - have become a specific and quantitative tool for studying the 
location and occupancy of single-particle states and correlation effects in the nuclear many- 
body system, as discussed in Refs. [156-159]. In a peripheral, sudden collision of the fast- 
moving projectile with mass A with a light target (typically Be or C), a single nucleon 
is removed from the projectile, producing projectile-like residues with mass A — 1 in the 
exit channel [158]. Typically, the final state of the target and that of the struck nucleon 
are not observed, but instead the energy of the final state of the residue can be identified 
by measuring the decay 7-rays emitted in flight by the excited projectile-hke residues (see 
Section 6.2.6). 



The momentum distributions of the projectile- like residues in one-nucleon knockout are a 
measure of the spatial extent of the wavefunction of the struck nucleon, while the cross 
section for the nucleon removal scales with the occupation amplitude, or probability (spec- 
troscopic factor), for the given single-particle configuration in the projectile ground state. 
The longitudinal momentum distributions are given by (see, e.g., Refs. [160,161]) 



rfOstr 



(27r)'2/ + 



1 ^io 



1 - \Sn {hn)[ 



dz exp [-ikzz] ipin, (r) 



(116) 



where represents the longitudinal component of (final momentum of the core of the 
projectile nucleus) and (C^S) is the spectroscopic factor, and ipim (r) is the wavefunction 
of the core plus (valence) nucleon system (c -f n) in a state with single-particle angular 
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momentum I, m. In this equation, r = (p, z, 0) = r„ — Tc, so that 



bc= \p-hn\ = J p'^ + bl - 2p bnCOS(p = Jr'^sm^ O + b^- 2rsm9 6„cos0. (117) 



S{b) arc the S-matrices for core-target and nucleon-target scattering obtained from the 
nuclear ground-state densities and the nucleon-nucleon cross sections by the relation [44] 

S{b) = exp [i(p{b)] , with 0jv(6) = / dq q pp (q) pt {q) /nn (q) Jo {qb) , (118) 

KNN Jo 

where pp^t (?) is the Fourier transform of the nuclear density of the projectile (nucleon or core) 
and the target nucleus, and fNN (q) is the high-energy nucleon-nucleon scattering amplitude 
at forward angles, which can be parametrized by [162] 

fNN (q) = + '^^'^) {-PNNq'^) , (119) 

where kNN is the nucleon-nucleon relative momentum, aNN is the nucleon-nucleon total cross 
section, and aNN and (3nn are parameters fitted to (free) nucleon-nucleon scattering data. 

The first term in the integral in Eq. (116), 1 — represents the probability for the 

knockout of the nucleon from its location at whereas the second integral carries the term 
jS'cp which is the probability of core survival at impact parameter be- At high energies, 
the (S'-matrices do not depend on the longitudinal direction. That is why the bound-state 
wavefunction is probed by the longitudinal Fourier transform in the last integral of Eq. (116). 
These results arise naturally by using eikonal scattering waves. For the transverse momentum 
distributions, the same formalism yields 
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Extensions of the nucleon knockout formalism include the treatment of final-state interac- 
tions were discussed in Ref. [160] where one has shown that the inclusion of the Coulomb 
final-state interactions are of relevance. They can be done by just adding the Coulomb phase 
4> — 4>N + 4>c in the eikonal phase described above [160]. Other higher-order effects have 
been included [163] and a theory for two-nucleon knockout [164-167] has been developed. 
Knockout reactions represent a particular case for which higher projectiles energies allow a 
simpler theoretical treatment of the reaction mechanism, due to the simplicity of the eikonal 
scattering waves and the assumption of a single-step process. 

Many new experimental approaches based on nucleon knockout reactions have been devel- 
oped and shown to reduce the uncertainties in astrophysical (rapid proton capture) rp-process 
calculations by the provision of nuclear structure data. One utilizes neutron removal from 
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a rare-isotope beam to populate nuclear states of interest in the knockout residue. In the 
first case studied with astrophysical background, e.g., in Ref. [f68,f69], excited states in "^'^Ar 
were measured with uncertainties of several keV. The 2-orders-of-magnitude improvement 
in the uncertainty of the level energies resulted in a 3-orders-of-magnitude improvement in 
the uncertainty of the calculated ^^Cl(p,7)^^Ar rate that is critical to the modeling of the 
rp process (see Section 6.2.6 for experimental details). This approach has the potential to 
measure key properties of almost all interesting nuclei on the rp-process path. 



5.5 Quasifree (p,pN) reactions 



Quasifree (p.pN) reactions (A^ = proton or neutron) represent one of the most direct ways 
to measure single-particle properties in nuclei. They have been used and extensively studied 
for over 4 decades (see [170,171] for reviews). In quasifree {p,pN) scattering, an incident 
proton of intermediate energy (200-1000 MeV) knocks out a bound nucleon. The only violent 
interaction of this process occurs between the incident particle and the ejected nucleon. The 
wavefunctions of incoming and outgoing nucleons are distorted while traversing the nucleus. 
From the measured energies and momenta of the emerging nucleons, direct information 
on single-particle separation energy spectra and nucleon momentum distributions can be 
obtained. Over the last four decades, quasifree scattering experiments have been performed 
with this basic purpose. 

A popular framework used in the analysis of quasi-free scattering is the so-called "impulse 
approximation" , which yields for the quasifree cross section 



K-^\Fm\ (121) 



dT^dflpd^N dfl 

where X is a kinematic factor, |F(Q)p is the momentum distribution of the knocked-out 
nucleon A^ in the nucleus and dcpN / dVt the free p-N cross section. In the knockout formalism, 
the off-shell p-N t-matrix is required, and the factorized form that appears in Eq. (121) is 
valid only if off-shell effects are not very important. In proton-induced knockout reactions, 
this is probably adequate since the energy variation of the free nucleon-nucleon cross section 
is small. As expected, there is an evident similarity between this equation and Eq. (104) 
used in the analysis of transfer reactions. 

Improvements of the above formulation can be easily made, assuming that (on-shell or 
off-shell) multiple-scattering effects are small. In the distorted-wave impulse approximation 
(DWIA), the scattering matrix for the A{p,pN)B reaction is given by 

Tv,pN = (xic^^xLt M Xk>o) (122) 

where xk'^ (Xk'^ ) distorted wave for an outgoing proton (knocked-out nucleon) in the 

presence of the residual nucleus B, Xk^ is the distorted wave for an incoming proton in the 
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presence of the target nucleus A, and 0o is the bound-state wavefunction of the knockcd-out 
nucleon; 0o inchidcs the spectroscopic amphtude. The p-N scattering matrix is denoted by 
TpN. The matrix element given by eq. (122) can be written as 

Tp,pN^^ J dXsd''r'j,Bd\AdhNBd^Ppd^PNd^Ppd^PNTiPp,PN;'P'p,'P'N) 
X xic^^*(P;)xici*(P'iv)xS(Pp)^o(r^B)5 + P'^ - P^ - P^v) , 

(123) 

where Xk(P) are the Fourier transforms of Xk(i'), which are normalized so that / d^rxk(r)Xk'(i') = 

The in-medium effects due to the propagation of the incoming proton and the outgoing 
particles in the nuclear interior can be included in the distorted waves Xi(r). The effective 
interaction TpN must also include medium and energy dependence effects. Spin dependence 
of the p-N t-matrix and of the scattering and bound states also needs to be included. Thus, a 
deviation from the simple formulation based on Eq. (122) is expected. By using inverse kine- 
matics, quasifree scattering will become a major tool for the investigation of the properties 
of rare isotopes with relevance for nuclear astrophysics [172]. 

5.6 Charge- exchange reactions 

Charge exchange induced in {p, n) reactions is often used to obtain values of Gamow- Teller 
matrix elements, B{GT), which cannot be extracted from /3-decay experiments. This ap- 
proach relics on the similarity in spin-isospin space of charge-exchange reactions and /3-decay 
operators. As a result of this similarity, the cross section a{p, n) at small momentum transfer 
q is closely proportional to B{GT) for strong transitions [173], i.e., 

^{q^O)^KNn\Jar\'B{a), (124) 

where K is a kinematical factor, Nd is a distortion factor (accounting for initial and final 
state interactions), Jg-r is the Fourier transform of the effective miclcon- nucleon interaction, 
and B{a = F, GT) is the reduced transition probability for non-spin-flip, 

B{F) ^ -^^\{f\\Y.Tt\\i)\\ 

* k 

and spin-flip, 

5(GT) = --^|(/||5:^fcrf^llOP, 

* k 

transitions. The Fourier transform of the effective interaction yields largest values for the 
I Jo-Tp component at energies around 100-200 MeV. This indicates that the energy range of 
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Fig. 22. Left: Electron captures by nuclei are dominated by Gamow- Teller (GT+) transitions. Right: 
The matrix elements for such transitions can be probed by using rare-isotope beams interaction 
with protons in inverse kinematics, as shown schematically for the reaction ^^Ni(p, n)^^Cu. (Adapted 
from work of D. Frekers and R. Zegers). 

100-200 MeV/nucleon will be appropriate to extract GT (F) matrix elements from studies 
charge-exchange experiments in inverse kinematics involving rare-isotope beams. The con- 
dition g ~ is met for very small scattering angles, such that 6 <^ 1/kR, where R is the 
nuclear radius and k is the projectile wavenumber. 

Eq. (124) is easily understood in the plane-wave Born-approximation. Then the charge- 
exchange matrix element is given by [176] 



Mexch{<i) 



(125) 



where q is the momentum transfer, "^^^'P intrinsic wavefunctions of nuclei a and b 

for the initial and final states, are the nucleon coordinates within a and b, and Vexch is 
the part of the nucleon-nucleon interaction responsible for charge exchange, which contains 
spin and isospin operators. For forward scattering and low- momentum transfers, q ~ 0, and 
the matrix element (125) becomes 



Mexchiq ~ 0) ~ t;2L(q ~ 0) Ma{F, GT) M,{F, GT) 



(126) 



where vf^^^^ is the spinless part of the interaction, and Aiexch{F, GT) = (^"^^a^b 11(1 cr)!"! |\E' 
are Fermi or Gamow- Teller (GT) matrix elements for the nuclear transition. The result above 
is also valid if, instead of plane waves, one uses the eikonal scattering waves for the nuclei. 
One can thus conclude that the ability to extract information on Fermi or Gamow- Teller 
transition densities in charge-exchange reactions in a simple way depends on the validity 
of the low-momentum transfer assumption in high-energy collisions. Recently, also for the 
(^He, t) reaction at 420 MeV, the proportionality between cross section and Gamow- Teller 
strength was shown to follow simple trends as function of the mass number [174]. 

Another assumption in Eq. (124), the validity of one-step processes, was proven to work 
rather well for (p, n) reactions (with a few exceptions). For heavy- ion induced charge- 
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exchange reactions the formula might not work so well as has been investigated in Refs. 
[175,176]. In Ref. [175] it was shown that multi-step processes involving the physical ex- 
change of a proton and a neutron can still play an important role up to bombarding energies 
of 100 MeV/nucleon. Refs. [176] use the isospin terms of the effective interaction to show that 
deviations from the Taddeucci formula are common under many circumstances. As shown 
in Ref. [177], for important GT transitions whose strengths are a small fraction of the sum 
rule, the direct relationship between (t(p, n) and B{GT) values also can fail. Similar discrep- 
ancies have been observed [178] for reactions on some odd-A nuclei including ^^C, ^^N, ^^Cl, 
and ^^K and for charge exchange induced by heavy ions [179]. It is still an open question if 
Taddeucci 's formula is vahd in general. 

Electron capture by nuclei in p/-shell plays a pivotal role in the deleptonization of a massive 
star prior to core-collapse [180]. During the period of silicon burning, supernova collapse oc- 
curs due to a competition of gravity and the weak interaction, with electron captures on nuclei 
and protons and /3-decay playing crucial roles. Weak-interaction processes become important 
when nuclei with masses A ~ 55 — 60 are most abundant in the supernova core. Weak inter- 
actions change and electron capture dominates, the value is successively reduced from 
its initial value ~ 0.5. Electron capture (EC) yields more neutron-rich and heavier nuclei, as 
nuclei with decreasing Z/A ratios are more bound for heavier nuclei. For densities p < 10^^ 
g/cm^, weak-interaction processes are dominated by Gamow- Teller and sometimes by Fermi 
transitions. This is shown schematically on the left panel of figure 22. Fuller, Fowler and New- 
man [181,182,22] made systematic estimates of EC-rates in stellar environments. But their 
calculations obtain only the centroid of the Gamow- Teller response function. More recently, 
i?(GT)-distributions have been obtained with modern shell-model calculations [183-186]. 
Some marked deviations from the previous rates [181,182,22] have emerged, e.g. Ye increases 
to about 0.445 instead of the value of 0.43 found previously [181,182,22]. 

Most of such theoretical developments need support and validation from experiments. A 
satisfactory understanding of how weak interaction rates influence stellar evolution will rely 
on experiments using charge-exchange reactions such as the traditional (^He,t) reaction and 
novel inverse-kinematics {p, n) reactions on rare-isotope projectile beams (right panel of Fig- 
ure 22). Reactions with unstable nuclei will provide crucial information on hitherto unknown 
B{F) and B{GT) values needed for astrophysical purposes and the validation of the handling 
of weak interactions by large-scale shell-model calculations. 

5.7 Central collisions 

At intermediate energies of Eiab ~ 100 — 1000 MeV/nucleon, the nucleons and the products 
of their collisions can be described individually and their propagation can be calculated by 
semiclassical equations. Hadronic transport theories have been quite successful in describing 
a multitude of measured particle spectra. 

The nuclear equation of state (EOS), e(p, 5), expresses the energy per nuclcon of nuclear 
matter as a function of the nucleon density p and the relative neutron excess 6. It is a funda- 
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mental quantity in theories of neutron stars and supernova explosions. The main measured 
quantities, which can provide information about the EOS are, for example, binding energies 
and other data of finite nuclei. As the finite nuclei are in states near the standard nuclear 
matter state with normal nucleon density po and zero neutron excess, 5 — 0, our knowledge 
about the EOS can be confirmed experimentally only in a small region around p ~ po and 
5^0. With very neutron-rich nuclei and energetic heavy-ion collisions, the nuclear EOS can 
be tested well beyond this region. 

If one assumes that the system of nucleons forms a dilute system (or gas) of particles, such 
that the total volume of the gas particles is small compared to the volume available to the 
gas, then na^ <^ 1, where n is the number density of particles and a is the (interaction) 
radius of a particle. Since the particles in a neutral gas do not have long-range forces like 
the particles in a plasma, they are assumed to interact only when they collide, i.e., when the 
separation between two particles is not much larger than 2a. The term "collision" normally 
means the interaction between two such nearby particles. A particle moves in a straight 
line between collisions. The average distance travelled by a particle between two collisions 
is known as the mean free path. The mean free path depends on the cross section a ~ a^, 
and is given by A = l/ncx. One consequence of the requirement that the gas be dilute is that 
A » a. If the gas is dilute, the probabihty of three-body coUisions is much lower than that 
for two-body coUisions and they can be neglected. 

Assuming that these conditions are valid, a practical equation describing nucleus-nucleus 
collisions in terms of nucleon-nucleon collisions can be deduced. A popular transport equation 
is the Boltzmann-Uehling-Uhlenbeck (BUU) equation, 

X {f[f^ [1 - /i] [1 - /2] - /1/2 [1 - /{] [1 -m, 

(127) 

where / is defined so that, if dN is the number of particles in the volume element d^r 
whose momenta fall in the momentum element d^p at time t, then the distribution function 
/ (r, p, t) is given by dN — f (r, p, t) d^rd^p. The left-hand side of the above equation is due 
to the variations of / by means of a mean field U. When the right-hand side is taken as zero, 
the resulting equation is known as the Vlasov equation. The right-hand side of Eq. (127) 
is the collisional term, which accounts for binary collisions between nucleons by using the 
nucleon-nucleon cross section aNN- This equation incorporates the Pauli principle through 
the (1 — /)-terms to avoid scattering into occupied states. 

Eq. (127) needs as basic ingredients the mean field U and the cross section a^N- Because 
these two quantities are related to each other, one should in principle derive them in a self- 
consistent microscopic approach, as the Brueckner theory, for example. However, in practice 
the simulations are often done with a phenomenological mean field and free nuclear cross 
sections. The most commonly used mean field is of Skyrme-type, eventually with a momen- 
tum dependent part [187]. Another important ingredient in the transport theory calculations 
is the compressibility K of nuclear matter, which refers to the second derivative of the com- 
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Fig. 23. Matter distribution at several time stages of a nearly central nucleus-nucleus collision. 
(Courtesy of P. Danielewicz). 

pressional energy E with respect to the density: 

- (§) . (128) 

This is an important quantity, e.g., for nuclear astrophysics. In fact, the mechanism of su- 
pernova explosions is strongly dependent on the value K. Supernova models might or might 
not lead to explosions depending on the value of K. The central collisions of heavy nuclei are 
one of the few probes of this quantity in the laboratory. The dependence of the calculations 
on K follow from the dependence of the mean field potential U {U ^ E/A+ kinetic energy 
terms) on the particle density p. 

Following an initial interpenetration of projectile and target densities, the A^A^ collisions 
begin to thermalize matter in the overlap region making the momentum distribution there 
centered at zero momentum in the cm. system. The density in the overlap region rises above 
normal and a disk of excited and compressed matter forms at the center of the system. More 
and more matter dives into the region with compressed matter that begins to expand in 
transverse directions. At late stages, when the whole matter is excited, transverse expansion 
dominates. A further view of the situation is illustrated in Figure 23. The measurement of 
the matter distribution in these collisions and the comparison with transport theories allows 
one to deduce the incompressibility of nuclear matter. With neutron-rich nuclei, one will 
also able to extract the EOS dependence on the asymmetry properties of nuclear matter. 
Measuring other observables, such as particle production and their kinematic properties, also 
have a valuable contribution from the use of neutron-rich projectiles. 

Using a transport theory based on the BUU equation together with available experimental 
data, Ref. [188] determined that maximum pressures achieve in central collisions in the range 
of P = 80 to 130 MeV/fm^ in collisions as 2 GeV/nucleon. This corresponds to 1.3 x 10^"^ to 
2.1 X 10^"^ Pa). A similar analysis for collisions at 6 GeV/nucleon yields the respective values 
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Fig. 24. Pressure for neutron matter as a function of density. The different curves are based on 
theoretical models. The shadow bands display the range of possibilities for the EOS based on the 
analysis of experiments using a soft (lower shadow region) and stiff (upper shadow region) EOS. 
(Adapted from Ref. [188]). 

of P = 210 to 350 MeV/fm^ (3.4 x 10^^ to 5.6 x 10^^ Pa, respectively). These correspond 
to pressures 23 orders of magnitude larger than the maximum pressure ever observed in the 
laboratory, being 19 orders of magnitude larger than pressures within the core of the Sun. 
They are, in fact, only comparable to pressures within neutron stars. The analysis seems to 
be consistent with bounds for K of Eq. (128) within the range K = 170 — 380 MeV [188]. 
Figure 24 shows the pressure for neutron matter as a function of the density. The different 
curves are based on theoretical models. The shadow bands display the range of possibilities 
for the EOS based on the analysis of experiments using a soft (lower shadow region) and 
stiff (upper shadow region) EOS (adapted from Ref. [188]). Rare-isotope beams have not yet 
been used for experimental work along these lines. 



5.8 Electron scattering on radioactive beams 



Electron vs. radioactive-beam colliders have been proposed in the framework of the physics 
programs of few rare-isotope beam facilities [189,190]. Since the reaction mechanism for 
electron scattering is very well-known, the extraction of valuable physical quantities such 
as charge distributions, transition densities, and nuclear response functions can be obtained 
with good accuracy. 

In the plane- wave Born approximation (PWBA), the relation between the charge density 
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and the cross section is given by 




PWBA 



1 + {2E/Ma) sin^ {9/2) 



Fchiq) 



(129) 



where ctm = (e^/4i^^) cos^ (^/2) sin^^ (^/2) is the Mott cross section, the term in the denom- 
inator is a recoil correction, E is the electron total energy. Ma is the mass of the nucleus 
and 9 is the scattering angle. 



The charge form factor Fch (q) for a spherical mass distribution is given by 



where q = 2ksm{9/2) is the momentum transfer, hk is the electron momentum, and E = 



t9k?c^ + m^c^. The low-momentum expansion of Eq. (130) yields the leading terms 



Thus, a measurement at low-momentum transfer yields a direct assessment of the root mean 

squared radius of the charge distribution, (r^^^)^^^. As more details of the charge distribution 
are probed, more terms of this series are needed and, for a precise description of it, the form 
factor dependence for large momenta q is needed as well. 



The systematics of the charge-density distributions, with the inclusion of nuclei having ex- 
treme proton-neutron asymmetry, forms a basis for investigations addressing both the struc- 
ture of nuclei and the properties of bulk nuclear matter. Information about the nuclear matter 
incompressibility [191,192], the nuclear symmetry energy [193], and the slope of the neutron 
matter equation of state in its dependence on density [194] can be extracted. Despite the 
importance of the nuclear equation of state (EOS), some of its features remain fairly obscure. 
The electron-scattering studies will allow new constraints on the EOS. 



Inelastic electron scattering is a powerful tool for studying properties of excited states of 
nuclei, in particular their spins, parities, and the strengths and structure of the transition 
operators connecting ground and excited states (e.g. Refs. [195,196]). Electron scattering is 
the only method which can be used to determine the detailed spatial distributions of the tran- 
sition densities for a variety of single-particle and collective transitions. These investigations 
provide a stringent test of nuclear many-body wavef unctions. 



In the plane- wave Born approximation (PWBA), the cross section for inelastic electron 
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scattering is given by [195,196] 
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where Jj (J/) is the initial (final) angular momentum of the nucleus, (E.p) and (£", p') are 
the initial and final energy and momentum of the electron, and (5'o,q) = {AE / he, Ap/h) 
is the energy and momentum transfer in the reaction. Fij {q; UL) are form factors for mo- 
mentum transfer q and for Coulomb (C), electric (E) and magnetic (M) multipolarities, 
U = C,E, M, respectively. Only for small momentum transfers and forward scattering, the 
matrix elements in electron scattering are the same as those for real photons and Coulomb 
excitation. The later are fixed by the energy and momentum transfer relation |Ap| = AE/c. 
Electron scattering allows probing the momentum and energy transfer response indepen- 
dently. 



Due to its strong selectivity, collective and strong single-particle excitations can be studied 
particularly well in electron scattering. Electric and magnetic giant multipole resonances are 
of special interest. Several of them have been discovered and studied using electron scattering 
(see Ref. [197] and references therein). Some of these giant resonances are related, e.g., to the 
different charge and matter radii of nuclei, quantities that are expected to vary strongly in 
dripline nuclei. Thus, the difference in radii of the neutron and proton density distributions 
is accessible via studies of the excitation of giant dipole resonances (GDR) or spin-dipole 
resonances. The cross section of these processes depends strongly on the relative neutron-skin 
thickness [198]. 

The reaction theory of electron scattering is very well understood. Even in situations for 
which higher order processes, such as radiative effects, are of importance, the theory can 
be handled with precision. That is why electron scattering has the potential to become a 
precious tool for studies with rare-isotope beams in the future. 

An example of calculated elastic electron-scattering properties for the scattering off exotic 
nuclei is given in Figure 25, where the charge form factor squared for elastic electron scat- 
tering off tin isotopes is displayed as a function of the momentum transfer [199]. The two 
curves are for the extreme values of the asymmetry parameter 5 — {N — Z)/A, that is 5 = 
(A^ = Z = 50), and 6 = 4/9 {N = 90). The curves form an envelope around other curves 
with intermediate values of S. The sensitiveness of the form factor obtained in electron scat- 
tering on the nuclear charge distribution properties is evident and depends essentially on 
the neutron skin. The size of the neutron skin will provide data on the volume symmetry 
energy which was used as input to get the estimates presented in Figure 24. Recent studies 
of electron scattering off exotic nuclei have been pubhshed in Refs. [201-205] . 

The electromagnetic response of light nuclei, leading to their dissociation, has a direct con- 
nection with the nuclear physics needed in several astrophysical sites. In fact, it has been 
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Fig. 25. Calculated charge form factor squared for elastic electron scattering off tin isotopes as a 
function of the momentum transfer. The two curves are for the extreme values of the asymmetry 
parameter S = {N - Z) /A, that is 6 = {N = Z = 50), and S = 4/9 {N = 90). The curves form 
an envelope around other curves with intermediate values of 6. (From Ref. [199]). 

shown [204] that the existence of pygmy resonances might have important implications for 
theoretical predictions of radiative neutron capture rates in the r-process nucleosynthesis, 
and consequently to the calculated elemental abundance distribution in the universe. Re- 
cently, the implications of inelastic electron scattering for studying the properties of light, 
neutron- rich nuclei have been discussed in Ref. [205]. 



6 Nuclear-astrophysics experiments with rare-isotope beams 

6. 1 Production of exotic ion beams 

The first challenge of any experiment with short-lived "exotic" nuclei is their production. 
Today, a broad range of rare isotopes is available for experiments in form of ion beams. Two 
main production and separation mechanisms have emerged as the workhorse techniques in 
rare-isotope beam production and are employed in nuclear physics laboratories around the 
world: 

• Beams of short-lived nuclei are produced, separated "in-fiight" via fragment separators 
and are directly used for experiments (in-fiight separation). 
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• Exotic nuclei are produced, stopped, ionized, and reaccelerated (isotope separation on-line 
-ISOL). 

A third approach for the production of rare-isotope beams is the batch mode, where a "batch" 
of radioactive material is made into a beam directly in an ion source. This technique is limited 
to nuclei that have a sufficiently long half-life and that can be produced in a chemically 
suitable form. This method was successfully used to produce exotic beams of ''Be [212,213], 
[214], [215], [216,217], ^^Ti [218], ^^Co [219], ^^Ni [219] and ^^Kr [220]. 

In the following subsections we will summarize in-flight separation and ISOL techniques and 
the underlying production and reaction mechanisms. 

6.1.1 Production of exotic nuclei - In-flight separation 

A variety of reaction mechanisms employing various projectile-target combinations over a 
broad range of collision energies are used to produce short-lived neutron-rich or proton-rich 
nuclei in-flight. Peripheral collisions at and above the Coulomb barrier create isotopes in 
the proximity of the target and projectile nuclei via nucleon transfer and nucleon exchange. 
The highest yield for direct one-nucleon transfer or exchange reactions is achieved when the 
momenta of the transfered or exchanged nucleons are matched in the initial and final state 
(5 - 30 MeV/nucleon energy range) [221]. Using stable beams and stable targets, one-nucleon 
transfers cannot lead very far from stability. Typical production reactions are d{^^0, ^^F)n, 
p(^Li, ^Be)n, and ^He(^Li, ^B)n induced by intense, low-energy beams of ^^O or ®'^Li im- 
pinging on light targets. The choice of inverse kinematics yields a forward focussing of the 
projectile-like reaction residues, which can then be guided to a secondary reaction target 
or detector setup via ion-optical transport systems like superconducting solenoids [222,223] 
or bending magnets [224]. This low-energy in-flight technique is presently being used at the 
TwinSol facility at the University of Notre Dame (ND) [225], at Argonne National Labo- 
ratory (ANL) [224], at Texas A&M University (TAMU) [226], at CRIB [227] at CNS/Tokyo 
(Japan) and at the RESOLUT facility at Florida State University (FSU). 

In the high-energy collision of a heavy-ion beam with a thin target (typically at energies ex- 
ceeding 50 McV/miclcon), exotic nuclei are produced via projectile fragmentation. At high 
energy, nucleons are removed from the projectile and the non-interacting part of the projec- 
tile nucleus proceeds as "spectator" at essentially the initial beam velocity and close to 0° 
with narrow angular and linear momentum distributions [228,229]. This spectator or prefrag- 
ment subsequently undergoes statistical deexcitation leading to the observed products [230]. 
The mass, charge and momentum distributions of the prcfragments have been described 
with macroscopic as well as microscopic models [231]. The macroscopic abrasion-ablation 
model [232,233] uses the geometric abrasion, with the size of the prefragment calculated from 
the geometry of the overlap zone in the projectile-target collision. Statistical models can pro- 
vide the proton-to-neutron ratio, assuming that both nucleon species are abraded indepen- 
dently of each other. Microscopic approaches employ intranuclear cascade models [234,235], 
where the interaction of target and projectile is formulated in terms of elastic and inelastic 
nucleon-nucleon scattering. The deexcitation is incorporated with evaporation or transport 



63 



models. All elements up to uranium can be produced in projectile fragmentation. 

A fragment separator behind the production target separates the wanted exotic nuclei from 
the primary beam and other reaction residues by a combination of magnetic elements 
(see [236,237] for reviews). A selection according to the ion's magnetic rigidity separates 
particles with the same momentum-to-charge ratio. A wedge-shaped degrader can remove 
the degeneracy as the momentum of each species is systematically shifted by the energy 
loss encountered in the degrader material [238,239]. This in-flight separation is chemistry- 
independent and a limitation on the lifetime of the exotic nuclei that can be produced 
for experiments is solely given by the ion's time-of-flight through the separator and to the 
experimental setup (typically less than 1 fis). The available exotic beams typically have ener- 
gies exceeding 30 MeV/nucleon and event-by-event particle identification allows for sensitive 
and efficient experiments also with cocktail beams. The drawbacks of secondary beams pro- 
duced by projectile fragmentation are the beam energies that are much higher than those 
required by certain experiments, a poor transversal and longitudinal beam emittance orig- 
inating from the fragmentation process and, in some cases, limited beam purity. Projec- 
tile fragmentation as a production mechanism was pioneered at the BEVALAC in Berke- 
ley [240,241] and has been the main technique for the production of fast rare-isotope beams 
at GANIL (Prance) [242], GSI (Germany) [243], NSCL/MSU (US) [244], and RIBP/RIKEN 
(Japan) [245]. 

Although projectile fragmentation is extensively used to produce neutron-rich nuclei, the 
maximum production yield is realized for neutron-deficient species. Nuclear fission is a pro- 
cess known to be a source of very neutron-rich nuclei [246]. When a fissioning nucleus is 
moving at energies large compared to the fission recoil, the fission products can be collected 
and separated using in-flight techniques similar to those described for projectile fragmen- 
tation [237]. In a pioneering experiment at GSI, it was shown that the fission fragments 
produced by a ^^^U projectile beam impinging upon a lead target at 750 MeV/nucleon can 
be separated and identified [247] . 



6.1.2 Isotope Separation On- Line (ISOL) techniques 

In the ISOL method, radioactive nuclei are produced in a target, thermalized in a catcher 
(often target and catcher are one and the same system) , ionized in an ion source - if extracted 
from the target/catcher in atomic form - separated and reaccelerated [248]. In the traditional 
ISOL approach, a beam of accelerated stable nuclei bombards a target to produce short- 
lived nuclei via target fragmentation, spallation, direct reactions, fusion or fission. The so 
produced atoms are extracted from the target /catcher system - typically with heat applied 
to accelerate the diffusion process. The nuclei are transported from the target/catcher to 
the ion source. After extraction from the ion source, the low-energy ion beam is mass- 
separated and reaccelerated to projectile energies required by the experiment. Coolers and 
charge breeders are used to improve the ion optical properties of the extracted beam and 
to increase the reacceleration efficiency, respectively [248-250]. When gas catchers are used, 
a fraction of the thermalizing nuclei will still be in ionic form so that ionization in an ion 
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source can be avoided. This scheme is referred to as Ion Guide Isotope Separation Onhne 
(IGISOL) [251,252]. 

Various projectiles over a broad range of energies impinging upon a variety of targets are 
used to produce the rare isotopes of interest. Light-ion induced fusion evaporation or direct 
reactions produce a hmited number of neutron-deficient species and rare isotopes close to 
stability at high conversion rates. For example, at Louvain-la-Neuve, a 30 MeV proton beam 
on a thick ^'^C target (>1 g/cm^) is used to produce a high-intensity beam of ^^N with 
the ^^C(p, n)^^N reaction [253]. Heavy-ion fusion evaporation reactions allow to access very 
neutron-deficient nuclei. However, energy loss considerations require comparably thin targets 
(about 1 mg/cm^) and thus result in much lower conversion rates for these reactions. 

Projectile fragmentation reactions (e.g. heavy-ion beam on a thick carbon target [254]) and 
target fragmentation (e.g. high-energy protons on heavy targets [249,250]) produce a variety 
of nuclei in the regions close the initial projectile and target nucleus, respectively, as well as 
very light nuclei. For example, ^^Li is produced with this technique. 

The fission of 235,238^^ 232rp]^ long-lived actinides is used to produce neutron-rich nuclei 
over a wide mass range. Low- and high-energy protons, fast heavy ions, fast or thermal neu- 
trons and electrons (photo fission) are employed to induce fission. Most fission products have 
a large range in the target material and thus, unlike in heavy-ion induced fusion evaporation 
reactions, thick targets can be used to increase the yield of the wanted rare isotopes. In 
the collision of the high-energy proton beam with a target, spallation, fragmentation and 
fission can occur [255-258] . In spallation reactions, a large numbers of protons, neutrons and 
a particles are ablated from the target nucleus, for example ^^^U. Spallation products are 
typically neutron deficient while fission gives access to neutron- rich nuclei. 

Due to the available high flux of protons or light ions, ISOL can provide very intense low- 
energy beams of rare isotopes with very high beam quality. However, the effectiveness of 
the thick-target ISOL technique strongly depends upon the chemistry of the wanted ele- 
ment, with large differences observed for different elements [259]. Refractory elements (e.g., 
vanadium, zirconium and molybdenum) , that have a very high melting point and low vapor 
pressure, cannot be released from an ISOL target in atomic form. Chemical evaporation 
techniques and the formation of molecular sidebands are being discussed to overcome this 
limitation for some refractory elements [260,261]. As the thick-target ISOL technique relies 
on the diffusion of the wanted atoms out of the target and their effusion into the ion source, 
decay losses significantly hmit the intensity for species with short lifetimes [262-264]. For nu- 
clei with lifetimes of miUiseconds or less, decay losses are often the most limiting factor. The 
gas-catcher based IGISOL scheme is characterized by shorter release times and is applicable 
to refractory elements and other species difficult to ionize. 

ISOL techniques are for example used at HRIBF/ORNL (US) [265], ISAC/TRIUMF (Canada) 
[266], ISOLDE/CERN (Switzerland) [267], SPIRAL/GANIL (France) [268], and Louvain-la- 
Neuve (Belgium) [269]. The IGISOL scheme is used at Jyvaskyla (Finland) [252], Louvain- 
la-Neuve (Belgium) [270], ANL (US) [271] and was recently developed for rare isotopes 
produced by projectile fragmentation at NSCL/MSU (US) [272]. 
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Fig. 26. Schematic view of the production of radioactive nuclear beams with in-flight separation and 
the ISOL technique. Also shown is the hybrid approach envisioned for the facility for rare isotope 
beams (FRIB) in the US, employing radioactive beams produced by projectile fragmentation or 
fission, thermalized in a gas stopper, for example, and reaccelerated. 

6.1.3 Outlook 



The Facility for Rare-Isotope Beams (FRIB) [273] planned in the US will provide both, 
fast exotic beams produced in-flight via projectile fragmentation or fission and reaccelerated 
rare-isotope beams at various energies. In an IGISOL-based scheme, short-lived isotopes pro- 
duced in a chemistry-independent way by projectile fragmentation or projectile fission will 
be separated in a fragment separator, thermalized in a gas catcher and quickly extracted 
(order of milliseconds) with the possibility of reacceleration to the beam energies required 
for experiments. The gas-catcher based thermalization of fast exotic beams produced by 
projectile fragmentation is presently operational at NSCL/MSU and a reacceleration capa- 
bility is presently being implemented [274]. At ANL, the CAlifornium Rare Isotope Breeder 
Upgrade (CARIBU) project [275] will produce neutron-rich isotopes from a "^^"^Cf source for 
reacceleration through the ATLAS facility and provide new opportunities for experiments 
with reaccelerated Cf fission fragments. 

In Europe, design studies for the next-generation European ISOL facility EURISOL are on- 
going [276]. The aim is to increase the variety of exotic nuclei produced and to increase the 
production yields by orders of magnitude compared to present-generation ISOL facilities. 
EURISOL would complement the Facility of Antiproton and Ion Research (FAIR) [277] in 
Germany, where high-energy rare-isotope beams (GeV/nucleon) will be produced by frag- 
mentation and fission, separated in-flight with the possibility of passive deceleration or cool- 
ing in storage ring structures. The SPIRAL2 and ISAC-II(I) upgrades are underway at 
GANIL [278] and TRIUMF [279], respectively. A major upgrade of the RIKEN facility (RI 
Beam Factory - RIBF) in Japan became operational recently [280] . 



6.2 Experimental techniques and applications 



This section describes experimental approaches and detection systems developed or adapted 
to probe properties of exotic nuclei that are important for nuclear astrophysics. Many ob- 
servables ranging from nuclear masses and ground-state half-lives to specific reaction cross 
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sections are nuclear physics input crucial for the understanding of a large variety of astrophys- 
ical processes and scenarios. Wc tried to focus here on the measurements and developments 
that were not covered in the review by Smith and Rehm from 2001 [281]. 

6.2.1 Targets 

Hydrogen- and helium-induced reactions are of great importance in many astrophysical sce- 
narios. A variety of hydrogen and helium (enriched) target schemes exist. Thin plastic tar- 
gets, for example, polyethylene - (CH2)n - foils with thicknesses between tens of jig/cvo? 
and several mg/cm^, have been used for proton- induced reactions. As helium does not form 
compounds, solid targets containing helium must be prepared by implantation. Helium- 
implanted aluminum targets have been developed at Louvain-la-Neuve [282]. The carbon 
in plastic targets and the heavier elements serving as carriers in implanted targets increase 
the energy loss of the beam passing through the target and induce intense background re- 
actions, typically dominated by elastic scattering. Aside from physical radiation damage at 
high beam intensities, hydrogen depletion in plastic targets poses problems as well. 

Gas cells have been developed to mitigate some of these problems. While gas cells with 
windows arc easy to handle, even thin windows can degrade the beam energy and produce 
unwanted background reactions. A He-gas filled chamber is used at Louvain-la-Neuve [283], 
while cryogenic helium and hydrogen gas cells have been built for experiments at ANL [218] 
and GANIL [284], for example. When excitation functions are of interest, the so-called thick- 
target inverse kinematics (TTIK) technique [285] can be employed. For this approach, the 
beam is slowed down and stopped in a thick target with the light reaction residues escaping 
and being detected. With this, a reaction can be probed over a wide energy range with one 
incoming beam energy. The TTIK technique for He-induced reactions on low-energy ^"^O 
beams has been used at CRIB in Tokyo [286] and at Texas A&M University [287] with 
thick He gas targets; a thick polyethylene target was used at HRIBF at Oak Ridge National 
Laboratory to study resonant reactions with ^^F and ^^F beams [288,289]. 

Window-less, differentially-pumped gas targets eliminate background reactions and increased 
energy loss arising from the presence of windows but require several stages of pumping at 
high pumping speeds to drop the pressure by orders of magnitude over short distances. At 
present, window-less gas targets are used in conjunction with the DRAGON setup [290] and 
at HRIBF/ORNL [291,292]. 

6.2.2 Direct measurement of (p,^) and (a, 7) radiative capture reactions 

Charged-particle induced radiative capture reactions such as (p, 7) and (a, 7) occur in many 
stellar environments, for example, in novae and X-ray bursts. In explosive environments - due 
to the high temperatures and short reaction times - capture reactions involving short-lived 
nuclei play an important role for energy generation and nucleosynthesis. 

Resonant capture reactions on stable and long-lived radioactive targets have been inves- 
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tigated with intense proton and a beams [293,294]. The most recent examples for direct 
measurements of radiative capture reactions are ^^Na(p, 7)^^Mg and ^'^0{p^^Y^Y studied at 
the Laboratory for Experimental Nuclear Astrophysics (LENA) at TUNL [295-297] and the 
3He(a,7)^Be [298-300] and ^^N(p,7)^50 reactions [301-304] at the Laboratory Underground 
Nuclear Astrophysics (LUNA) facihty [305,306] located in the Gran Sasso Laboratory. 

For the resonant capture on short-lived nuclei, the experiments have to be performed in 
inverse kinematics. Recoil separators are typically used to detect and identify the recoil- 
ing reaction products and to reject the direct beam. The separators that are presently 
used for inverse-kinematics capture reactions with stable and rare-isotope beams include 
DRAGON [290] at TRIUMF (Ganada), ARES [307] at Louvain-la-Neuve (Belgium), the 
Daresbury Recoil Separator (DRS) [291] at HRIBF/ORNL (US), the FMA [308] at ANL 
(US) and, for stable beams only, ERNA [309] at Bochum (Germany). The STrong Gradient 
Electro-magnetic Online Recoil separator for capture Gamma ray Experiments (St. George) 
is under construction at the University of Notre Dame and will be used for inverse-kinematics 
(a, 7) capture reactions induced by a He-jet target [310]. 

The pioneering measurement with a rare-isotope beam was the first direct determination of 
the ^^N(p, 7)^'^0 reaction cross section using a radioactive ^^N beam produced at Louvain- 
la-Neuve [311,312]. A ^^N beam with an intensity of 3 x 10^ particles per second and an 
energy of 8.2 MeV impinged upon a 180 //g/cm^ polyethylene target. A surface-barrier diode 
located at 17° with respect to the beam axis detected the scattered ^^N and the carbon and 
hydrogen target recoils. Assuming pure Rutherford scattering on the ^^C in the target, the 
number of incident ^'^N was determined. The resonant radiative proton capture was tagged 
via the 5.173 MeV capture 7-ray originating from the ground-state decay of the first excited 
1~ state of ^^O measured with a large volume Ge diode. The peak-to-total in the 7-ray 
spectrum was improved by vetoing cosmic rays with Gherenkov detectors positioned above 
the 7-ray detector setup. The F^ width of the dominant ^"^0 1" level was determined for 
the first time from a direct measurement. Gompared to the previously adopted ^^N(p, 7)^^0 
reaction rate, the results of the measurement by Delbar et al. suggested that the reaction 
actually proceeds twice as fast [313,312]. After studies of the non-resonant capture cross 
section [314], the 7 width of the 5.173 MeV state in -"^^O now poses the largest uncertainty 
of this reaction rate [315]. 

The later study of the ^^Ne(p, 7)^°Na reaction at Louvain-la-Neuve used the resolving power 
of the recoil separator ARES to identify ^°Na events from capture through the 2.643 MeV 
level [316]. An upper limit of the resonance strength of the 2.643 MeV level was determined 
with this approach that bypasses low-efficiency 7-ray detection. It was concluded that an 
increased transmission and efficiency of ARES would be required to be able to determine 
more than an upper limit [316]. 

Some of the most recent resonant capture reactions induced by radioactive and stable beams 
in inverse kinematics have been measured at the DRAGON facility at TRIUMF. DRAGON 
is a recoil mass separator at ISAG dedicated for the measurement of low-energy reactions of 
astrophysical interest. Radioactive and stable beams with energies between 0.15 and 1.5 MeV 
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are delivered by ISAC. The DRAGON facility consists of a recirculating, differentially- 
pumped window-lcss gas target surrounded by a 30-crystal BGO 7-ray detection array, an 
electromagnetic separator and a heavy- ion recoil detection system [290]. At the DRAGON 
facility, seven resonances in ^^Mg have been characterized in the reaction ^^Na(p, 7)^^Mg 
at center-of-mass energies between Ecm — 200 — 1103 keV [317,318], the strength of the 
Ecm — 184 kcV resonance in the ^^^Al(p, 7)^^81 reaction was determined for the first time 
from a direct measurement in inverse kinematics [319] and the '^°Ca(a, 7)'^'^Ti reaction was 
probed with a ^'^Ca beam of 0.6 - 1.15 MeV/nucleon covering the relevant temperature range 
for the a-rich freeze-out during a core-collapse supernova [320]. 

At HRIBF/ORNL, the cross section of the ^'''F(p, 7)^^Ne reactions was measured directly 
for the first time [321]. A cocktail beam of ^''O and ^'''F was produced at HRIBF and im- 
pinged upon ORNL's differentially pumped, windowless hydrogen target. The Daresbury 
Recoil Separator was tuned to transmit the ^^Ne recoils to the focal plane where they were 
identified in the gas-filled ionization chamber. Two normalization schemes for counting the 
incoming beam were used, (i) the detection and counting of elastically scattered protons in 
surface-barrier detectors and (ii) beam-current measurements facilitated by a metal plate 
and two plastic scintillator paddles; the decays of ^^F nuclei implanted into the metal plate 
that was inserted into the beam from time to time and retracted to sit between the scintilla- 
tor paddles were measured [321]. The strength of the 3+ resonance in ^^Ne (about 600 keV 
above the proton separation energy corresponding to a 10.83 MeV ^^F beam energy for 
the resonant-capture measurement) was determined. At an off-resonance ^''F beam energy 
(800 kcV above the proton threshold), an upper limit for the direct capture away from the 
dominating resonances could be determined [321]. Compared to previous nova nucleosynthe- 
sis calculations [322], the abundance of the important galactic 7-emitter ^^F, the /?-dccay 
daughter of ^^Ne, increased in the hottest zone of a 1.35 solar mass white dwarf nova by a 
factor of 1.6 over the previously adopted rate. The uncertainty of the measured rate leads 
to a spread of a factor of 2.4 in the final abundance of ^^F, while varying the previous res- 
onant contribution by a factor of 10 resulted in a spread of as much as 15-16 times in the 
final abundance of ^^F [321]. Further measurements, for example the determination of the 
^^F(p, 7)^^Ne direct capture cross section and the ■'^^F(p, q;)^^0 reaction cross section, are 
needed to further reduce uncertainties in the production of ^^F in novae [321]. The impact 
on X-ray bursts is under investigation; there are indications that the abundances of ^^O and 
^^F in X-ray bursts can be modified by more than an order of magnitude at a reduction in 
uncertainty by a factor of 20 compared to previous network calculations [321]. 

The cross section for radiative capture has been observed to show the characteristic reso- 
nances of the reaction on top of a background that slowly varies with beam energy [25,293]. 
This smooth background is attributed to the direct capture process - a transition for the 
projectile from an initial continuum state to a final bound state via the electromagnetic in- 
teraction - that depends on the properties of the bound states of the compound nucleus [25] . 
This direct non-resonant capture can become important for stellar reaction rates when the 
nuclear level density is low in the region of the Gamow window. Also interference effects be- 
tween resonant and non-resonant capture can occur. For capture reactions involving stable 
nuclei important in hydrostatic burning, the contributions from direct capture have been 
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estimated from capture cross sections measured away from resonances [323] . The total cross 
section for direct proton capture is of order [293] in the relevant energy range and except 
for the case of ^Be(p, 7)*B [324,32], direct measurements of direct capture with radioactive 
beams have not been performed due to the high beam intensity requirements. 



6.2.3 Coulomb dissociation - the time-reversed approach to radiative capture 

The dissociation of a fast- moving projectile in the Coulomb field of a high-Z nucleus - 
Coulomb excitation of continuum states - was proposed as alternative method to determine 
radiative capture cross sections [126,325]. For the Coulomb dissociation approach, the resid- 
ual nucleus B of the radiative capture reaction of interest, A{x,'-f)B, is Coulomb excited in 
the electromagnetic field of a high-Z nucleus to an unbound state that decays into A + x. 
Coulomb dissociation is mediated by the absorption of a virtual photon, B{'jvirt,x)A, which 
is related to photodisintegration via the virtual photon method [326,327]. As photodisinte- 
gration is the time reverse process of radiative capture, the cross sections of the two processes 
can be related via the detailed balance theorem (see Section 5.2.2 for the theoretical formal- 
ism of Coulomb dissociation). 

In Coulomb dissociation measurements at intermediate beam energies, the relative energy 
of the residues A and x is determined from their invariant mass, which is deduced from the 
measured velocity vectors of A and x. The relative energy corresponds to the center-of-mass 
energy of the A{x,j)B capture. Experimental advantages of the Coulomb breakup method 
are threefold, (i) the typically large number of virtual photons leads to a big Coulomb dis- 
sociation cross section, (ii) the photoabsorption cross section is favored compared to the 
radiative capture cross section by two to three orders of magnitude (phase space), and (iii) 
thick targets can be used in the regime of fast beams (higher luminosity) [328]. A challenging 
requirement for experiments, however, is the need for the measurement of angular distribu- 
tions. The multipolarities have to be disentangled since the electromagnetic multipoles con- 
tribute with different strengths to Coulomb excitation and capture processes. Particularly 
important are the non-negligible contributions and interference of electric dipole (£^1) and 
electric quadrupole {E2) excitations. The applicability of the Coulomb dissociation approach 
is limited by the fact that preferentially El modes are excited and that the detailed balance 
theorem cannot be used for the extraction of the radiative capture cross section in case the 
capture feeds bound excited states that decay by 7-ray emission with unknown branching 
ratio [328]. The main body of experimental work utilizing Coulomb breakup of rare isotopes 
for nuclear astrophysics has been aimed at exploring radiative proton capture reactions. 

An extensive experimental program that utilizes Coulomb breakup reactions for the extrac- 
tion of radiative capture cross sections was initiated first at RIKEN/ Japan by Motobayashi 
et al. [329,330]. In the pioneering experiment, the electromagnetic transition strength for the 
excited 1~ resonance in ^"^0 was determined in the Coulomb dissociation, ^"^O ^^N + p, in- 
duced by 87.5 MeV/nucleon ^^O impinging upon a 350 mg/cm^ Pb target [329]. The Coulomb 
dissociation cross section was determined from coincidence spectroscopy of protons and ^^N 
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in position-sensitive detection systems consisting of telescopes of 24 position-sensitive Si de- 
tectors for AE measurements backed by CsI(Tl) scintillators for E measurements [331] and 
plastic-scintillator hodoscopes with x — y position sensitivity. The complete kinematics - 
total energy and relative momentum vectors of the p + ^^N system - were determined [329]. 
The radiative width F-y was deduced and found in agreement with the direct measurement 
performed at Louvain-la-Neuve [311] (see Section 6.2.2 for the direct measurement). 

Subsequent experiments at RIKEN were aimed at the Coulomb dissociation of [127,332,333] 
to extract the reaction rate of the crucial ^Be(]9, 7)^B reaction, which is of great importance 
for the neutrino production in the Sun through the (5 decay of ^B. In the first experiment in 
1994, the breakup of a ^B beam of 46.5 MeV/nucleon incident on a 50 mg/cm^ ^°^Pb target, 
used a two-layered, x — y position-sensitive plastic-scintillator hodoscope with an active area 
of 1 X 0.96 m^ for the detection of the coincident ^Be + p breakup products. The hodoscope 
with a AE layer of 5 mm thick plastic (10 segments) and the E plane of 60 mm thick plastic 
(16 segments) was positioned 5 m downstream of the Pb target. The energy of the breakup 
fragments was deduced from the time-of-fiight (TOP) over the 5 m flight path, the angle 
from the location of the hits in the segmented hodoscope. The particle identification was 
performed with the AE — E method and from TOF-E information. A helium-filled bag was 
inserted between the target and the hodoscope to reduce the background due to breakup in 
air [127]. In a second experiment [332,333], the hodoscope was modified for higher count-rate 
capability and moved closer to expand the angular range covered from 6° to 10°. Also, the 
Nal(Tl) scintillator 7-ray detection array DALI was positioned around the target to tag the 
fraction of the breakup leading to the 429 keV bound excited state in ''Be. The larger angu- 
lar range covered by the hodoscope allowed for angular-distribution measurements and an 
assessment of the E2 contribution to the breakup process [332,333]. Similarly, the reactions 
^B(p, 7)^0, ^■'^C(p, 7)^^N and ■'^^N(p, 7)^^0, important for the hot pp mode nuclear buring 
in hydrogen-rich, very massive objects [211], were studied via Coulomb breakup at RIKEN 
to extract the reaction rates relevant to explosive hydrogen burning [330]. Subsequent mea- 
surements on sd-shell nuclei aimed at the study of the breakup of ^^Al [334] and ^^P [335] 
into ^^Mg -|- p and ^^Si -|- p, respectively, continued to use 7-ray spectroscopy with DALI 
and employed a position-sensitive Si detector telescope in front of the plastic hodoscope for 
the identification of the projectile-like breakup residues. The experimental results in the sd 
shell are relevant to the reaction path in Ne novae [336], where at high temperature and 
density many nuclear reactions involving rare isotopes contribute in the hot CNO cycle and 
the NeNa- and MgAl-cycles [13]. A specific signature is the nucleosynthesis of long-hved 
galactic 7 emitters such as ^^Na and ^^Al, and nucleosynthesis up to the silicon and calcium 
range [337]. Assuming temperature and density conditions given by nova models [338], the 
radiative width obtained in the RIKEN study indicates that the main reaction path favors 
the j3 decay rather than the proton capture on ^^Mg [334]. For ^^Si + p, the preliminary re- 
sult of the 7-decay width of the first excited state in ^^P is ten times smaller than the value 
estimated based on a shell-model calculation in [339]. This indicates that the ^^Si(p,7)^^P 
reaction does not contribute significantly to the amount of Al produced in novae [335] . 

At GANIL, the Coulomb dissociation cross sections of ^^O + p [340] and ""^^N 

^^C -|- p [341] were measured. Radioactive beams of ^''O at 70 MeV/nucleon and ^^N at 
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65.5 McV/nucleon were provided by the ALPHA spectrometer and guided onto 100 mg/cm^ 
and 120 mg/cm^ ^°^Pb targets, respectively. In both experiments, the heavy breakup residues 
were characterized with the SPEG spectrometer and its focal-plane detection system con- 
sisting of two drift chambers and a plastic scintillator. The protons were detected in CsI(Tl) 
detectors downstream of the Pb target [340,341]. 

Much higher beam energies are available for Coulomb dissociation experiments at GSI. The 
breakup of ^Be -\- p was measured in two experiments at 254 MeV/nucleon ^B beam 

energy [38,40]. Both experiments utilized the KaoS spectrometer [342]. Beam particles were 
identified with the TOF-AE' method using a plastic timing detector upstream of the tar- 
get and a large-area scintillator wall close to the focal plane of KaoS. Momentum vectors 
of the reaction residues were analyzed by trajectory reconstruction with KaoS using posi- 
tion sensitive microstrip detectors and two 2d multi-wire proportional counters (MWPC). 
The first experiment used a 199 mg/cm^ 208p]^ target, the second experiment used a Pb 
target of 50 mg/cm^ thickness. Tracking of the angle of the incoming beam with parallel- 
plate avalanche counters allowed to measure angular distributions to disentangle contributing 
multipolarities in the second experiment [40,343]. 

At MSU, inclusive and exclusive measurements of the ^B —>■ ^Be -|- p Coulomb dissoci- 
ation were performed [39,344]. Inclusive measurements used 44 and 81 MeV/nucleon ^B 
beams provided by the A1200 fragment separator at NSCL. The S800 spectrograph [345] 
was used to detect the ^Be residues from the Coulomb dissociation of ^B impinging upon 
27 mg/cm^ Ag and 28 mg/cm^ Pb targets. The longitudinal momentum distributions of ^Be 
were derived from trajectory reconstruction using the position information provided by the 
position-sensitive cathode readout drift chambers of the spectrograph's focal-plane detec- 
tion system [346] . Particle identification was achieved with the AE — E method employing 
the S800 ionization chamber for AE and the first focal-plane plastic scintillator for E mea- 
surements [346]. In comparison to model calculations, a rather high E2 contribution to the 
Coulomb breakup was deduced, which has not been confirmed by any other measurement to 
date. An exclusive measurement at 83 MeV/nucleon on a 47 mg/cm^ Pb target in front of a 
1.5 T dipole magnet was performed in addition. Two multi-wire drift chambers (MWDCs) 
each were used to track the flight paths of the ^Be and the proton after passage through 
the magnetic field. An array of plastic scintillators was placed behind the MWDCs. Particle 
identification was performed with the energy loss measured in the plastic scintillator array 
and the TOP information. The invariant mass method was used to calculate the relative 
energy of the breakup residues [39,344]. 

In a recent comprehensive study, M. Gai [347] summarized the results of the various ^B 
Coulomb dissociation measurements and the extracted astrophysical 5*17(0) factor in com- 
parison to the results of direct capture measurements. 

Recently, the neutron capture reaction on ^^C has been explored using Coulomb dissociation 
of ^^C ^ ^^C + n at 69 MeV/nucleon in the field of a 224 mg/cm^ Pb target at RIKEN [348]. 
The ^^C residues were analyzed with a magnetic spectrometer that used a drift chamber and 
plastic scintillator hodoscopes. The momentum vector of ^"^C was determined by combining 
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the tracking data and time-of-flight information. Particle identification was acliieved witli 
Ai? - TOF and tracking information. The breakup neutron was detected by two layers of 
neutron hodoscope arrays with an active area of 2.14 x 0.92 m^ positioned almost 5 m 
downstream of the Pb target. The neutron momentum vector was extracted from the hit 
position of the neutron and its TOF. The spectral shape and amphtude was found consistent 
with similar data taken at higher beam energies at GSI [348]. The deduced direct capture 
cross section agrees with the most recent direct capture measurement and demonstrates 
that the Coulomb breakup with a neutron in the exit channel is a useful tool to derive 
the radiative capture cross section. A program aimed at Coulomb breakup of neutron-rich 
nuclei, not exclusively with an astrophysical background, is also ongoing at GSI at higher 
beam energies [349-351,132]. 



6.2.4 Direct ip,a) and {a,p) measurements 

Proton-induced reactions are of great importance in many astrophysical scenarios as the 
corresponding reaction cross sections are high due to the small Coulomb barrier. The pre- 
dominant reactions in nova outbursts, for example, are of type {p, 7) and {p, a). A significant 
fraction of the energy produced during nova explosions stems from proton-induced reactions 
on pre-existing carbon, nitrogen, and oxygen seed nuclei. The {p, a) reactions play a crucial 
role in cycling reaction flow back to lower masses. At sufficiently high temperatures, breakout 
sequences can occur and heavier nuclei are synthesized through reactions on neutron-deflcient 
nuclei. In X-ray bursts, for example, a sequence of {a,p) and (p, 7) reactions (ap process) 
leads to nuclei up to mass A = 40, eventually followed by rapid proton capture reactions (rp 
process) and (3 decays resulting in the formation of nuclei with A pa 100 [6]. 

The study of reactions with charged particles in the exit channel advanced tremendously with 
the advent of large-area silicon strip detectors. Si strip detectors can be single-sided - pro- 
viding segmentation in one dimension - or double-sided for two-dimensional segmentation. 
Detector thicknesses between 50-1000 /xm are typically available. A high degree of segmenta- 
tion also allows for higher total counting rates (without pileup) than could be achieved with 
a single detector of the same total area. When operated in transmission geometry, silicon 
detectors can be stacked to enable Z identification of a broad range of particles by measuring 
energy loss and total energy {AE — E method). 

The Louvain-Edinburgh Detector Array (LEDA) [352] was one of the pioneering charged- 
particle arrays for the detection of reaction residues in low-energy nuclear astrophysics exper- 
iments. For LEDA, the strips are curved in a circular pattern to realize an annular geometry. 
The CD-like detector consists of eight sectors with 16 strips (5 mm pitch) per sector. For such 
a configuration, the range of scattering angles and the angular resolution are determined by 
the target-detector distance. The removal of two sectors from the complete annular LEDA 
detector allows the remaining six sectors to be arranged as a six-sided cone (referred to as 
the LAMP configuration) providing very large solid-angle coverage at poorer angular reso- 
lution [352]. ORNL's Silicon Detector Array (SIDAR) [353] consists of single-sided silicon 
strip detectors of the LEDA- type (128 strips each) and a smaller CD-like detector with 64 
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strips in total for the detection and identification of the heavy recoil [354,353]. The TRIUMF 
UK Detector Array (TUDA) facility at ISAC/TRIUMF is a cylindrical scattering chamber, 
where large-area silicon detectors of the LEDA-type or similar CD-shaped detectors can be 
mounted [355]. 

For the TUDA facility at TRIUMF, a large solid angle active target for the detection of 
low-energy charged particles in measurements of astrophysical reaction rates is being devel- 
oped. The TRIUMF Annular Chamber for Tracking and Identification of Charged particles 
(TACTIC) [356] is a cylindrical ionization (time projection) chamber with segmented anode 
strips for the energy loss and total energy determination of light charged particles emerging 
from reactions. Measurements of the drift time provide trajectory reconstruction. Gas elec- 
tron multipliers (GEM) will amplify the signals. Digital readout and pulse shape analysis 
of the signals is envisioned [356]. The chamber gas will serve as target similar to the active 
target device MAYA [357]. 

The ^^F{p,ay^O reation has received much attention as it is of importance for nova 7-ray 
astronomy [358]. The 7-ray spectrum emitted from novae in the first few hours after expan- 
sion is dominated by the positron annihilation from the f3~^ decay of radioactive nuclei. One 
of the main contributors is ^^F with a /3-decay half-life of 110 min. Therefore, its production 
and destruction - via ^^F(p, q;)^^0 - is of great importance. Direct measurements with ^^F 
radioactive beams have been performed at LLN [359-362], ANL [363] and ORNL [364-366]. 
Most of the early measurements [359,361,363,364] where performed at center-of-mass ener- 
gies exceeding 550 keV, which dominate the reaction rate at temperatures above 0.4 GK. 
Graulich et al. in 1997 [360] found evidence for a 3/2~ resonance at Ecm = 330 keV but 
with insufficient statistics to determine the resonance strength. In 2002 Bardayan et al. [365] 
used a ^^F beam provided by the Holifield facility at ORNL with an intensity of 2 x 10^ 
particles/s interacting with a CH2 target to measure the ^^F(p, a)^^0 reaction. The 
coincidences were unambiguously identified with the SIDAR setup and the 3/2^ resonance 
strength at Ecm = 330 keV could be determined for the first time. This resonance domi- 
nates the reaction rate over a range of temperatures important for ONeMg novae [365] . The 
region of lower center-of-mass energies, which is relevant for typical nova temperatures (~ 
0.25 GK), has remained unaccessible to direct measurements. Recent work at LLN [362] and 
ORNL [366] constrains the interference of 3/2+ resonances located just above the proton 
separation energy from new measurements at Ecm > 400 keV and Ecm = 663 — 877 keV, 
respectively, and from the wealth of previous data in comparison to 7?.-matrix calculations. 

In explosive hydrogen burning up to a temperature of 0.2 GK, the burning with carbon takes 
place through a series of reactions known as the hot CNO cycle. At higher temperatures of 
about 0.4 GK, the ^''O waiting point can be bypassed by a chain of reactions initiated 
by ^^0(q;,p)^^F. At still higher temperatures, the breakout from the CNO cycle becomes 
possible. The reaction that dominates the leak rate is ^^0(q;, 7)-'^^Ne [367]. However, in the 
regime of high temperature, alternative breakout routes have been suggested, in particular 
the -^^Ne(a, j9)^^Na reaction [368]. While a direct measurement of the key rate ^^0(q;, 7)^^Ne 
has not been feasible so far, the -'^*Ne(Q;,p)^-'^Na [369,370] and ^^0(a,p)^^F [286] reactions have 
been probed directly in experiments induced by ^^Ne and ^^O radioactive beams at LLN and 
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CRIB- Tokyo, respectively. The experiments at Louvain-la-Neuve covered the energy region 
of = 2.04 - 3.01 MeV [369] and E,,„, = 1.7 - 2.9 MeV [370]. The measurements of the 
^®Ne(a,p)^^Na reaction rate were performed with a scattering chamber that consisted of a 
gas target and a vacuum side. The incoming ^^Ne beam interacted with the 500 mbar He 
gas target. Two double-sided sihcon strip detector telescopes in the gas volume for A£' — E 
measurements were used to detect the protons and to reconstruct their trajectories. The 
normalization of the incoming ^^Ne beam rate was based on the detection of clastically 
scattered ^^Ne off a gold foil in surface-barrier detectors on the vacuum side of the scattering 
chamber [283,369,370]. Groombridge et al. [370] identified eight states in the compound 
nucleus. Calculations of the enhanced stellar reaction rate using the new resonances as input 
show a good agreement with theoretical predictions by [368] . The experimental reaction rate 
represents a lower limit and causes a breakout from the CNO cycle via ^^Ne(Q;,p)^^Na to 
be delayed by several hundred ms when compared to calculations based on previous rates; 
further measurements are required at lower energies to map the lower-lying resonances above 
the a threshold [370]. 

At the CRIB facility, a ^''O beam interacted with a He gas target operated at a temperature 
of 30 K, enhancing the density compared to room temperature by a factor of 10 [286]. The 
^'^0(a,p)^^F reaction products were identified with the A£' — E methods in sihcon detector 
telescopes. The reaction cross section was measured with the thick-target technique over 
an energy range of Ecm = 0.8 — 3.8 MeV. This constituted the first direct measurement 
of the ^'^0(a,p)^^F reaction. The measured cross section differs from conclusions based on 
an indirect measurement and from a direct measurement of the time-reversed reactions as 
the contributions of the ^^F excited states could not be taken into account [286] . The result 
seems to suggest an increase of 50% for the ^^0(q;,p)^^F reaction rate which might impact 
the ignition phase of X-ray bursts [286]. 



6.2.5 Resonance properties from elastic and inelastic scattering experiments 

Some reaction rates of astrophysical interest are totally or partially dominated by the contri- 
bution of resonances. Important experimental methods to study the properties of resonance 
is resonant elastic and inelastic scattering. In particular, proton resonance scattering mea- 
surements have provided an extensive amount of data on unbound states in proton-rich 
nuclear systems relevant to reaction rates in explosive burning scenarios. 

In inverse-kinematics resonant scattering, a rare-isotope beam bombards a proton-rich tar- 
get - typically (CH2)n, rarely cryogenic H2 targets. The spectrum of the scattered protons 
depends sensitively on resonances present in the compound nucleus. When a resonant state 
is scanned in the appropriate energy range, the scattered-proton spectrum shows a distinct 
structure that allows to extract resonance properties - energy, width, spin and parity - in 
comparison to 7?--matrix theory. The energy of the scattered proton is related to the reso- 
nant energy and the shape of the spectrum allows to determine the width and the spin of 
the resonant state [371]. The protons are typically detected in position-sensitive Si AE — E 
telescopes or with annular Si strip detectors of the LEDA type. Important for the extraction 
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of resonance properties from the shape of the proton recoil spectra in comparison to theory 
is a complete understanding of the experimental resolutions that contribute to the shape of 
the detected proton spectrum [371]. 

Two different experimental approaches are used, the thick-target and the thin-target tech- 
nique. In the thin-target approach, see for example [372,353,373], the energy region of interest 
is scanned by using beams at different energy. At ORNL, resonances in the compound nucleus 
-^^Ne were studied with a ^^F beam at 15 different energies between 10 and 14 MeV [353] 
irradiating a 35 ixg/cw? CH2 target. The scattered protons were detected in SIDAR and ^^F 
passing through the target was identified in an ionization chamber [353] . In the thick-target 
approach, the beam is slowed down and stopped in the target and the clastic proton scatter- 
ing can be performed over a range of energies that depends on the choice of target thickness 
and beam energy. This technique relies on the fact that the protons escape from the target 
as their energy loss is negligible compared to that of the heavy ion. This technique was used 
for example by Angulo et al. at Louvain-la-Neuve [374]. A ^^Ne beam at 28 MeV impinged 
upon a 520 fig/cw? CH2 target. The recoil protons were detected at 20 different angles with 
LEDA and the first excited state of ^^Na could be characterized [374]. 

Inverse-kinematics resonant proton scattering measurements have been performed, for ex- 
ample, with the thick-target technique at the TUDA facility at TRIUMF/ISAC, ^^Na + 
p [375,376] and ^^Na + p [377]; with the thin and thick-target approach at HRIBF/ORNL, 
"^'''F-l-p [372,373] and ^^F + p [353,378]; at Louvain-la-Neuve with the thick-target technique, 
i3N,C + p [379], i^Ne + p [380], + p [361],^SNe + p [374], and ^Be + p [381]; at CRIB in 
Tokyo with the thick-target technique, "C + p [382], ^^Mg + p [383] ^^N+p [384,385], and 
^Be+p [386]; at the BEARS facihty in Berkeley, ^^0+p [387], ^^C+p [388] and ^^0+p [389]; as 
well as at Texas A&M [390,388], Notre Dame [391,392], ANL [393] and Spiral/GANIL [394]. 

Inelastic resonance proton scattering has been employed to probe properties of states in 
the compound nucleus in the cases where particle decay to excited states is energetically 
possible. Recently, proton emitting states in -"^^Ne were studied with the inelastic scattering 
reaction H(^^Ne,p)^^Ne*(p)^^F at Louvain-la-Neuve. Resonance energies and widths were 
assigned from the shape of the scattered proton spectrum, while spins were assigned using 
proton-proton angular correlations between recoil and decay protons [395]. 

Inelastic scattering gains importance when rates are derived from measurements of the in- 
verse reaction. For example, when the reaction A{p, a)B is measured, A is typically in its 
ground state and the rate of B{a,p)Ags can be deduced from the detailed balance. How- 
ever, contributions from B{a,p)A* leading to an excited state in A are not accessible from 
the inverse reaction with A in the ground state. Inelastic scattering A + p ^ A* can then 
reveal important excitations that may contribute. One example is the ^^0(q;,p)^^F reaction, 
which was studied from measurements of the inverse reaction ^^F(p, q;)^^0 [396,397]. Proton 
inelastic scattering ^'^F + p were measured at ANL [393] and ORNL [373] to characterize res- 
onances of importance for the ^^0(q;,p)^^F reaction rate. For the measurement of inelastic 
scattering, thin targets are typically used to separate elastic and inelastic channels from the 
spectroscopy of the protons [373] or 7-ray spectroscopy is used to tag bound final states, see 
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for example [386], where Nal detectors were used to measure ^Be(p,p')^Be*. 

6.2.6 Nucleon transfer reactions, nucleon knockout and population of excited states in frag- 
mentation and projectile fission 

For stable nuclei, light-ion induced transfer reactions in normal kinematics have been widely 
used to explore stellar reaction rates by extracting spectroscopic information on resonances 
close to the threshold that dominate the reactions of interest. The parameters of a resonance 
- energy, orbital angular momentum, partial and total width, spectroscopic factors and decay 
modes, for example - can be determined. In general, the sensitivity of transfer reactions to the 
single-particle degree of freedom continues to provide important data to benchmark nuclear 
structure models, in particular the nuclear shell model which provides important input to 
the modeling of many astrophysical processes (see Section 5.3 for details on the theoretical 
treatment of transfer reactions) . 

Neutron capture cross sections are important in the r process, in which heavier nuclei are 
formed from seed elements by consecutive neutron capture reactions and (5 decays. In an 
environment of high neutron density, tens of neutron captures may occur until the /3 decay 
half-lives become shorter than the half-life against neutron capture. Consequently, the process 
proceeds off the valley of stability towards neutron-rich nuclei. For most neutron capture cross 
sections, statistical models can be applied. However, near closed shells - in a regime of low 
level density - direct capture becomes important [398]. For direct capture calculations, level 
properties (excitation energies, spins and parities) have to be known accurately unlike for 
statistical calculations where averages over resonances are considered. 

Following pioneering inverse-kinematics {d^p) one-neutron transfer experiments induced by 
stable Xe beams impinging upon deuterated Ti targets performed at GSI [399] , a program of 
low-energy (d, p) transfer experiments in mass regions relevant to the r process was started 
at HRIBF at ORNL. The experimental study of (i(^^^Sn,p) in inverse kinematics at energies 
close to the Coulomb barrier proved that Q-value spectra and angular distributions can be 
extracted for low-energy inverse-kinematics transfer reactions with heavy beams [400]. This 
experiment used a deuterated polyethylene target (CD2) of 100 //g/cm^ thickness that was 
angled to achieve an effective thickness of 200 fig/cm^ and to enable proton detection under 
6iab = 90°. Two silicon telescopes were positioned covering angles of 6iah = 70—102° and 6iab = 
85 — 110°, respectively. The silicon detector array SIDAR was mounted in half-lampshade 
configuration covering backward angles of 9iab = 130 — 160°. The states were determined 
from the Q- value spectrum and the angular distributions were used to determine the orbital 
angular momentum (£ value) of the transferred neutron. An absolute normalization of the 
cross sections was performed using the elastically scattered dcuterons [400]. The results were 
found in agreement with normal-kinematics ^^^Sn((i,p) and demonstrated that this important 
mass region for the r process can be accessed with low-energy, inverse-kinematics transfer 
reactions. 

The first {d,p) transfer on a neutron-rich r process nucleus, ^H(^^Ge,p)^^Ge at 330 MeV on 
a 430 /ig/cm^ (CD2)n target, was subsequently performed at HRIBF [401]. Protons were 
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detected with SIDAR, subtending 6iah = 105 — 150°. A segmented ion chamber was em- 
ployed to determine the Z of the projectile-hke reaction residues. The measured reaction 
Q-value provided an indirect measurement of the mass of ^^Ge. The extracted neutron sep- 
aration energy Sn ~ 3.69 MeV was found low suggesting that the ^^Ge(n, 7)^^Ge reaction 
rate has a significant direct neutron capture component. The excitation energy of the 1/2+ 
excited state followed the falling trend as semi-magic ^^Ni is approached [401]. ^^Se was 
the main contaminant in the beam and the {d^p) neutron transfer to ^^Se was studied at 
the same time [402]. Direct-semidirect contributions to the neutron capture cross sections 
were computed [402]. Recently, neutron capture extracted from transfer reactions was 
benchmarked for ^^Ca(d,p)^^Ca at different deuteron energies [403]. 

One of the key reactions in nova 7-ray astronomy, ^^F(p, q;)^^0, has been explored indi- 
rectly by exploiting ^H(^^F,p)^^F neutron transfer reactions and the ^^Ne-^^F mirror sym- 
metry [404-406]. Spectroscopic information on excited states just above the proton separation 
energy in the compound nucleus ^^Ne is needed for the determination of the ^®F(p, q;)^^0 
reaction rate at nova temperatures. Some of the levels in ^^F + p cannot be reached in reso- 
nance scattering experiments since the resonances are well below the Coulomb barrier. With 
the availabihty of a ^^F rare-isotope beam, neutron spectroscopic factors of the analog levels 
in the mirror nucleus of ^^Ne, ^^F, can be studied with (rf, p) reactions in the corresponding 
energy region and - invoking mirror symmetry - serve as input to calculate the proton width 
in ^^Ne. At Louvain-la-Neuve, a 14 MeV ^^F beam impinging upon a 100 ng/cvo? CD2 target 
was used to induce the {d,p) neutron transfer reaction [404]. The relevant levels are above 
the a separation energy leading to ^^F* — > ^^N -|- a. The measurement employed the LEDA 
detector for the detection of ^^N (downstream) and the LAMP detector (backward) for pro- 
ton detection. Spectroscopic factors deduced in comparison to DWBA calculations were used 
to put new limits to the contribution of low-energy resonances to the ^®F(p, q;)^^0 reaction 
rate [404]. A similar experiment was performed at HRIBF at ORNL at a much higher en- 
ergy with a 108.49 MeV beam of ^^F on a 160(10) //g/cm^ CD2 target [405,406]. The protons 
emerging from the ^H(^^F,p)^^F transfer reaction were detected with the silicon strip de- 
tector array SIDAR at Oiab = 118 — 157°. Another silicon strip detector at the focal plane 
of the Daresbury Recoil Separator was used to tag particle-stable A = 19 reaction residues 
and a-decaying final states in coincidence with the protons detected in SIDAR. Neutron 
spectroscopic factors for eight [405] and 13 [406] analog levels of astrophysical relevance in 
the mirror ^^Ne were determined. The results implied significantly reduced ^^F(p, 7)^^Ne and 
^^F(p, q;)^^0 rates compared to what was reported previously [405]. 

The {d, n) transfer reaction has been used at the RESOLUT facility at Florida State Univer- 
sity as surrogate for proton capture. The reaction ^^Al((i, n)^^Si* — > p + ^^Al was induced by 
a 91.5 MeV beam of ^^Al impinging upon a 1.66 mg/cm^ thick CD2 target [407]. The decay 
protons were detected under forward angles with the AE — E method using a hybrid ion 
chamber and a double-sided silicon strip detector backed by a second, unsegmented silicon 
detector. The lowest £ = proton resonance in ^®Si could be identified. This shifts the main 
experimental uncertainty in the synthesis of ^^Al to the destruction of ^e™^} through the 
(p, 7) reaction, which is planned to be studied at RESOLUT in the future [407]. 
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The use of heavy-ion induced, peripheral, proton-adding transfer reactions of type (-""^Nj^^C) 
[408,409,314], (iOB,9Be) [410,411,408] and {^^C,^^C) [412] in inverse kinematics as well as 
(^He, d) in normal kinematics [413-415] has been pioneered at Texas A & M to indirectly 
assess the non-resonant direct-capture contribution in (p, 7) reactions. Asymptotic normal- 
ization coefficients (ANCs) [42,413,408] are deduced from the measured angular distributions 
in transfer reactions at energies above the Coulomb barrier where the cross sections are orders 
of magnitude higher than for direct capture processes at astrophysical energies (see Section 
5.3.2). Radioactive beams of ^Be (84 MeV), ^^C (10 MeV/u), and (11.8 MeV/u) were 
produced at Texas A & M and separated with the MARS spectrometer. ^^B(^Be, ®B)^Be, 
i%(^Be,«B)i3C, i4N(13n,i^O)13C, and ^'^NC^C,^^Ny^C transfer reactions induced by i°B and 
melamine {C^NqHq) targets, respectively, were used to determine the ANCs for ^B ''Be+p, 
^"'O -"^"^N + p, and ^^N ^^C + p to deduce the non-resonant capture rates for the 
^Be(p, 7)'^B, ^'^N(p, 7)^'^0 and ^^C(p, 7)^^N reactions. In these experiments, the elastically 
scattered projectile beam and transfer products are observed simultaneously in two detector 
telescopes consisting of position-sensitive 16-strip Si detectors backed by Si E detec- 
tors. Particle identification is then done via the Ai? — E method. The reaction telescopes 
were cooled down to -6°C to reduce thermal noise. A slightly different approach was taken 
at HRIBF at ORNL to study the ^^F(p,7)iSNe direct capture cross section [416]. A beam 
of 170 MeV ^^F impinged upon a melamine target inducing the ^^N(^^F,^^Ne)^^C reaction. 
Recoiling ^^Ne nuclei were detected via energy loss in position-sensitive Si strip detectors 
covering center-of-mass angles from 2 to 9°. To be able to distinguish the transfer to indi- 
vidual excited states, 7 rays emitted by ^'^Ne were detected in CLARION in coincidence and 
allowed to tag the final state. Preliminary results are reported [416]. 

The ANC technique was also apphed to a sub-Coulomb a-transfer reaction ^^C(^lA,d) per- 
formed at FSU with a ^^C beam at 60 MeV to study the astrophysical reaction rate of 

^^C(Q;,n)^^0 which - at low energies - is dominated by a 1/2"*" sub-threshold resonance in 
^^O [417]. The application of the ANC method also to one-nucleon removal and breakup 
reactions induced by weakly-bound projectile nuclei and light as well as heavy target nuclei 
was proposed [418,419]. 

A different experimental approach that utilizes the high luminosity of fast-beam experiments 
at fragmentation facilities was developed at NSCL. The one-neutron removal from a projec- 
tile beam of ^^Ar upon collision with a CH2 target at 84 MeV/nucleon was used to populate 
excited states in ^^Ar [168]. In-beam 7-ray spectroscopy using the segmented germanium 
array SeGA [420] provided precision information on the excitation energy of states in the 
proximity of the proton separation energy and for states above Sp that have a 7 branch. High- 
resolution 7-ray spectroscopy allowed to determine the excitation energy with uncertainties 
of several keV. The 2-orders of magnitude improvement in the uncertainty of the level energy 
translated into an 3-orders of magnitude improvement for the rate of ^^Cl(p, 7)^^Ar [168] (see 
Figure 27) . Similar measurements have been performed to precisely determine relevant level 
energies m 30s [421]. Pioneered at NSCL, ^Be-induced one-nucleon removal reactions from 
fast exotic beams have proven to be a crucial tool to study the detailed evolution of nuclear 
shell structure in exotic nuclei [158,422]. Their reach to the most exotic nuclei and the sen- 
sitivity to the single-particle degree of freedom continue to probe the predictability achieved 
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Fig. 27. Total astrophysical ^^Cl(p, 7)^^Ar ground-state capture rate as function of Tg [169]. The 
black region shows the capture rate restricted by the one-neutron knockout measurement of Clement 
et al. [168] relative to the previously possible range of rates (gray band). The total estimated 
uncertainty of the rate ranges from a factor of 3 to a factor of 6. (Figure adapted from Hendrik 
Schatz and Ref. [169]). 

by large-scale shell model calculations, which provide crucial input for many astrophysical 
processes (see also Section 5.4). 

A less selective method is the population of isomeric excited states in projectile fragmentation 
or fission and the spectroscopy of the level structure fed by the long-lived state(s). At GSI, 
7-ray decays of excited states in the = 82 r-process waiting-point nucleus ^'^'^Cd were 
observed unambiguously for the first time [423]. An 8^ isomeric state in ^^"Cd with a two- 
quasiparticle structure was populated in the fragmentation of a^^^Xe beam and in the fission 
of ^^'^U. In the two parts of the experiment, heavy, neutron-rich beams were produced by 
fragmentation of a ^^^Xe beam at 750 MeV/nucleon impinging upon a 4 g/cm^ ^Be target 
and by a 650 MeV/nucleon ^^^U beam fissioning upon collision with a 1 g/cm^ thick ^Be 
target. The produced ^^°Cd nuclei were separated with the GSI fragment separator (FRS) and 
identified via the measured energy loss, time of flight, magnetic rigidity, and various position 
measurements. The nuclei were implanted into a passive stopper surrounded by 15 large- 
volume Ge cluster detectors [424]. The new 2f energy of 1395 keV for ^^°Cd is not consistent 
with the previous tentative assignment of the 2+ state at 957 keV [425]. The evidence for the 
existence of an 8+ isomer in ^^°Cd is in line with expectations based on the valence analog 
^®Cd. These new results on the level scheme of ^^'^Cd provided no evidence for a reduction 
of the N = 82 shell gap in the vicinity of ^^^Sn contrary to what was implied by the results 
of a previous experiment [426]. In a similar measurement at GSI, excited states in the magic 
nucleus ^°^Pt (neutron number = 126) were populated in the fragmentation of ^°^Pb at 
1 GeV/nucleon beam energy [427]. Medium spin isomers with half-lives between 150 ns and 
55 IIS and the level structures fed by their decays were identified. The data suggested a 
revision of the two-body interactions for A^ = 126 and Z < 82, which is important for the 
evolution of nuclear structure toward the waiting points of the r process in this region [427]. 
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6.2.7 Weak-interaction strength 



Supernovae are a major source of nucleosynthesis. For both types of supernovae - core- 
collapse and thermonuclear - weak interaction rates play are crucial role [181,182,18]. In 
particular electron capture rates for many nuclei in the /p-shell and beyond are among 
the important nuclear physics ingredients needed for supernova models [428]. The weak- 
interaction rates are largely provided by nuclear structure calculations, as for example by 
large-scale shell-model calculations [183,185]. Important experimental benchmarks necessary 
for the reliable modeling of weak-interaction rates by nuclear theory are provided by mea- 
sured Gamow- Teller (GT) strength distributions (see Section 5.6 for details on the theoretical 
formalism) . 

GT transition strength [B{GT) values) can be obtained from Gamow- Teller (3 decays within 
excitation-energy limitations given by decay Q-value window. Charge-exchange reactions, 
however, can map the GT strength distributions over a wider range of excitation energies. A 
variety of charge- changing reactions can be employed, for example, '^Z(p, n)'^(Z+l) [429,173], 
^Z(n,p)^(Z-l) [430],^Z(rf,2He)^(Z-l) [431], ^Z(3He, t)^(Z+l) [432-435], and ^Z(t, 3He)^(Z-l) 
[436-439]. In charge-exchange reactions at beam energies exceeding 100 MeV/nucleon, the 
cross section at low momentum transfer (small angles) is proportional to the B{GT) tran- 
sition strength [440-442,173,174]. GT transitions have been probed extensively with the 
charge- exchange reactions induced by stable projectiles on stable targets (see references 
above). Charge-exchange reactions on short-lived exotic nuclei in inverse kinematics are 
proposed to be developed into spectroscopic tools at FAIR (Germany) [443] and at NSCL 
(US) [446] . At FAIR, the utilization of charge-exchange reactions is proposed by the EXL [447] 
and by the R^B [448] collaborations. At R^B, (p, n) charge-exchange reactions will be induced 
by rare-isotope beams at relativistic energies impinging upon a liquid hydrogen target with 
the emerging slow neutrons measured by plastic scintillators surrounding the target [448]. 
The EXL collaboration envisions the study of charge exchange via (p, n), (^He, t), and (ci, ^He) 
reactions with unmatched luminosity induced by rare isotopes circulating in a storage ring 
passing through a gas jet target [447] . Two avenues are taken at NSCL, the first program 
utilizes the '^Z(^Li, '^Bc)'^(Z-l) induced charge-exchange reaction where the spin transfer is 
tagged with the 430 kcV 7-ray decay (1/2^ 3/2^) in the target-like reaction residue 
^Be [444]. For a second program, the Low- Energy Neutron Detector Array (LENDA) [445] 
is presently constructed for slow-neutron detection in (p, n) charge-exchange reactions in- 
duced by rare-isotope beams bombarding a hydrogen target. In all NSCL experiments, the 
projectile-like fragments will be tracked and identified with the S800 spectrograph [446] . 



6.2.8 Beta-decay half-lives and (5 -delayed neutron emission probabilities 

Isotopes of elements beyond iron are almost exclusively produced in neutron-capture pro- 
cesses. The two main neutron-capture processes for astrophysical nucleosynthesis arc the 
slow s process and the rapid r process (see [4] for a review). The attributes slow and rapid 
refer to the timescale between subsequent neutron captures relative to the competing time 
scale for /5 decay. While the s process stays close to the valley of stability, the r process leads 
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to very neutron-rich nuclei. To determine the path of the r process, among others, /3-decay 
properties are crucial. Beta decay competes with the neutron capture and the drives the path 
toward stable nuclei and fission determines the heaviest nuclei produced. On the proton-rich 
side of the nuclear chart, /3-decay half-lives are an important input for nucleosynthesis mod- 
els that describe the reaction flow in the rp process and the flnal abundances. While much 
data has been available on the proton- rich side of the nuclear chart, some information in key 
regions, as for example around the nuclei ^^'^^Mo and ^^'^^Ru, which are found in the solar 
system in unexpectedly large abundances, could only be gathered recently with technical 
improvements crucial for experiments on proton-rich nuclei. 

Two approaches have been taken to measure the /9-decay half-lives and /9-delayed neu- 
tron emission probabilities. At the Jyvaskyla (IGISOL) and the CERN/ISOLDE facilities, 
neutron-rich isotopes of interest are produced by fission, mass-separated and implanted into 
a collection tape, which moves at pre-set times to suppress the /? and 7 activity from longer- 
lived isobaric contaminants. Time spectra are recorded for 7-ray, neutron and (3 events 
where the half-lives and production rates can be derived from the growth and decay curves 
of the collected activities during pulsed beam mode. At fragmentation facilities, tracking and 
event-by-event particle identification become possible in the regime of fast beams [449]. The 
fast- moving beam particles are, for example, implanted into a double-sided silicon strip de- 
tector (DSSD) with subsequent (3 decays detected on an event-by-event basis and correlated 
with the position of a specific implantation. All events are typically time-stamped and the 
decay curve is constructed from the time difference between implant and correlated decay. 
This method is applicable at rates of less than 0.1 particles per second. 

In early experiments at Jyvaskyla, the /?-decay half-lives and /3-delayed neutron emission 
branching ratios of neutron-rich Y, Nb, and Tc isotopes produced in the proton-induced 
fission of 238^ were determined [450,451]. The mass-separated nuclei of interest were im- 
planted onto a collection tape and their decays were measured inside the Mainz 47r neutron 
long counter used for the detection of /3-delayed neutrons. The /3-particle detection and 7-ray 
detection was performed with a thin plastic scintillator and a planar germanium detector, 
respectively, positioned inside the Mainz counter. The neutron long counter consisted of 42 
^He ionization chamber tubes, arranged in two concentric rings embedded in a polyethy- 
lene matrix surrounding the implantation point [450,451]. At CERN/ISOLDE, the /?-decay 
properties of i30,i3i,i32qj ^^^^ studied [452,426]. In 2000, the heavier i31'132qj isotopes were 
produced from proton-induced uranium fission and purified with resonance-ionization laser 
ion source (RILIS) and subsequent separation of the extracted nuclei with the General Pur- 
pose Separator (GPS). Beta-decay half-lives and /^-delayed neutron decay branching ratios 
Pn were determined using the j3n spectroscopy with the Mainz neutron long counter and 
a (3 plastic scintillator [452]. Superior separation was achieved for the study of ^^°Cd 
and its decay daughter ^^°In, in the vicinity of the Z = 50, = 82 r-process waiting 
point. ^^°Cd was produced by fission induced with fast reaction neutrons - generated from 
a 1 GeV proton beam impinging upon a Ta or W rod - interacting with the uranium car- 
bide/graphite production target. Highest chemical selectivity was achieved with RILIS and 
separation using the High- Resolution Separator (HRS). P'y and P'-ff spectroscopy around 
the moving-tape-collector was performed with AE — E (3 telescope and four large-volume 
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high-purity germanium detectors in close geometry [426] . 

On the proton-rich side of the nuclear chart and important for ONc novae, ^^Al was studied 
at Texas A&M [453]. ^^Al was produced via the p(^'^Mg, ^'^Al)2n reaction, separated by the 
MARS separator and collected with a moving-tape system to a detector setup for (3 and 
(3^ spectroscopy consisting of a thick plastic scintillator for [3 detection and a high-purity 
germanium detector for 7-ray detection [453]. Beta branching ratios and log ft values for 
transitions to levels in ^^Mg were extracted. From this, the ground-state spin and parity 
5/2+ of ^'^Al was determined unambiguously. This excludes the large increases in the radiative 
proton capture cross section for the reaction ^^Mg(p, 7) at astrophysical energies, which were 
implied by earlier claims that the spin and parity of the ^^Al ground state are 1/2+ [453]. 

Rare-isotope beams produced by projectile fragmentation are used at NSCL to measure the 
/3-decay properties of exotic nuclei of importance to nuclear astrophysics. In-fiight tracking 
and event-by-event particle identification become possible in the regime of fast beams. The 
fast-moving beam particles are implanted into a double-sided silicon strip detector (DSSD) 
with 40 horizontal and vertical strips resulting in a pixelation of 1600, the center piece of 
NSCL's Beta Counting Station (BCS) [454]. Subsequent (3 decays are detected on an event- 
by-event basis and correlated with the position of a specific implantation. Downstream of 
the DSSD, six single-sided Si strip detectors (SSSDs) and two Si PIN detectors comprise the 
13 calorimeter part of the BCS. All events in the BCS carry an absolute time stamp and the 
time difference between implant and decay serves to map the /3 decay curve. Coincident 7- 
ray spectroscopy and /3-delayed neutron spectroscopy can be performed with the Segmented 
Germanium Array (SeGA) [420] and the Neutron Emission Ratio Observer (NERO) [455] 
- a neutron long counter with ^He and BF3 gas counters. As the 7-rays are emitted by 
nuclei at rest, the segmentation of SeGA is not used. NERO consists of 60 proportional 
gas-counter tubes - 16 filled with ^He and 44 filled with B3F - embedded in a polyethylene 
moderator. The counters are arranged parallel to the beam direction in three concentric 
rings around the vacuum beam line that accommodates the BCS. A DSSD-based /3 counting 
system similar to the BCS is being developed for experiments at GSI with the RISING 
setup [456]. A challenge in implant and decay setups is the large difference in the energy 
response to fragments that are stopped and implanted in the DSSD (several GeV) and the 
energy signals of charged-particle decays (hundreds of keV). The /3-decay detection system 
for RISING will use logarithmic preamplifiers to cover the required dynamic range [456] 
while NSCL's BCS is using preamplifiers with dual gain capability, low gain for the fragment 
implantation and high gain for the detection of decays [454] . 

The half-life of the doubly-magic nucleus ^^Ni is important input for r-process model predic- 
tions of the nucleosynthesis around A = 80 [457]. At NSCL, ''^Ni was produced via projectile 
fragmentation from a beam of ^^Kr at 140 MeV/nucleon impinging upon a beryllium target 
of 376 mg/cm^ thickness. Each nucleus in the secondary beam was individually identified in 
fiight by measuring energy loss and time of flight (TOF). The ions were continuously im- 
planted into the DSSD. In the regime of very low beam rate for this exotic nucleus, the typical 
total implantation rate for the entire detector was under 0.1 per second. The DSSD registered 
the time and position of the decays following the implantation. This allowed the correlation 
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of 7 implanted nuclei with subsequent decay events. The measured half-life indicated a rather 
short time scale for the buildup of heavy elements beyond = 50 compared to some earlier 
predictions and thus unraveled an acceleration of the r process in this region [457]. 

For a similar experiment near the N — 82 neutron shell closure, neutron-rich Tc, Ru, Rh, Pd, 
and Ag isotopes were produced at the NSCL by fragmentation of a 121.8 MeV/nucleon ^^^Xe 
primary beam impinging onto a 206 mg/cm^ Be target [458]. The secondary beam was then 
transmitted to the experimental vault, where a second plastic scintillator provided the stop 
signal for the TOF measurement (relative to a timing signal from a scintillator located at the 
intermediate image of the A1900 fragment separator) before the beam was implanted into the 
NSCL's BCS for the study of the subsequent /3 decay. The average implantation rate in the 
DSSD was 0.4 Hz. Implantation and decay events were time stamped and correlated via their 
pixel locations. Coincident /3-delayed neutrons were detected with NERO. /9-delayed neutron 
emission branchings (or upper limits) for the neutron-rich nuclei ^^^~^^°Rh, i20-i22pj^ 
^^^Ag were determined as well as /3-decay half-hves for neutron rich ^^^~^^^Tc, ^^^~^^^Ru, 
116-121 j^j^^ and ii9-i24p^ j^gg^j. ^Yie proposed the r-process path [458]. A little lower in Z, the 
/5-decay properties of Y, Zr and Mo nuclei around A = 110 were measured [459]. Fast beams 
of neutron-rich Y, Zr, Nb, Mo, and Tc isotopes were produced at NSCL by fragmentation of 
a 120 MeV/nucleon ^^^Xe beam on a thick Be target. The species of interest were separated 
and identified with the A1900 fragment separator using the Bp — AE — Bp method and 
implanted in the BCS which was surrounded by NERO. New half-lives for -"^^^Y, io^.iot^j;^ 
and ^^^Mo, along with new P„ values for ^O'^Y and ^^'^'^^^Mo and P„ upper limits for ^^-^-^^'^'Zy 
and ^^'^'^^^Mo were determined. Analysis of the measured T1/2 and P„ values in the framework 
of QRPA calculations brought new insights in terms of deformation and shape coexistence, 
compatible with the hypothesis of a quenched N — 82 shell gap invoked by some to explain 
the r-process abundances around A — 110. An experimental complication at high Z is 
the presence of charge states in the secondary beam [459]. Measurements of the ion's total 
kinetic energy with the PIN detectors and the DSSD of the BCS are used to separate the 
fully stripped ions from charge states; isomer tagging with germanium detectors of the SeGA 
array is typically performed to confirm the particle identification. 

The = Z nuclei ^''Mo and ^°°Sn and nuclei in the vicinity of = Z = 50, ^^Cd and ^^In, 
that lie along the rp-process path became accessible for /?-dccay measurements at NSCL 
only after implementation of the Radio- Frequency Fragment Separator (RFFS). With the 
energies available at NSCL, secondary beams optimized on neutron-deficient species are of 
low purity due to extended low-momentum tails of higher-rigidity fragments which overlap 
with the fragment momentum distribution of interest. The RFFS applies a transverse RF 
electric field which defiects nuclear species in the secondary beam based on their phase 
difference with the primary beam. Secondary beams containing ^^Cd, ^®In, ^''"Sn [460] and 
^^Mo [461] were produced by projectile fragmentation of ^^^Sn and ^^^Xe primary beams, 
respectively, and purified by the RFFS. The BCS was used in both cases to measure the 
/3-decay half-lives. The new result for ^^Mo resolved a previously reported deviation from 
theoretical predictions [461] and the measurements around ^^Cd revealed that the rp process 
in X-ray bursts is not the main source of the unexpectedly large amount of ^®Ru in the solar 
system [460]. 
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More technical details and applications to nuclear structure research can be found in the 
recent review articles by Mantica [449] and Rubio and Gelletly [462]. 



6.2.9 Mass measurements 

Atomic masses of rare isotopes across the nuclear chart are among the key input data for 
large-scale reaction network calculations that quantify the nucleosynthesis, for example, via 
the rp process, up process and the r process (see for example [463]). Experimental methods 
for the determination of atomic masses basically fall into two categories. Approaches that 
measure the Q values in decays or reactions make use of Einstein's mass-energy equiva- 
lence; mass measurements that are based on the deflection of ions in electromagnetic fields 
determine the mass-to-charge ratio. The highest precision in mass spectrometry today is 
obtained by frequency measurements. The revolution or cyclotron frequencies of ions in a 
magnetic field are measured to deduce the mass-to-charge ratio in Penning trap experiments, 
in time-of-fiight measurements involving cyclotrons and in storage rings. The different ex- 
perimental methods are optimized for different beam energy regimes and are thus linked to 
the production and separation schemes of the rare isotopes. 

Penning traps use the three-dimensional confinement of ions with static magnetic and elec- 
tric fields. Mass measurements in a trap are based on the determination of the cyclotron 
frequency of the stored ions in a magnetic field. The ion motion in the combined magnetic 
and electric fields is the superposition of three (ideally) independent harmonic motions with 
eigenfrequencies lJz (axial motion), a;_|_ (modified cyclotron motion) and a;_ (magneton mo- 
tion). Their frequencies are: 

'^^"Vw^' ""^^y VT~y' ^c-^++c^- (133) 

with the cyclotron frequency ujc — qB/m and the characteristic trap dimension r. The sum 
c<j_ + c<j_|_ = qB/m is directly proportional to the charge-to- mass ratio. Inside the trap, the 
ions are excited by an azimuthal quadrupolar RF field with frequencies around uOf.. The 
ions are released from the trap and their kinetic energy is determined from a time-of-fiight 
measurement. In case of a resonant excitation ujc, the radial kinetic energy of the released ion 
will be maximal and its time of flight minimal. Unknown masses are then determined relative 
to calibration measurements of nuclei with precisely known mass. Relative mass uncertainties 
of order 10^^ are routinely achieved. The Penning trap setups used for mass spectrometry 
of exotic nuclei are ISOLTRAP [464,465] at CERN/ISOLDE (Switzerland), CPT [466,467] 
at ANL (US), SHIPTRAP [468-470] at GSI (Germany), JYFLTRAP [471,472] at Jyvaskyla 
(Finland), LEBIT [473-475] at NSCL (US) and TITAN [476] at TRIUMF (Canada). 

At GSI, rare-isotopes are produced in-flight, selected in the Fragment Separator (FRS) and 
subsequently injected into the Experimental Storage Ring (ESR) [477]. For two ion species 
circulating in the ring the relative difference in the mass-to-charge ratio m/qoi the revolving 
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ion species is expressed to first order as: 



^ = _J^«^fl_2!)^, (134) 

/ It rn/q \ -ff J u 

where /, m/q,v and 7 are the mean values of frequency, mass-to-charge ratio, velocity and 
Lorentz factor, respectively, and Af = /i — /2, A{m/q) = {mi/qi) — {m2/q2) and Av = V1—V2 
are the corresponding differences for the two ion species. The transition point 7^ characterizes 
the point where the revolution frequency becomes independent of the energy for each species 
with fixed m/q. This quantity can be varied within certain limits by adjustments of the 
ring's ion optics. The ESR can be operated in two different optics modes to deduce mass- 
to-charge ratios. For longer-lived nuclei (T1/2 of the order of seconds), the velocity spread 
Av/v of the ions can be reduced by electron coohng (Schottky method). Then the term of 
eq. (134) that depends on the velocity spread becomes negligible and the m/q can be directly 
determined from the frequency of revolution. The frequency is determined via the Schottky 
noise pickup technique. Many species can be stored in the ring simultaneously, with ions of 
precisely known mass serving as calibrants. Relative mass uncertainties are typically of order 
10~^ [478,477]. For short-lived nuclei, the ring can be operated in isochronus mode. Here, the 
Lorentz factor is chosen to match 7^, where the revolution time becomes independent of the 
velocity and the mass-to-charge ratio can be deduced. Typically, thin-foil timing detectors 
are used to measure the time of flight. Since this approach does not require cooling, short- 
lived ions (Ti/2 of order 10 /is) can be studied [479,477]. Future measurements at the ESR 
will probably also address the masses of neutron-rich nuclei, in particular those involved 
in the r process of stellar nucleosynthesis, as these can be produced efficiently by uranium 
projectile fission with subsequent separation in-flight by the FRS at GSI. 



Another frequency-based mass measurement concept is implemented at ISOLDE with the ra- 
diofrequency transmission mass spectrometer MISTRAL [480,481], where in a homogeneous 
magnetic field a controlled and coherent manipulation of the ion trajectory is performed to 
determine the mass via the cyclotron frequency. 



Time-of-flight measurements in combination with magnetic spectrometers or cyclotrons are 
used in various experimental schemes to determine the mass-to- charge ratio of short-lived 
ions. An ion's motion in a magnetic field can be characterized by the magnetic rigidity Bp 
which is connected to the mass-to-charge ratio via the ion's velocity v or cyclotron frequency 

UJc' 

-fmv vB 

Bp= = — , (135) 

q 

where 7 is the Lorentz factor. The magnetic rigidity can be measured with high-resolution 
magnetic spectrometers, the velocity can be determined from the acceleration potential at 
low energy [482] (ISOLDE) or directly from time-of-fiight measurements at higher beam 
energies [483,484] at GANIL and NSCL, for example. For a given detector time resolution, the 
resolution of the time-of-fiight measurement is limited by the total fiight time. An increased 
fiight time can be realized with the use of cyclotrons [485]. For two ions with masses m and 
m + dm, accelerated simultaneously in a cyclotron, the time difference ST, after ut turns is 
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given to first order by 5T /T = 5m/ m. Relative to a well-known mass m, the mass m + 5m 
can be derived from the flight times of the ions. 

Q-values of radioactive decays can provide accurate mass differences between parent and 
daughter nuclei. For a long time, P decay was one of the main sources of mass determinations 
{f3^ decay and electron capture (EC) for proton-rich nuclei and f3~ decay on the neutron-rich 
side of the nuclear chart). If the decay proceeds to excited states in the daughter nucleus, the 
mass difference has to be adjusted for the excitation energy, meaning that in addition to the 
determination of the (5 end-point energy, spectroscopy of the daughter has to be performed 
in coincidence (see [486] for a review). Qa or Qp measurements are less complicated owing to 
the fact that the a particle or proton do not share their energy with a neutrino, thus making 
the energy difference a more direct observable for ground-state to ground-state decays [487]. 
Mass differences may also be derived from the energy balance in reactions. In a two-body 
reaction, A{B, C)D, the mass excesses A are related via the amount of released energy Q: 

Q = Aa + Ab-Ac-Ad (136) 

If three of the masses are well known, the determination of the Q- value from a measurement of 
the reaction kinematics will allow the mass of the remaining reaction partner to be extracted. 
This methods also works for extracting the mass of unbound systems when the decay particle 
is detected (invariant mass method) [488]. For both approaches based on decay or reaction Q 
values, the mass differences that enter the energy balance must be linked to known masses. 

Measured masses are evaluated in the 2003 Atomic Mass Evaluation (2003AME) [489,490]. 
Below, we summarize some of the more recent mass measurements of exotic nuclei with 
relevance to nuclear astrophysics. 

The rapid neutron-capture process (r process) is responsible for the synthesis of roughly 
half of the elements beyond iron. However, its site and reaction path is not known with 
certainty. The abundance patterns of various r-process models show a significant dependence 
on the underlying nuclear physics (see [11,3,491,12,9,492] for reviews). High precision mass 
measurements of 7im,72-8i2n were performed with ISOLTRAP at CERN [493] . Prom ^^Zn and 
^°Zn, the neutron separation energy and neutron-capture Q-value of ^°Zn were determined 
for the first time. Depending on the stellar environment, the r-process path either includes 
the slow (3 decay of ^''Zn making it a major waiting point or it proceeds rapidly to ®^Zn and 
beyond via neutron capture [493]. The new results improved the mass-related uncertainties 
for r-process calculations compared to the previously used 2003AME. In the temperature 
regime below 1.5 GK, where the (^,7) ^ (7,^) equilibrium begins to break down, the 
previously dominating mass uncertainties to the reaction flow in this region have become 
negligible with the new results. High-precision Penning trap mass measurements of i32,i34gj^ 
carried out with the ISOLTRAP setup revealed a 0.5 MeV deviation of the binding energy 
of ^^^Sn compared to the previously accepted value [494]. The new deduced value for the 
N = 82 shell gap in ^^^Sn was found to be larger than the = 28 shell gap of the stable 
doubly magic nucleus ^^Ca. The A^ = 82 shell gap is thought to influence the flssion cycling 
in the r process [495] since, in the presence of a shell gap, the r process slows down and more 
neutrons are created by photodisintegration inducing more flssion [494] . 
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X-ray bursts are initiated when the temperature and the density in the accreted layer on a 
neutron star become high enough to allow a breakout from the hot CNO cycle. The nuclei 
start to capture protons and proceed along a capture chain via the rp process. The capture 
rate is favorable compared to the /?-decay rate of the nuclei involved until a nucleus with a 
small Q value for the proton-capture reaction is encountered. There, an equilibrium between 
the 7) proton capture and (7, p) photodisintegration develops and the rp process is delayed 
until the subsequent /9-decay of this so-called "waiting-point" nucleus [16,6]. The effective 
lifetime of a waiting-point nucleus depends strongly on its proton separation energy which 
can be directly determined from mass differences. The effective lifetime of a nucleus in an 
X-ray burst is exponentially dependent on the proton separation energy and thus accurate 
mass values - preferably with less than 10 keV uncertainty - are desired for X-ray burst 
model calculations [496,497]. The selfconjugate nuclei ^^Ge, ^^Se and ^^Kr are believed to be 
important, long-lived waiting point nuclei [6]. 

Mass measurements of ^^Se, ^^As and ^^Ge with the Canadian Penning Trap (CPT) confirmed 
^^Se to be a waiting point nucleus, causing a considerable delay in the rp process [498]. 
Penning trap mass spectrometry of ^^'^^Ga, ®^'^^'^^Ge, 66,67,68^g ^^^^ 69gg performed at 
the LEBIT facility at NSCL [499]. By using theoretical Coulomb-shift energies and measured 
masses of ^^'^^Ge, ^''As, and ^^Se, mass values for ^^As, ^^'^'^Se, and ^^Br were predicted with 
an uncertainty of about 100 keV [499]. The results together with the measured masses were 
used to calculate improved effective lifetimes of the rp-process waiting-point N — Z nuclei 
^^Ge and ^^Sc. It was found that ^^Ge is less of a waiting point while ^'^Sc poses a larger delay 
in the rp process than previously thought [499]. Recently, the masses of ®^Se, ^°Se, ™"*Br 
and ^^Br were measured with experimental uncertainties of 0.15-15 keV with the LEBIT 
facility at NSCL [500] . The new and improved data were used in conjunction with Coulomb 
displacement energies as input for rp-process calculations. An increase of the effective lifetime 
of the waiting point nucleus ^^Se was found and more precise information on the luminosity 
of type I X-ray bursts and final element abundances after the burst were obtained [500]. At 
ISOLTRAP, the mass of the selfconjugate nucleus ^^Kr was measured for the first time with 
Penning trap mass spectrometry yielding a precision of 8 keV [501]. The masses of '''^-'''^Kr 
were determined with an order of magnitude improvement in accuracy compared to previous 
results. ^^Kr was found to be an important waiting point in the rp process during X-ray 
bursts, delaying energy generation with at least 80% of its /3+-decay half-life [501]. 

Of interest is the possible existence of an endpoint to the rp process. Reaction network 
calculations suggest that the rp process in X-ray bursts terminates at tellurium and is lim- 
ited by the SnSbTe cycle [502]. The nuclei io4-io8g^^ loe-iiogb, los.iog^e, and ^^^I, in the 
vicinity of the expected endpoint of the astrophysical rp process, were produced in fusion- 
evaporation reactions with a ^^Ni beam irradiating a natural Ni target at the IGISOL facility 
in Jyvaskyla [503]. Their mass values were precisely measured at the JYFLTRAP Penning 
trap facility. For io6,io8,iiog|-,^ these were the first direct mass measurements. One-proton 
separation energies were derived and the low value for ^'^^Sb excludes a strong SnSbTe cy- 
cle starting at ^°^Sn [503]. Masses of 34 neutron-deficient nuclei along the rp-process path 
were measured with SHIPTRAP at GSI with an uncertainty of 10 keV [470]. The nuclei 
relevant for the rp-process path above A = 90 were produced in fusion-evaporation reactions 
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induced by ^°Ca, ^''Cr and ^^Ni beams impinging upon a ^^Ni target. Preliminary results 
were reported for 101,102,103,104^^^ 102,103,104,105 i05,i06gj^^ 99,101,103^5^ 112,111,110,109^6, 
113,112, iiij^ ii3^g^ i07,i09,iiig|-, Except for ^"^In and several tellurium isotopes, the measured 
masses were found to agree well with the 2003AME. The rp process proceeds along the Sn 
isotopic chain until the proton separation energy of Sb isotopes becomes large enough to 
delay the competing photodisintegration and hence, causes the rp process to proceed up to 
tellurium where SnSbTe cycles prevent the production of heavier elements. With the masses 
measured at SHIPTRAP, more accurate proton separation energies can be deduced in order 
to clarify which SnSbTe cycle dominates [470]. The Q values for the proton emission from 
^°^Sb and ^°^Sb were also recently deduced indirectly from a-decay measurements carried 
out at HRIBF (ORNL) [504]. 

The up process occurs in explosive environments when proton-rich matter is ejected in the 
presence of strong neutrino fluxes [505]. This includes the inner ejecta of core-collapse su- 
pernova [506], for example. Neutron-deficient nuclei from yttrium to palladium of relevance 
to the up process were produced in fusion-evaporation reactions and studied at SHIPTRAP 
(GSI) and JYFLTRAP (Jyvaskyla) [507]. Precise masses of 21 nuclei were determined with 
Penning trap mass spectrometry, the masses of ^^Tc, ^''"^^Ru, ^^~^"^Rh, and ^'^'^^Pd for the 
first time. The data was used for i/p-process nucleosynthesis calculations to probe their im- 
pact on reaction flow and final abundances. In particular the reaction flow around ^^Tc was 
found to be significantly modified as the new measurements give a proton separation energy 
that differs from the previously used 2003AME by 1 MeV. However, the final abundances 
were basically unchanged [507]. 

Direct mass measurements of the short-lived nuclei ^^V, ^^Mn, ^-"^Ti and '^'^Cr were performed 
with the isochronous mass spectrometry at the experimental storage ring (ESR) at GSI [508]. 
An accuracy of about 100-500 kcV was achieved. The results largely confirmed the previously 
employed theoretical mass predictions and thus lead only to small changes in the rp process 
in astrophysical X-ray burst models [508] . 

The mass and electron capture Q value of ^^Si was measured with JYFLTRAP Penning trap 
facility and a ten times improved precision was achieved [509] . This leads to an improvement 
of the reaction Q value for the ^^Al(p, 7)^^81 proton capture. The new result changed the 
stellar production rate of ^^Si at nova ignition temperatures by about 10% and makes the 
limiting factor for higher precision of the reaction Q value the mass of ^^Al [509] . 



7 Outlook 

Tremendous progress has been made over the past years in experimental nuclear astrophysics 
as well as in the developments of the multi-faceted theoretical framework that is crucial 
to link the study of nuclear structure and reactions to astrophysical processes in stars and 
explosive scenarios. The future for nuclear physics and nuclear astrophysics looks bright with 
next-generation rare-isotope facilities on the horizon like GSI/FAIR in Germany and FRIB 
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in the United States and major upgrades on the way at GANIL in Prance and TRIUMF in 
Canada. 

Short-term, the RIBF facihty in JAPAN, which became operational recently, and significant 
capability upgrades at NSCL, TRIUMF, ANL, HRIBF/ORNL, and CERN/ISOLDE (see 
Section 6.1) will continue to provide a multitude of exciting opportunities to advance the 
experimental research in nuclear astrophysics at all fronts. For example, the nuclear equation 
of state (EOS) - crucial for an understanding of neutron stars and supernova explosions - can 
be probed with fast-beam heavy-ion collisions induced by projectiles with extreme N /Z ratios 
utihzing the present NSCL fast-beam capability as well as the different energies available at 
GSI and RIBF. Pioneering programs at NSCL aimed at charge-exchange reactions on exotic 
nuclei will measure electron-capture rates for exotic nuclei and generate crucial input for the 
nucleosynthesis in supernovae. The CARIBU at the ATLAS facility at ANL will provide low- 
energy and reaccelerated beams of neutron-rich ^ssQf fission products and enable experiments 
in the proximity of r- process path in several mass regions [275]. The ReA3 reaccelerator 
addition underway at NSCL will deliver for the first time reaccelerated beams of rare isotopes 
produced by projectile fragmentation and fission with final energies ranging from 0.3 to 
3 MeV/nucleon for ions with a charge-to-mass ratio of 0.25, and 0.3 to 6 MeV/nucleon for 
ions with a ratio of charge to A of 0.5, which are of great interest for direct and indirect 
nuclear astrophysics measurements across the chart [274,446]. 

The important scientific questions to be addressed both experimentally and theoretically in 
nuclear physics of exotic nuclei with relevance for astrophysics comprise: (a) How do loosely- 
bound systems survive and what are the general laws of their formation and destruction? 
(b) Are new types of radioactivity possible? (c) Are new types of nuclear symmetry and 
spatial arrangement possible? (d) What are the limits of nuclear existence? (e) How do the 
properties of nuclear matter change as a function of density, temperature and proton-to- 
neutron ratio? (f) How do thermal and quantum phase transitions occur in small systems? 
(g) What determines the shape and symmetry properties of an exotic nucleus? (h) How does 
quantum tunneling of composite particles occur in the process of reactions and decay? (i) 
What are the manifestations of fundamental forces and symmetries in unusual conditions? 
(j) How were the elements heavier than iron formed in stellar explosions? (k) How do rare 
isotopes shape stellar explosions? (1) What is the role of rare isotopes in neutron stars? 

These questions provide extreme challenges for experiments and theory. On the experimental 
side, producing the beams of radioactive nuclei needed to address the scientific questions has 
been an enormous challenge. Pioneering experiments have established the techniques and 
present-generation facilities have produced first exciting science results, but the field is still 
at the beginning of an era of discovery and exploration that will be fully underway once the 
range of next-generation facilities will become operational. The theoretical challenges relate 
to wide variations in nuclear composition and rearrangements of the bound and continuum 
structure, sometimes involving near-degeneracy of the bound and continuum states [510]. 
The extraction of reliable information from experiments requires a solid understanding of 
the reaction process, in addition to the structure of the nucleus. In astrophysics, new observa- 
tions, for example the expected onset of data on stellar abundances, will require rare-isotope 
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science for their interpretation. 
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